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INTRODUCTION 


)f 

Under  these  contracts,  MPL  was  engaged  in  a  wide  variety  of  Navy  related 
activities  which  are  grouped/categorized  as  follows: , 


1.  Augmented  Exploratory  Research  and  Development:  Expansion  or  acceleration 
of  ongoing  effoi  .s,  principally  in  underwater  acoustics,  signal  processing, 
hydrophone  arrays.  Also,  in  optics  with  the  folding  of  Visibility  Laboratory, 
activities  within  MPL,  enhanced  Deep  Tow  capabilities, 

2.  Support  of  Research  and  Development:  FLIP/ORB  support  and  stable  platform 
related  issues, 

3.  Support  of  MPL  Sponsor  related  needs: 

Operating  CNR  tower.  Equipping  of  USS  DOLPHIN  with  research  gear.  Surface 
Effect  Ship  (SES),  SUBDEVGRPONE,  modifications  of  existing  equipment  or 
equipment  orders  for  MPL  or  MPL  sponsor  related  needs  (:.e ,  NAVOCEANO, 
NAVAIR). 
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4.  Special  Projects:  Quick  response  MPL  support  of  NAVY  mission  related  needs. 
This  w'juid  be  utilizing  MPL  assets  such  as  Deep  Tow  to  meet  unforseen 
requirements. 

5.  Exploratory  and  Advisory:  Directly  funded  participation  of  MPL  scientific, 

engii leering,  and  other  personnel  in  various  Na-y  panels  and  committees,  etc. 
This  is  for  participation  beyond  the  normal  Exploratory  and  Advisory  in 
connection  with  the  overall  MPL  program. _ 


These  activities  are  within  the  framework  of  our  Navy  mission  which  is: 

To  investigate  and  apply  knowledge  about  the  ocean,  its  boundaries 
and  their  surrounding  media  to  the  solution  of  the  Navy’s  problems  in 
anti  submarine  and  pro-submarine  warfare; 

To  provide  research  training  of  students  in  areas  of  oceanography 
and  ocean  technology  which  have  application  to  Navy  requirements. 

To  maintain  certain  special  ocean  engineering,  research  and  develop¬ 
ment  capabilities  which  are  essential  to  the  Navy. 


Categories  1  and  2  are  the  most  important  from  the  standpoint  of  developing 
the  research  skills  and  facilities  with  which  to  pursue  research  f.oals  within  the 
context  of  the  University  and  our  Navy  mission.  Capabilities  developed  as  a  result 
of  the  above  have  led  to  many  instances  of  being  asked  by  the  Navy  for  quick 
response  efforts  which  are  in  Categories  3  and  5.  Similarly,  the  skills  and 
experience  acquired  as  a  result  of  activities  in  Categories  1  and  2  led  to  requests  by 
the  Navy  to  serve  on  i  inels  and  committees  on  a  sustained,  funded  basis;  that  is 
Category  5,  which  is  aoove  and  beyond  our  usual  Exploratory  and  Advisory 
activities. 
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SUMMARY  OF  ACTIVITIES 


A  brief  description  of  these  activities  is  included  below  along  with  a  detailed 
program  summary  listing  by  categories.  A  bibliography  of  reports  regarding  these 
activities  is  also  included  in  Appendix  A: 


Category  1 

It  is  appropriate  to  summarize  activities  where  there  is  a  common  thread. 

One  of  MPL’s  major  research  efforts  directly  related  to  our  mission  is  in  the 
area  of  understanding  and  exploiting  ambient  noise  and  sound  propagation  in  the 
ocean.  The  work  of  Hodgkiss,  Hildebrand  and  Fisher  was  concerned  with  the 
development  and  deployment  of  hydrophone  arrays,  and  water  column 
measurements  of  ambient  noise  and  signals  in  various  experiments.  This  included 
joint  experiments  with  various  Navy  Centers,  NRL  and  NORDA,  as  well  as 
contractors.  Dr.  Hodgkiss’s  v/ork  with  signal  processing  and  array  processing  is 
directly  related  to  this  area.  At  lower  frequencies,  the  work  of  Drs.  Anderson,  Webb, 
Dorman,  Hildebrand  and  Hodgkiss  made  use  of  different  sensors  that  are  bottom 
mounted  to  itudy  ambient  noise  and  sound  propagation. 

For  acoustic  measurements  in  the  water  column,  two  major  experiments  were 
successfully  completed  in  which  unique  new  measurements  of  the  vertical 
directionality  of  the  ambient  noise  field  were  made  from  FLIP:  one,  a  study  of  the 
effect  of  distance  from  coastal  shipping  and  local  wind  on  the  vertical  distribution 
of  ambient  noise  in  the  75-300  Hz  region  (Broadband  CUARP);  and,  two,  the  first 
high  resolution  measurements  of  the  vertical  directionality  of  signal  and  noise  in 
the  50-200  Hz  region  (ONT  High  Gain  Initiative  SVLA  experiment).  A  third  major 
effort  in  this  area  was  preparation  for  the  DVLA  experiment  which  was  a  follow-on 
to  the  SVLA  experiment.  Whereas  the  SVLA  experiment  made  use  of  a  900  meter 
vertical  array  originally  designed  for  horizontal  deployment,  its  center  frequency 
was  100  Hz,  too  high  for  the  goals  of  the  DVLA  experiment  scheduled  for  July, 

1989.  Construction  of  a  3000  meter,  200  element  (50  Hz)  vertical  array  was 
completed  in  time  for  the  DVLA  experiment  and  it  was  successfully  operated  for  the 
duration  of  the  experiment.  The  9G0  and  3000  meter  arrays  produce  12  bit  digital 
data  from  each  hydrophone  and  make  use  of  acoustic  navigation  for  array  element 
location  to  make  beamforming  and  matched  field  processing  possible. 

In  conjunction  with  the  above  efforts,  freely  drifting  Swallow  float  arrays  with 
3  axis  geophones  were  deployed  by  Dr.  Hodgkiss  to  obtain  ambient  noise  and 
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signal  data  uncontaminated  by  flow  noise  to  compare  with  data  obtained  with 
arrays  deployed  from  FLIP.  The  Swallow  floats  have  also  been  the  key  assets  in 
pursuing  Slack  Line  array  concepts,  sea-floor  measurements,  and  led  to  the  need  for 
an  acoustic  modem  so  that  data  at  a  high  rate  (9600  baud)  can  be  transmitted  to  a 
surface  vessel  from  various  Swallow  float  arrays.  These  Swallow  floats  can  be 
deployed  in  vertical  or  horizontal  configurations  and  have  an  array  element 
location  capability  built  into  them  so  data  from  them  can  be  used  in  post-cruise 
beamforming. 

Array  research  with  the  MPL  200  meter  array  was  concerned  with  the  vertical 
distribution  of  ambient  noise  in  the  100-400  Hz  region  and  received  support  from 
NAVAIR  in  conjunction  with  the  Broadband  CUARP. 

in  all  the  above  research,  data  processing  and  analysis  are  a  key  to  monitoring 
experiments  in  progress  and  interpreting  the  results.  MPL  has  been  involved  with 
developing  the  computer  and  recording  systems  for  our  own  research  and,  in  so 
doing,  we  have  rlso  been  tasked  to  build  high  speed  data  acquisition  systems  for 
other  uses  by  our  sponsors,  notably  the  ACSAS  system,  as  well  as  applying  the 
dynamic  beamformer  (built  for  the  STRAPP  experiment)  to  towed  arrays. 

Doppler  sonars  mounted  on  FLIP  were  used  by  R.  Pinkel  to  study  upper  ocean 
dynamics  in  about  a  cubic  kilometer.  These  sonars  have  been  used  to  measure 
currents  down  to  1  cm/sec  relative  to  FLIP  and  have  revealed  Langmuir  cells,  and 
surface  slope  distributions  in  addition  to  delineating  internal  waves  and  directions 
of  energy  flux. 

Multispectral  imaging  of  the  marginal  ice  zone  and  whole  sky  imaging  with 
respect  to  cloud  cover  (E/O  System  5  units)  are  parts  of  an  overall  program  to  build 
a  database  from  which  R.  Johnson  and  his  colleagues  can  ''mplement  autonomous 
pattern  recognition  algorithms  for  quantitative  remote  sensing  of  cloud  cover  and 
visibility.1 


^  Upon  dissolution  of  the  Visibility  Laboratory  of  the  Scripps  Institution  of 
Oceanography,  a  part  of  the  "Whole  Sky  Imagery"  program  (Richard  johnson/Principal 
Investigator)  came  under  the  Marine  Physical  Laboratory. 
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Category  2 


Facilities  support  for  FLIP  and  ORB  (6.5)  from  NAVSEA  05R12  includes  support 
for  personnel  crewing  these  platforms  and  regular  maintenance  and  inspection  such 
as  periodic  drydocking.  For  example,  with  respect  to  FLIP,  whereas  it  was  originally 
designed  for  a  twenty-year  life  based  on  marine  corrosion  allowances,  excellent 
corrosion  protection  has  allowed  us  to  continue  operating  FLIP  with  a  more  intense 
maintenance  program  v/ith  respect  to  monitoring  fatigue  cracking.  Strain  gage 
recordings  during  sea  operations  and  towing  provide  a  quantitative  measurement  of 
stress  cycling.  After  the  equivalent  of  a  normal  years  exposure  to  cyclic  stresses, 
FLIP  is  drydocked  for  inspection  and  maintenance.  These  platforms  have  been  used 
in  a  wide  variety  of  tasks,  mostly  generated  within  MPL;  several  outside  users  have 
also  used  these  platforms  on  mostly  Navy  related  programs. 

Results  of  a  workshop  conducted  by  MPL  led  to  studies  of  advanced  research 
platforms,  including  a  larger,  more  capable  version  of  FLIP  and  the  conversion  of  an 
existing  large  offshore  oil  exploration  semi-submersible  platform  for  multi¬ 
disciplinary  research  programs.  The  new  FLIP  has  been  proposed  for  FY-92  support, 
and  the  semisubmersible  has  been  acquired  by  the  Department  of  Defense  for  use 
in  the  drug  enforcement  program,  with  opportunities  for  research  support. 


Category  3 

In  this  category  is  included  a  wide  range  of  activities  in  which  MPL  has 
participated  directly  in  supporing  Navy  needs  or  has  acted  as  cognizant  monitor  of 
various  Navy  projects  within  Navy  activities  or  at  contractor  companies. 

Several  tasks  were  assigned  to  MPL  under  this  contract  in  order  to  make 
special  talents  and  capabilities  at  MPL  available  to  Navy  sponsors,  particularly  with 
regard  to  advanced  capability  platforms.  These  included  the  successful 
rehabilitation  and  operation  of  an  offshore  oceanographic  research  tower  for  the 
Chief  of  Naval  Research,  until  it  was  destroyed  in  a  major  storm  in  January  1988. 
The  research  submarine  USS  DOLPHIN  (AGSS-555)  was  equipped  with  two  major 
research  systems,  one  involving  an  off-board  acoustic  and  optical  sensor  system 
and  the  other  involving  the  installation  of  a  precision  photo-optical  system  for 
seafloor  investigations.  Also  included  in  this  category  wrere  a  number  of  tasks 
devoted  to  increasing  the  capability  of  Navy  deep  submergence  assets  to  support 
research  and  survey  programs,  both  at  Submarine  Development  Group  ONE  and  at 
the  Naval  Oceanographic  Command.  Investigations  of  advanced  capability 
platforms  included  studies  of  the  producibility  of  high-technology  surface  effect 
ships  under  contracts  with  several  major  shipyards  (SES  Hullform),  and  tests  and 
evaluation  of  a  new  closed-cycle-diesel  propelled  submarine  for  potential  diver 
lock-out  capability. 
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Category  4:  Quick  Response 

The  loss  of  a  high  priority  space  vehicle  initiated  an  urgent  Navy  request  for 
MPL  to  conduct  a  prompt  and  detailed  search  of  a  seafloor  area  in  the  East  Pacific 
Ocean.  This  was  accomplished  using  SEABEAM  sonar  survey  data  to  delineate  the 
most  probable  search  areas  for  MPL’s  Deep  Tow  instrument  system,  and  for  the 
coordinated  search  by  Navy  deep  submergence  assets. 

Several  times  in  the  past,  urgent  Navy  needs  have  required  the  services  of  MPL 
in  operations  at  sea.  During  the  pe.  iod  of  this  contract,  F.  Spiess  and  P.  Lonsdale 
participated  in  sea-floor  search  activities  with  the  MPL  Deep  Tow  system. 


Category  5:  Exploratory  and  Advisory 

Normally,  MPL  principal  investigators  participate  in  various  workshops  and 
advisory  panels  as  part  of  their  regular  activities.  For  sustained  efforts  along  these 
lines,  particularly  when  a  lot  of  travel  is  involved,  funding  in  this  category  is 
provided. 
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CONTRACT  PROGRAM  SUMMARY 


TABLE  I 

Category  1:  Augmented  Fxploratory  Research  and  Development 
(Array  Deployments,  Deep  Tow  Research) 


Date  of 
Modification 

Sponsor 

Program 

Title 

Principal 

Investigator 

Mod  00 

NAVSEA  63D 

Appl.  of  Dynamic  Beamf.  to  Towed  Line 
Arrays 

Hodgkiss 

Mod  01 

NAVSEA63R 

(VIM)  Vibration  Isolo'.ion  Module 

Andrews 

Mod  03 

NAVSEA63R14 

Adaptive  Signal  Processing 

Hodgkiss 

Mod  04 

NAVSEA63R 

Reverberation  Dynamics 

Hodgkiss 

Mod  04 

NAVSEA63R 

(VIM)  Vibration  Isolation  Module 

Andrews 

Mod  05 

NAVSEA63D1 

Appl.  Dynamic  Beamforming  to  Towed 

Line  Arrays 

Hodgkiss 

Mod  05 

NAVSEA63R 

Vibration  Isolation  Module 

Andrews 

Mod  05 

NOSC  7213 

Active  Acoustic  Cancellation 

Abarbanel 

Mod  07 

NAVSEA63R14 

Adaptive  Array  Processing 

Hodgkiss 

Mod  07 

NAVSEA63R12 

Devel.  &  Evel.  Vibration  Isolation  Module 

Andrews 

Mod  08 

NAVSEA63R12 

Devel.  &  Eval.  Vibration  Isolation  Module 

Andrews 

Mod  13 

NOSC743 

Intelligent  Beamformer 

Hodgkis 

Mod  13 

NOSC512 

Range-extended  Sonar  Targets 

Alter 

Mod  13 

NOSC541 

SEA  SNAKE  Surveillance  System 

Hodgkiss 

Mod  14 

PMS395 

Modularized  Towed  Deep  Sea  Vehicle 

Spiess 

Mod  16 

ONR112 

Marginal  Ice  Zone  Investigations 

Maynard 

Mod  19 

NAVAIR 

Ambient  Noise  Vert.  Directionality 

Hodgkiss 

Mod  19 

ONR112 

Spectral-Properties  at  Marginal  Ice  Zones 

Maynard 

Mod  19 

NOSC522 

Plankto’iic  Growth  cn  Ph  Seawater 

Dickson 

Mod  19 

NAVA1R633 

Vertical  Correlation  Measurements 

Fishes 

Mod  23 

ONR112 

Spec.  Properties  Marginal  Ice  Zones 

Maynard 

Mod  28 

NAVA1R633 

Vertical  Correlation  Measurements 

Fisher 

Mod  30 

ONR/AF 

Multi-Station  WSI  Stations 

Johnson 

Mod  30 

NOSC 

SEA  SNAKE  Surveillance 

Hodgkiss 

Mod  30 

NUSC 

VLF  Acoustic  Experiment 

Hodgkiss 

Mod  31 

ONT 

High  Gain,  High  Resolution  Program 

Fisher/Hildebrand 

Mod  32 

ONR/AF 

Multi-Station  WSI  Network 

Johnson 

Mod  32 

NOSC 

SEA  SNAKE  Surveillance 

Hodgkiss 

Mod  32 

NOSC 

Reverberation  Simulation 

Hodgkiss 

Mod  32 

ONR/AF 

Multi-Station  WSI  Station 

Johnson 

Mod  32 

NOSC 

SEA  SNAKE  Surveillance 

Hodgkiss 

Mod  33 

NUSC 

Parametric  Sonar  Support 

Hodgkiss 

Mod  33 

NOSC634 

Effects  of  Monomolecular  Films 

Anderson 
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Marine  Physical  Laboratory 


MPL-U-03/89 


Date  of 
Modification 

Sponsor 

Program 

Title 

Principal 

Investigator 

Mod  35 

ONT/ Booth 

Vertical  Line  Array  Measurements  of 
Ambient  Noise 

Fisher/Hodgkiss 

Mod  35 

ONR 

Upper  Ocean  Effects 

Pinkel 

Mod  35 

CMT/NRL 

Program  in  Marine  Physics 

Watson 

Mod  36 

ONR/NRL 

Program  in  Marine  Physics 

Watson 

Mod  36 

ONT/Booth 

MPL  Digital  Array  Design  Studies 

Hildebrand/Fisher 

Mod  37 

ONT/Booth 

MPL  Digital  Acoustic  Array  for  the  DVLA 
Experiment 

Hildebrand/Fisher 

Mod  38 

ONR/NRL 

Program  in  Marine  Physics 

Watson 

Mod  38 

White  Sands 

E/O  System  5  Camera  Spare  Parts 

Johnson 

Mod  39 

ONR 

Upper  Ocean  Effects 

Pinkel 

Mod  39 

NRL/SPAWAR 

Certical  Sea  Test 

Fisher/Hodgkiss 

Mod  41 

ONT/NRL 

Program  in  Marine  Physics 

Watson 

Mod  41 

AF 

E/O  System  5 

Johnson 

Mod  42 

DARPA 

ACSAS  Program 

Hodgkiss 

Mod  42 

Army 

E/O  System  5 

Johnson 

Mod  43 

ONR/SDI 

E/O  System  5 

Johnson 

Mod  45 

ONT/230 

Vertical  Array-Hi  Gain 

Hildebrand/Fisher 

Mod  45 

ONT/230 

Vertical  Array  Data  Analysis 

Hodgkiss/Hiluebrand 

Mod  46 

SDI/'Army 

E/O  Sys5  Unit  #3 

Johnson 

Mod  46 

SDI/Army 

Night-Time  Prototype 

Johnson 

Mod  48 

ONT 

Array  Analysis 

Hodgkiss 
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MPL-U-03/89 


TABLE  II 

Category  2:  Support  of  Research  &  Development  -  Facilities  Support 


Date  of 

Modification  Sponsor 


Program  Principal 

Title  Investigator 


Mod  02 

NAVSEA05R 

Mod  04 

NAVSEA05R 

Mod  07 

NAVSEA05R 

Mod  10 

NNAVSEA05R12 

Mod  24 

NAVSEA05R 

Mod  35 

NAVSEA05R 

Mod  42 

NAVSEA05R 

Mod  46 

NAVSEA05R 

Mod  48 

NAVSEA05N 

R/P  FLIP  and  ORB  Bishop 

R/P  FLIP  and  ORB  Bishop 

R/P  FLIP  and  ORB  Bishop 

R/P  FLIP  &  ORB  Support  Bishop 
R/P  FLIP  and  ORB  Bishop 

R/P  FLIP  and  ORB  Bishop 

ORB  Thruster  Beck 

FLIP/ORB  Support  Beck 

FLIP/ORB  Beck 
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MPL-U-03/89 


TABLE  III 

Category  3:  Various  MPL  Sponsor  Research  Activities 


Date  of 
Modification 

Sponsor 

Program 

Title 

Principal 

Investigator 

Mod  03 

NAVOCEANO 

Crane  Support 

Bishop 

Mod  05,08,  13 

PMS395 

Techn.  Support  for  USN  Deep  Subm. 

Bishop 

Systems 

Mod  05 

NAVOCEANO 

Techn.  Support  for  USN  Deep  Subm. 

Boegeman 

Tow 

Mod  05 

NAVOCEANO 

Transmissometer  Repair  &  Maintenance 

Edwards 

Mod  05 

NOSC1603 

ACSASS  Calibration  &  Analysis  Support 

Vislab/ Austin 

Mod  06 

NAVOCEANO 

Slack  Tensions  inspection/Repair 

Bishop 

Mod  09 

NOSC943 

Sample  Recovery  Basket 

Bishop 

Mod  09 

NOSC844 

NOSC  Radiometric  Facility 

Austin 

Mod  10 

NAVAIR 

Pacific  Bathymetry 

Spiess 

Mod  13,  16,  25 

NAVSEA05R12 

Wide-Area  Imaging  System 

Bishop/Ballard 

Mod  14 

NAVO7210 

Multi-Spectral  Optical  Probe 

Austin 

Mod  14 

NAV07111 

Transmissometer  Modifications 

Austin 

Mod  19 

NAVSEA0133 

Passive  Sonar  Improvements 

Maynard 

Mod  19,  32,  35,  48 

PMS395 

Eval.  Acoustic  &  Optical  Sensor  Systems 

Bishop 

Mod  19 

NOSC641 

Support  for  USS  DOLPHIN 

Bishop 

Mod  19 

PMS395 

Seacliff  Support 

Bishop/Ballard 

Mod  26 

NAVSEA05R12 

Uplooking  Imaging  System  (Seacliff) 

Bishop/Ballard 

Mod  27 

NAVOCEAN071 1 1 

Repair  &  Modification  of 

Edwards 

Transmissometer 

Mod  29 

NAVSEA05R 

SES  Hullform  Technology 

Bishop 

Mod  30 

NAVOCEANO 

Camera  Spares-NAVOCEANO 

Johnson 

Mod  31,  32 

NWC 

NWC  Use  of  Tower 

Bishop 

Mod  33 

NAVSEA05R47 

Floating  Platform  Study 

Bishop/Fisher 

Mod  33 

NAVSEA05R47 

SES  Supplment 

Bishop 

Mod  33 

NAVAIR933 

Broadband  88  (Sonobuoys) 

Fisher 

Mod  33 

NOSC/NRL 

Undersea  Equipment  Support 

Watson 

Mod  34 

NAVOCEANO 

MERT  Support 

Austin 

Mod  34 

ONR1112D 

Tower  Support 

Bishop 
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MPL-U-03/89 


TABLE  III  (corn) 

Category  3a:  Various  MPL  Sponsor  Research  Activities 


Date  of 

Modification  Sponsor 


Program 

Title 


Principal 

Investigator 


Mod  35,  40.  48 

PMS395 

Evaluation  of  Diver-Lockout  Subm. 

Bishop 

Mod  35 

PMS395 

Eval.  of  Acoust.  &  Oftical  Sensor  Systm. 

Bishop 

Mod  36 

NAVOCEANO 

Deep  Towed  Vehicle  Cable  Tensioning 
System 

Bishop 

Mod  38 

PMS395 

Wide-Area  Wide  Area  Imaging  System 

Bishop 

Mod  38 

NAVOCEANO 

MERT  Calibration 

Nolten 

Mod  39 

NADC 

Mod  of  AN/SSQ-538  Sonobuoys 

Fisher 

Mod  40 

ONR1112D 

Tower  Add  On-Saxon 

Watson 

Mod  40 

PMS395 

Wide-Area  Wide  Imaging  System 

Watson 

Mod  43 

NOSC 

FLIP  Use 

Beck 

Mod  44 

ONR/PMS395 

AOS-Hydroproducts 

Bishop 

Mod  44 

DARPA 

CM/CCM 

Anderson 

Mod  45 

ONR/112D 

Saxon  Support 

Beck 

Mod  46 

NAVSEA/DTRO 

SCAT  Prototype  Comp. 

Beck 

Mod  48 

DARPA 

CM/CCM 

Anderson 

Mod  48 

NADC 

Mod  of  AN/SSQ  Sonobuoy 

Fisher 

Mod  48 

DTRC 

SCAT  Prototype  Com. 

Beck 

Mod  48 

PMS385 

AOS  Evaluation 

Bishop 

Mod  48 

NEL 

AOS  Evaluation 

Bishop 

Mod  49 

NADC 

Mod.  of  AN/SSQ  Sonobuoys 

Fisher 

Mod  49 

DARPA 

SCAT 

Beck 

Mod  50 

NAVSEA 

SCAT 

Beck 

Mod  51 

NAVOCEANO 

MERT  Calibration/Software 

Nolten 

Mod  51 

NAVOCEANO 

Diver  Lockout 

Bishop 
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MPL-U-03/89 


TABLE  IV 

Category  4:  MPL  Quick  Response  Special  Projects 


Date  of 
Modification 

Sponsor 

Program 

Title 

Principal 

Investigator 

Mod  10 

NAVSEAOOc 

Seafloor  Search  &  Survey 

Lonsdale/Spiess 

Mod  13 

NAVSEAOOc 

Seafloor  Search  &  Survey 

Lonsdale/Spiess 

Mod  15 

NAVSEAOOc 

Seafloor  Search  &  Survey 

Lonsdale/Spiess 

Mod  17 

NAVSEAOOC 

Seafloor  Search  &  Survey 

Lonsdale/Spiess 

Mod  18 

NAVSEAOOc 

Seafloor  Search  &  Survey 

Lonsdale/Spiess 

TABLE  V 

Category  S:  Exploratory  &  Advisory 


Date  of 
Modification 

Sponsor 

Program 

Title 

Principal 

Investigator 

Mod  04 

NAVSEA63R 

MOST  Support 

Anderson 

Mod  04 

NOSC 

Design  Study  -  Night-time  K  System 

Visla’u/ Austin 

Mod  05 

NAVSEA63R 

Expl.  &  Dev.  Undersea  Warefare 
Techn.  Contrib.  Study 

Andrews 

Mod  33 

DARPA 

Mid-Frequency  Array  Study 

Anderson 

Mod  36 

NOSC 

FDS  Review  Panel 

Hodgkiss 

Mod  38 

DARPA 

Very  Large  Mid-Frequency  Array 

Anderson 

Mod  51 

NORDA 

Deep  Tow  Consult/Shop  Services 

Spiess 
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Abstract — The  problem  of  beam  formation  from  a  towed  line  array 
whose  shape  has  been  distorted  is  considered.  Emphasis  is  placed  on 
the  beam  broadening  and  range  estimation  effects  of  array  shape 
perturbations  and  how  the  resulting  losses  can  be  regained  if  the 
actual  element  positions  are  known.  Specific  illustrations  are  provided 
for  various  levels  of  shape  distortion.  For  example,  a  15-m  bow  in  a 
232.5-m-long  array  broadens  tke  beamwidth  by  a  factor  of  3  at  50  Hz. 
As  another  example,  a  6-m  bow  in  an  800-m-long  array  leads  to  a 
20-percent  range  underestimation  at  10  km  for  a  100-500-Hz  broad¬ 
band  source. 


I.  INTRODUCTION 

rT,OWED  LINE  ARRAYS  provide  a  means  of  physically 
separating  acoustic  sensors  from  the  relatively  high  noise 
environment  of  the  towing  platform  and  the  aperture  size 
limitations  imposed  by  the  length  of  the  platform.  Against 
these  advantages  are  the  disadvantages  inherent  to  an  array 
processing  situation  where  the  sensor  locations  vary  dynami¬ 
cally  due  to  various  hydrodynamic  forces  and  towing  platform 
maneuvers. 

The  subject  of  this  paper  is  the  problem  of  beam  formation 
from  a  towed  line  array  whose  shape  has  been  distorted. 
In  Section  II,  plane  wave  beamfciming  will  be  considered. 
Correspondingly,  Section  III  will  discuss  curved  wavefront 
or  focused  beamforming.  In  both  sections,  the  emphasis  is 
on  the  effects  of  array  shape  perturbation  and  how  the  re¬ 
sulting  losses  can  be  regained  if  the  actual  array  element  posi¬ 
tions  are  known.  Towed  arrays  have  been  the  subject  of  con¬ 
siderable  attention  over  the  past  decade,  although  much 
of  the  literature  is  not  readily  accessible.  Articles  which  can 
be  found  in  the  open  literature  in  the  areas  of  beamforming 
•and  both  bearing  and  range  estimation  include  references 

[1]-[13]  • 


II.  BEAMFORMING 

The  advent  of  long  towed  arrays  has  improved  significantly 
the  ability  to  detect  and  track  distant  sources  A  line  array 
has  a  fundamental  limitation  in  that  it  exhibits  a  right- 
left  ambiguity.  One  method  of  resolving  this  ambiguity  is 
for  the  towing  platform  to  undergo  a  maneuver  thus  altering 
the  relative  spatial  orientation  of  the  source  and  platform. 

Manuscript  received  January  17.  1983;  revised  May  9,  1983.  This 
work  was  supported  by  the  Office  of  Naval  Research  under  Contract 
N00014-80-C-0220.  The  paper  is  a  contribution  of  the  Scripps  Institu¬ 
tion  of  Oceanography,  new  series. 

The  author  is  with  the  Marine  Physical  Laboratory,  Scripps  Institu¬ 
tion  of  Oceanography,  San  Diego,  CA  92152. 
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E  n  nth  element  position  vector 
S  ,  |th  source  direction  unit  vector 
B  m  mth  beam  direction  unit  vector 


Fig.  1.  Coordinate  system  definition  and  vector  notation. 


The  major  difficuky  with  this  approach  is  that  conventional 
beamformers  can  not  continue  accurate  beamforming  during 
the  maneuver  due  to  array  deformation.  A  certain  amount  of 
time  must  elapse  (from  several  minute*  to  tens  of  minutes 
depending  on  the  array  length  and  platform  speed)  for  the 
array  to  straighten  out.  Since  the  source  track  most  likely 
will  have  been  lost  during  this  period,  it  must  be  reacquired 
at  the  completion  of  the  maneuver. 

The  mathematics  of  beamforming  is  developed  below  along 
with  several  examples,  illustrating  the  influence  of  severe 
array  distortion  on  conventional  beamformer  output.  Related 
literature  includes  [  1  ] -[3] . 

A.  Plane  Wave  Beamforming 

Consider  a  coordinate  system  defined  as  in  Fig.  1.  The 
beamforming  task  consists  of  generating  the  waveform  bm(t) 
for  each  desired  steered  beam  direction  Bm.  Each  bm(t)  con¬ 
sists  of  the  sum  of  suitably  time  delayed  replicas  of  the  indi¬ 
vidual  element  signals  en(t).  The  time  delays  compensate 
for  the  assumed  differential  travel-time  differences  between 
sensors  for  a  signal  from  the  desired  beam  direction. 

Let  the  output  of  an  element  located  at  the  origin  of 
coordinate  due  to  the  /th  source  by  s^t).  Under  the  assump¬ 
tion  of  plane  wave  propagation,  a  source  from  direction  S, 
produces  the  following  sensor  outputs 

/  £' 
e«(0  =  S/(f +  — 

where  c  is  the  (constant)  speed  of  propagation  (independent 
of  location)  and  E„'  indicates  the  true  nth  element  position 
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Fig  2.  Array  geometries: 

which  may  differ  from  the  assumed  location  E„.  Appropri¬ 
ately  delaying  the  individual  element  signals  to  point  a  beam 
in  the  direction  Bm  yields  the  beamformer  output 


where  the  w„  are  weights  which  have  been  applied  to  each 
element  signal. 

B.  Examples 

Beamforming  for  the  six  array  geometries  depicted  in  Fig. 
2  will  be  used  in  examples  to  supplement  the  foregoing  mathe¬ 
matical  derivation..  Each  array  has  32  elements  with  a  constant 
interelement  spacing  of  7.5  m  thus  yielding  a  total  array 
length  of  232.5  m.  The  cases  considered  form  a  progression 
of  increasing  array  bow,  starting  with  an  exactly  straight 
array  and  continuing  in  15-m  increments  up  to  75-m  bow, 
which  corresponds  roughly  to  a  half-circle  geometry.  These 
shapes  have  been  chosen  for  the  sake  of  simplicity  and  are 
not  intended  to  represent  any  particular  operational  scenario. 
For  the  50-Hz  analysis  frequency  considered  and  an  assumed 
sound  speed  at  1500  m/s,  X  =  30  m  and  the  array  has  an  inter- 
elemeftt  spacing  of  X/4,  total  length  of  7  3/4  X,  and  bow  incre¬ 
ments  of  X/2  from  0  X  (case  1)  to  2  1/2  X  (case  6) 

The  remaining  figures  illustrate  four  source  direction  ->itua- 
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CASE  A  "  BO* 


(b) 


MCT&b 

CASl  S  60  •»*  BO* 

(d) 


»£1CAS 

CASE  6  75  m  BO* 


(f) 

elements,  7.5-m  spacing. 

tions:  0  =  90°,  <f>  =  0°  (Fig.  3);  6  =  45°,  0  =  0°  (Fig.  4):  6  = 
135°,  <P  -  0°  (Fig.  5);  and  both  6  =  45°  and  6  =  90  ,  £  - 
0°  (Fig.  6).  Each  figure  investigates  all  six  cases  of  array  bow. 
Two  bearing  response  plots  are  provided  for  each  case-  1)  as¬ 
sumed  element  positions  same  as  actual  element  positions 
(£n  =  Em')  and  2)  assumed  element  positions  always  those  of 
case  1  (straight  array).  The  two  bearing  response  plots  cor¬ 
respond  to  (3)  where  wn  =  1,  Sj(f)  =  exp  (/cur),  and  the  set 
{£  }  consists  of  128  beam  vectors  uniformly  spaced  in  spatial 
angle  from  0  =  0°  through  6  =  180°.  All  of  the  plots  have 
been  normalized  by  their  own  individual  maximum  response 
value  versus  bearing. 

Figs.  3-5  investigate  situations  involving  a  single  source. 
Ir,  all  six  cases  for  each  Figure,  the  first  bearing  response  plot 
is  that  for  a  beamformer  which  knows  the  array  element 
positions  exactly.  As  expected  for  uniform  weighting  of  the 
element  signals,  the  first  sidelobe  is  at  -13  dB  for  case  1 
(straight  array)  and  -7  dB  for  case  6  (half-circle  array).  The 
second  bearing  response  plot  is  that  for  a  beamformer  which 
assumes  the  element  positions  are  those  of  the  straight  array 
irrespective  of  their  true  locations.  The  effect  of  not  properl) 
incorporating  the  actual  array  element  locations  into  the 
beamformer  process  is  a  progressive  degradation  in  the  direc¬ 
tivity  or  resolution  capability  of  the  beamformer  while  moving 
from  case  2  through  case  6. 

Fig.  6  investigates  a  situation  involving  two  sources  spaced 
45°  apart  in  spatial  angle.  Comments  regarding  the  two  bearing 
response  plots  for  the  six  cases  essentially  are  the  same  as 
noted  for  the  single  source  situations  in  Figs.  3-5. 
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III.  PASSIVE  RANGING 

Triangulation  and  wavefront  curvature  are  two  popular 
approaches  to  the  passive  ranging  problem  [9],  [11],  [13]. 
Although  differing  from  a  systems  engineering  viewpoint, 
they  are  quite  similar  mathematically.  Essentially,  triangula¬ 
tion  considers  two  arrays  where  the  intersection  of  their 
respective  bearing  lines  provides  source  location  information. 
Three  arrays  are  employed  in  wavefront  curvature  processing 
to  estimate  the  two  delays  from  the  forward  array  to  the 
midarray  and  from  the  midarray  to  the  after  array.  Wave- 
front  curvature  can  be  viewed  as  triangulation  by  envision¬ 
ing  the  triple  aperture  and  one  system  as  being  two  arrays, 
one  composed  of  the  forward  and  midarrays  and  one  com¬ 
posed  of  the  mid  and  after  arrays.  Both  triangulation  and 
wavefront  curvature  are  extremely  sensitive  to  locaujn  un¬ 
certainty  of  the  constituent  subarrays.  In  particular,  bowing 
or  noncollinearity  is  more  serious  than  array  displacement 
and  array  rotation  [9] . 

From  an  estimation  theory  point  of  view,  the  optimal 
solution  to  the  passive  ranging  problem  consists  of  the  maxi¬ 
mum  likelihood  (ML)  processor.  When  the  sensor  (element) 
positions  are  known,  the  ML  estimate  of  bearing  and  range 
is  obtained  by  focusing  the  individual  element  time  delays 
at  many  hypothesized  range  and  bearing  pairs  and  by  select¬ 
ing  the  pair  which  yield  the  largest  system  time-delay-and- 
sum  output  [4]-[8],  [1 0]  — [1 2] .  ML  processing  (focused 
beamforming)  is  robust  to  element  location  uncertainty 
when  the  error  is  independent  from  element  to  element. 
The  mathematics  of  focused  beamforming  is  developed  below 
along  with  several  examples  illustrating  the  influence  of  vari¬ 
ous  geometrical  parameters  on  processor  output. 


0  200 
Om  BOW 


400 

meters 

(a) 


600  8C0 


tb) 


6  m  BOW 


meters 
(c) 

Fig.  7.  Array  geometries:  17  elements. 


ately  delaying  the  individual  element  signals  to  focus  a  beam 
in  the  direction  Em  and  range  rm  yields  the  focused  beam- 
former  output 


N  /  r 

£m(0=  £  wnen(f--- 
n=l  \  . 


A.  Curved  Wavefront  Beamforming 

Consider  again  the  coordinate  system  defined  in  Fig.  1.  The 
curved  wavefront  or  focused  beamforming  task  consists  of 
generating  the  waveform  bm{t)  for  each  desired  steered  beam 
direction  Bm  and  range  of  focus  rm.  Each  bm{t)  consists 
of  the  sum  of  suitably  delayed  replicas  of  the  individual 
element  signals  en{t).  The  time  delays  compensate  for  the 
assumed  differential  travel  time  differences  between  sensors 
for  a  signal  from  the  desired  beam  direction  and  range. 

Let  the  output  of  an  element  located  at  the  origin  of 
coordinates  due  to  the  /th  source  by  sft).  A  source  from 
direction  St  and  range  r,  produces  the  following  sensor  out¬ 
puts: 


=  2 

n=  1 

(  rm-(rm2  +  \En\2-7rmBm  *£„)I/2\ 

•'■V - 7 - ) 

(5) 


-2 


HV/O  +  r,~rm  +  (r,„2  +  | E2  I 


-2 rmBm  •  En)ll2  -  (r2  +  \En'\2 

-2r,S,  *  En')ll2/r)  (6) 


en(0  =  S/  \  t  + 


(t  + 


rJ  ~  lr,S,  ~E„'  I 


r,  ~(r,2  +  |E„' I2  — 


2  rfi,  -£nV/2j 


(4) 


where  c  is  the  (constant)  speed  of  propagation  (independent 
of  location)  and  En'  indicates  the  true  nth  element  position 
which  may  differ  from  the  assumed  location  E„.  Appropri- 


where  the  vv„  are  weights  which  have  been  applied  to  each 
element  signal. 

B  Examples 

Curved  wavefront  or  focused  beamforming  for  the  six 
array  geometries  depicted  in  Fig.  7  will  be  used  as  examples 
to  supplement  the  foregoing  mathematical  derivations.  The 
three  arrays  in  Fig.  7  each  have  17  elements  with  a  constant 
interelement  spacing  of  50  m  thus  yielding  a  total  array- 
length  of  800  m.  The  second  and  third  arrays  differ  from  the 
first  in  that  the  elements  take  on  a  gradual  bow  reaching  a 
peak  offset  for  the  center  array  element  of  3  m  (Fig.  7(b)) 
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Fig.  8.  17-element  anay:  0-m 

or  6  m  (Fig.  7(c))  from  a  straight  line  connecting  the  first 
and  last  elements.  As  with  the  examples  in  Section  II,  these 
shapes  have  been  chosen  for  the  sake  of  simplicity  and  are  not 
intended  to  represent  any  particular  scenario. 

For  each  array  geometry,  four  ranges  to  a  broads5  '*  source 
(®i  =  90°,  0  =  0°)  are  investigated:  1)  rt  =5  km-  .)  r,  -  10 
km;  3)  rt  —  20  km;  and  4)  rt  =  65  km.  A  wide-band  source 
signal  was  simulated  by  the  incoherent  (power)  addition  of 
the  response  of  the  focused  beamformer  to  five  tonals  uni¬ 
formly  spread  across  the  band  of  100-500  Hz.  Three  types 
of  plots  are  presented  illustrating  the  response  of  the  focused 
beamformer  range  with  6m  =  0 ,,  azimuth  with  rm  =  rh  and 
the  three-dimensional  range-azimuth  surface. 

The  range  and  azimuth  plots  represent  orthogonal  slices 
through  the  range-azimuth  surface  which  intersect  at  the 
source  location  rt  and  For  each  source  range  rh  three 
range  and  azimuth  plots  are  provided:  1)  focused  beamformer 
response  for  the  straight  array  (£„  =  £„');  2)  focused  beam- 
former  response  for  the  bowed  array  when  the  elements  in¬ 


«  10  km 


(b) 


■  86  km 

(d) 

and  3-m  bow  (6m  •  g,  »  90°). 

correctly  are  assumed  straight  (£„  ¥=  £„');  and  3)  focused 
beamformer  response  for  a  bowed  array  whose  element 
positions  are  correctly  incorporated  into  the  focused  beam¬ 
forming  process  (£„  =  £„').  Also  provided  for  each  source 
range  rt  are  three-dimensional  range-azimuth  surfaces  for  the 
first  case  noted  above  (straight  array  with  £„  =  £„'). 

Figs.  8-12  illustrate  the  focused  beamformer  response 
for  the  1 7-element  arrays.  The  range  plots  indicate  a  distinct 
response  peak  when  the  array  element  positions  are  conectly 
incorporated  into  the  focused  beamforming  process.  When  the 
array  actually  is  bowed  but  is  assumed  straight  for  the  focused 
beamforming  operation,  a  significant  distortion  of  the  range 
response  occurs.  The  distortion  consists  of  both  a  shift  in  the 
peak  response  such  that  the  range  is  underestimated  as  well 
as  a  slow  decrease  in  response  for  ranges  beyond  the  peak 
rendering  a  determination  of  the  exact  range  where  the  peak 
occurs  more  difficult.  The  surface  of  the  focused  beamformer 
reponse  as  a  function  of  both  range  rm  and  azimuth  6m 
is  illustrated  nicely  in  Fig.  12. 
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r  ^  =  65  km 
(6) 

Fig.  12.  17-element  array:  0-m  bow  (9|  =*  90°). 


IV.  SUMMARY 

The  subject  of  this  paper  has  been  the  problem  of  beam 
formation  from  a  towed  line  array  whose  shape  has  been  dis¬ 
torted.  Essentially,  the  following  has  been  demonstrated  via 
these  examples. 

1)  Perturbations  from  an  assumed  perfectly  straight  shape 
have  serious  consequences  with  regards  to  bearing  and  range 
estimation. 

2)  When  the  actual  perturbed  array  shape  is  incorporated 
into  the  signal  processing,  desired  system  performance  is  re¬ 
gained. 

Conventional  array  processing  systems  do  not  have  the 
capability  to  account  for  the  noncollinearity  of  the  array 
elements.  Ar  approach  to  eliminate  the  undesired  loss  of 
system  performance  would  be  to  carry  out  the  beamforming 
dynamically.  As  an  example,  the  MPL  Dynamic  Beamformer 
permits  the  incorporation  of  slow  changes  in  element  positions 
and  beam  steering  directions  while  the  beamformer  carries 


out  the  real-time  formation  of  1300  beams  from  32  input 
sensors  [14],  The  element  coordinates  are  updated  periodi¬ 
cally  so  that  the  required  time  delays  for  each  element  signal 
can  be  computed  for  pure  time-delay-and-sum  beam  formation 
in  any  specified  direction.  The  original  motivation  behind 
the  fabrication  of  a  dynamically  programmable  beamformer 
by  MPL  was  the  coherent  processing  of  data  from  a  drifting 
sonobuoy  array.  In  that  program,  the  spatial  distribution  of 
the  sonobuoy  field  was  determined  by  an  active  array  element 
location  system.  For  towed  array  work,  an  alternative  to  an 
active  element  location  system  would  be  utilization  of  an 
algorithm  which  would  predict  the  element  positions  based 
on  parameter  inputs  (e.g.,  tow  point  trajectory)  which  char¬ 
acterize  the  time-evolving  shape  of  the  array. 
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ABSTRACT 

The  application  of  adaptive  least-squares  lattice  struc¬ 
tures  to  problems  in  underwater  acoustics  is  explored.  The 
focus  of  attention  is  on  the  active  sonar  problem  where  a 
strong  component  of  surface  reverberation  exists.  Both  single 
and  dual  channel  problems  are  considered.  The  prewhitening  of 
acoustic  reverberation  data  is  used  as  an  illustration  of  the 
former  while  sea  surface  reverberation  rejection  is  used  as  an 
illustration  of  the  latter. 


1.0  INTRODUCTION 

Numerous  applications  exist  which  require  a  linear  filter¬ 
ing  operation.  Often,  the  nature  of  that  filtering  task  is 
time  varying  in  some  nondeterministic  fashion  due  to  nonsta- 
tionarity  of  the  underlying  time  series.  In  such  situations,  a 
filter  which  can  adapt  to  a  changing  environment  is  needed. 

Here,  the  application  of  adaptive  least-sqnares  lattice 
structures  to  problems  in  underwater  acoustics  is  explored. 
Both  single  and  dual  channel  problems  are  considered.  The 
prewhitening  of  acoustic  reverberation  data  is  used  as  an 
illustration  of  the  former  while  sea  surface  reverberation 
rejection  is  used  as  an  application  of  the  latter. 
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Acoustic  reverberation  data  has  a  complex  and  highly  vari¬ 
able  power  spectrum.  The  spatial  transfer  function  charac¬ 
teristics  of  both  the  transmitting  and  receiving  transducers 
have  a  significant  influence  on  the  time-evolving  shape  of  the 
corresponding  range-Doppler  map.  Furthermore,  depending  on  the 
locatior  of  the  transducers  in  the  water  column,  the  ocean's 
boundaries  (both  surface  and  bottom)  can  add  major  contribu¬ 
tions  which  characteristically  have  a  sudden  onset.  Figure  1 
illustrates  graphically  the  origin  of  the  three  major  contribu¬ 
tors  (surface,  volume,  and  bottom)  to  the  reverberation  time 
series. 


As  an  illustration  of  an  interesting  application  area 
where  reverberation  is  both  the  signal  of  interest  as  well  as  a 
potential  contaminant,  consider  the  Doppler  sonar  problem. 
Doppler  sonars  use  the  perceived  shift  in  carrier  frequency 
between  the  outgoing  pulse  and  the  returning  echo  to  make  an 
estimate  of  scatterer  radial  velocity  at  several  ranges  of 
interest.  Scatterers  suspended  neutrally  buoyant  in  the  water 
column  are  used  as  tracers  for  the  remote  sensing  of  water  mass 
motion.  Essentially,  the  problem  is  one  of  high  resolution 
(distinguishing  small  differences  in  velocity)  spectral  estima¬ 
tion  across  short  data  segments  (range  bins).  Adaptive  spec¬ 
tral  estimation  techniques  can  be  used  to  follow  the  time- 
evolving  (corresponding  to  range)  spectral  characteristics  of  a 
returning  echo.  A  simplistic  model  of  the  corresponding 
range-Doppler  map  (a  three-dimensional  surface)  consists  of  a 
single  symmetrical  hump  whose  track  as  a  function  of  time  is 
indicative  of  the  radial  velocity  of  scatters  in  successive 
range  cells.  However,  a  more  careful  consideration  of  the 
sonar  transducer's  spatial  response  characteristics  suggests 
that  surface  reverberation  can  be  a  significant  contaminant  in 
the  range-Doppler  map. 
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2.0  HIE  COMPLEX  ADAPTr.E  JOINT  PROCESS  ALGORITHM 

The  general  problem  of  interest  is  dual  channel.  As  will 
be  seen,  the  single  channel  algorithm  is  embedded  in  the  dual 
channel  structure.  The  direct  form  implementation  of  the  joint 
process  structure  is  illustrated  in  Figure  2.  The  current  and 


Figure  2.  Direct  form  realization  of  the  joint  process  structure. 

most  recent  p  samples  of  the  reference  channel  process 
{x(n) ,x(n-l) , . . . ,x(n-p)}  are  linearly  combined  to  form  an  esti¬ 
mate  of  -d(n).  Of  particular  interest  in  noise  cancelling 
applications  is  the  residual  ed  (n)  obtained  by  subtracting  the 
filtered  reference  channel  fronrthe  primary  channel  16].  The 
corresponding  transfer  function  of  the  reference  channel  filter 
is  given  by 

HU)  =  2  h  z~k  (1) 

k=0 


Note  that  if  h_=C  and  d(n)-x(n),  the  filter  in  Figure  2 
becomes  the  one-step  forward  linear  predictor  illustrated  in 
Figure  3(a).  Removing  the  predictable  components  from  x(n) 
yields  the  forward  prediction  error  sequence  ed  (n).  The 
corresponding  transfer  function  of  the  forward  prediction  error 
filter  is  given  by 

A(z)  =  1  a.z_k  ,  a  =  1  .  (2) 
k=0  *  0 
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Figure  3.  (a)  One-step  forward  prediction  error  filter. 

(b)  One-step  backward  prediction  error  filter. 

A  companion  to  Figure  3(a)  is  the  one-step  backward  linear 
predictor  (coefficients  b  ,b.,...,b  j)  shown  in  Figure  3(b) 

along  with  the  backward  prediction  erfor  sequence  r  (n). 

P 

The  forward  and  backward  linear  predictors  can  be  realized 
equivalently  in  the  form  of  a  lattice  structure.  Shown  in  Fig¬ 
ure  4,  the  lattice  parameters  K*\  and  are  known  either  as 
reflection  coefficients  or  partial  correlation  coefficients. 
The  e.(n)  and  ri(n)  are  the  ith  order  forward  and  backward 
prediction  error  sequences,  respectively.  Thus,  the  pth  order 
linear  predictor  is  created  on  a  stage-by-stage  basis  a  single 
order  at  a  time. 

As  shown  in  Figure  S,  the  lattice  realization  of  the  for¬ 
ward  and  backward  linear  predictors  can  be  embedded  in  the 
joint  process  structure.  Additional  parameters  K  .  weight  the 
backward  prediction  error  residuals  r.(n).  As  with  the  lattice 
itself,  the  pth  order  joint  process  structure  is  created  on  a 
stage-by-stage  basis  a  single  order  at  a  time. 


A  continuously  adapting  approach  to  the  realization  of  the 
time  varying,  complex  joint  process  structure  is  summarized 
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«  ,  In) 

Figure  4.  (a)  Forward  and  backward  prediction  error  filters, 

(b)  The  ith  stage  of  the  lattice. 


Figure  S.  Lattice  realization  of  the  joint  process  structure. 

below.  The  joint  complex  least-squares  lattice  (JCLSL)  is 
based  on  a  recursive  in  time  solution  for  the  filter  whose 
transfer  function  is  ll(z).  Bounded  by  [0,1],  a  parameter  a  is 
found  in  the  algorithm  which  adjusts  the  time  constant  of  its 
exponentially  decaying  memory  of  past  data.  Large  a's  imply 
short  time  constents.  Additional  material  related  to  this 
algorithm  can  be  found  in  [1-5]  and  the  references  cited 
therein. 

The  least  squares  lattice  algorithm  is  summarized  as  fol¬ 
lows  (see  Figure  5): 
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Initialization  (i  =  0,  !••••  >  P) 


r^-1)  =  0  ,  i*P 

E'(-i)  =  *JCLSL  »  eJCLSL 

A.(-l)  =  0  ,  i*0 


0.001  and  i*p 


b(i)(_D  =  o  ,  0<k£i-l  .  i^O.  and  i*P 

=0  , 

K*f  ( i-1 )  =  0 
1 

Ti-ic  ..date  (n  2  0) 


(3a) 

(3b) 

(3c) 

(3d) 

(3e) 

(3f ) 


c0(n)  =  r0(n> 

Ee (n)  =  E*(n> 

o  ° 

r_j(n)  =  o 
e“1(n)  =  d(n) 


x(n) 

(1 '  *jclsi.>  Eo(n_1)  +  ix(n,i2 


( 3  g ) 


<3h) 


(3i) 


( 3  j ) 


Order  update  (i  v  0,1 


P> 
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Lattice 

A^n) 

ei-l 
1  - 

K*(n) 

Kf(a) 

ei(n) 

^(n) 

E*(n) 

E*(n) 

Vl(n) 

A?(n) 

Kj(a) 


_  (1  aJCLSL)  Ai(n-1) 


(3k) 


(nJrj.jCn-l) 


i#0 


=  A*(n)  /  ,  i^O 

=  Ai(n)  /  Ef.jtn-l)  ,  i#0 

=  e^tn)  +  KjlUJr^U-l)  ,  i#0 

=  r^U-l)  +  KjCnJe^U)  ,  i#0 

=  E^U)  -  lA^n)]2/  E^n-l)  ,  i#0 

=  ET_1(n-l)  -  jAitn)  j2  /  E®__1(n)  ,  i#0 

-  Yi_2<B>  +  |ri-l(n>|2  '  EI-i(n)  •  l*° 

e?  , (n)r?(n) 

=  (1  -  °jaSL>Ai(n_1)  -  1  1  - 


(31) 

(3m) 

(3n) 

(3o) 

(3p) 

(3q) 

(3r) 

(3s) 


A?(n> 

e5(») 


(3t) 
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ej(n)  =  eJ_j(n)  +  Kj(n)r.(n)  (3u) 

Predictors 

aj^ta)  =  Ka(n)  (3v) 

b^U)  =  Z°(n)  (3w) 

o  1 

‘ki)<n)  =  aiil>(a)  +  Ki(a)  +  ik-i1>(a_1)  (3x) 

iikii-i 

b<i}(n)  =  b^'^U-l)  +  Kf(n)  a* i_1) <n>  (3y) 


3.0  APPLICATION  EXAMPLES 


3.1  Prewhitening 

As  an  example  of  a  single  channel  problem,  the  prewhiten¬ 
ing  of  acoustic  reverberation  will  be  considered.  The  results 
of  processing  a  ping  selected  from  an  experimental  ocean  rever¬ 
beration  data  set  will  be  used  for  illustration.  Of  interest 
are:  (1)  the  inputs  to  the  prewhitening  filters,  (2)  the 
prewhiteners  themselves,  and  (3)  the  outputs  of  the  prewhiten¬ 
ing  filters.  The  first  0.3  s  of  the  complex  basebanded  ping 
under  study  has  been  set  to  zero  so  that  transmit  pulse  leakage 
into  the  receiver  at  the  time  of  data  collection  would  no.'* 
influence  the  processing.  Since  the  goal  of  a  prewhitening 
filter  is  to  produce  an  output  with  a  'white*  spectrum,  it  fol¬ 
lows  that  the  inverse  of  the  spectrum  of  the  prewhitening 
filter  is  an  estimate  of  the  spectrum  of  the  input  time  series. 
Such  time-evolving  spectral  plots  will  be  compared  with  conven¬ 
tional  range-Doppler  maps  of  the  input  time  series. 
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Figure  6. 


Plots  characterizing  the  performance  of  the  complex 
least-squares  lattice  prewhitening  structure  are  provided  in 
the  following  order.  First,  Figure  6  is  a  conventional  range- 
loppler  map  of  the  ping  under  study.  Then,  a  set  of  plots 
explores  the  characteristics  of  the  prewhitening  structure 
(Figures  7-12).  The  plots  are  presented  in  horizontal  pairs  - 
the  plot  on  the  left  being  the  inverse  of  the  prewhitening 
filter  (an  adaptive  spectral  estimate)  and  the  plot  on  the 
right  being  the  conventional  range-Doppler  map  of  the  output  of 
the  prewhitening  filter. 

Trade-offs  are  illustrated  as  a  function  of  the  adaptation 
rate  parameter  (a).  The  values  of  a  were  chosen  as  follows. 
First,  a  mid-range  value  of  a  was  determined.  Then,  values  an 
order  of  magnitude  larger  and  an  order  of  magnitude  smaller 
were  added  to  the  list  of  cases  to  process. 

3.2  Surface  Reverberation  Rejection 

As  an  example  of  a  dual  channel  problem,  sea  surface 
reverberation  rejection  will  be  considered.  In  an  active  sonar 
system,  transmitted  energy  reflected  off  the  sea  surface  typi¬ 
cally  is  discriminated  against  by  the  sidelobe  characteristic 
of  the  receiving  transducer.  However,  even  though  attenuated, 
this  energy  still  can  be  a  significant  contributor  to  the  noise 
background  level.  Here,  an  adaptive  joint  process  structure  is 
used  which  rejects  boundary  reverberation  by  taking  advantage 
of  the  spatial  separation  between  returning  echoes  of  interest 
and  transmitted  energy  reflected  ofx  the  sea  surface  [6-7]. 

The  results  of  processing  another  ping  selected  from  the 
the  experimental  ocean  reverberation  data  set  mentioned  above 
now  will  be  presented.  The  sonar  transducer  consisted  of  a 
number  of  individual  elements.  For  this  ping,  row  sums  of  the 
elements  were  available.  The  prim  ry  channel  signal  (d(n)  in 
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Figure*  2  and  5)  wa*  obtained  as  a  aiaple  ana  of  these  row  sun* 
(i  e  a  forward  looking  beam) .  The  reference  ehannel 
(x*(n)  in  Figure*.  %  and  5)  wa*  obtained  by  differenoing  two 
adjacent  row.V.e.  a  be«  with  a  null  in  the  forw.rd  look  n, 
direction  and  large  lobe*  towards  the  surface  and  the  bottoa). 


*doc’>»*  •prc''ot  *»» i-o'*  1  U  iNo-mH.I 

cut  rt  •cUi-t  Mf? 
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Figure  7- 


Figure  8. 
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Figure  9. 
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Figure  10. 


APPLICATIONS  01  ADAPTIVI.  ARRAY  PROCI.SSING 


457 


The  corresponding  input  and  output  range-Doppler  naps  for 
the  joint  process  structure  are  provided  in  Figures  13-15. 
Both  the  primary  and  reference  channel's  range-Doppler  maps  are 
heavily  contaminated  by  surface  reverberation.  In  contrast, 
the  output  range-Do; pier  map  shows  a  substantial  reduction  in 
this  surface  reverberation  contamination. 


4.0  SUMMARY 

The  application  of  adaptive  least-squares  lattice  '  struc¬ 
tures  to  problems  in  underwater  acoustics  has  been  discussed. 
The  prewhitening  of  acoustic  reverberation  data  was  used  as  an 
illustration  of  a  single  channel  problem.  Then,  sea  surface 
reverberation  rejection  was  used  as  an  illustration  of  a  dual 
channel  problem. 
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DISCUSSION* 

Comment  J  W.  R.  Griffiths 

Would  you  comment  on  the  computational  complexity  ot  adaptive  lattice  struc¬ 
tures7 

Reply  W  S  Hodgkiss 

The  lattice  structure  naturally  falls  out  of  the  solution  to  the  least-squares  problem 
At  the  heart  of  things  is  Levinson's  algorithm  which  has  computational  complexity 
on  the  order  of  pJ  for  a  p-stage  filter. 

Comment  C  van  Schooneveld 

What  is  the  loss  of  preformance  due  to  pursuing  this  scheme  of  adaptive  beam¬ 
forming  as  apposed  to  an  approach  which  is  "optimal"  from  a  decision-theoretic 
standpoint7 

Reply  W  S  Hodgkiss 

Your  question  points  to  an  area  deserving  of  additional  work.  In  general,, compari¬ 
sons  have  not  beer-  -J e  of  performance  between  an  adaptive  structure  (such  as 
reported  here)  wi,  nears  to  be  doing  a  good  job  "locally"  and  a  globally  opti¬ 
mal  processor  woo  the  same  problem. 

Comment  G.C.  Carter 

It  is  apparent  that  the  overall  cancellation  ratio  (conventional  beam  output  power 
divided  by  adaptive  beam  output  power)  is  rarely  more  than  a  few  dB.  Is  the 
gain  worth  the  effort7 

Reply  W  S  Hodgkiss 

The  benefits  of  adaptive  beamforming  are  selective.  In  the  results  presented,  not 
all  regions  of  the.  range-Doppler  map  were  affected  equally  For  example,  a  large 
spike  of  reverbcation  coming  in  a  sidelobe  of  the  transducer  was  removed  as  well 
as  a  significant  amount  of  the  boundary  reverberation  beyond  1.2  s  in  range 
Whether  the  gam  is  worth  the  effort  or  not  is  critically  dependent  on  the  actual 
problem  being  worked  on 
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Boundary  reverberation  rejection  via  constrained  adaptive 
beamforming 
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The  problem  of  selective  reverberation  cancellation,  whereby  both  “signal”  and  “noise”  are 
constituent  components  of  the  received  reverberation  process,  is  the  focus  of  this  article.  The 
proposed  solution  involves  the  application  of  a  constrained  adaptive  beamforming  technique. 

The  “prewindowed”  deterministic  least-squares  lattice  filter  is  used  as  the  central  adaptive 
element.  Constraints  are  in  the  form  of  simple  spatial  filtering  prior  to  adaptation.  The  spatial 
correlation  characteristics  of  volume  and  boundary  reverberation  are  found  to  be  directly 
applicable  in  a  reverberation  cancellation  context.  Experimental  verification  is  offered  by 
processing  reverberation  data  from  a  shallow-water  deployment  of  an  active  sonar  system.  It  is 
shown  that  the  boundary  reverberation  components  can  be  effectively  suppressed  while 
preserving  the  volume  return.  Computer  simulations  of  the  experiment  offer  additional  insight 
into  the  adaptation  process. 

PACS  numbers:  43.60.Gk,  43.30.Gv,  43.30. Vh 


INTRODUCTION 

The  ocean  abounds  with  objects  that  can  intercept  and 
reradiate  acoustic  energy.  Suspended  sediment,  organic  de¬ 
tritus,  air  bubbles,  plankton,  fish,  and  minute  discontinuities 
in  the  thermal  structure  are  all  capable  of  redirecting  sound. 
Irregularities  of  the  sea  surface  and  the  ocean  floor  are  also 
significant  contributors  to  this  reradiation  of  sound  known 
as  scattering.  The  composite  echo  from  all  scatterers  is 
known  as  reverberation.  When  the  objective  is  signal  detec¬ 
tion,  reverberation  is  clearly  a  form  of  noise.  However,  rever¬ 
beration  can  also  be  a  valuable  information-bearing  signal;  it 
offers  information  about  the  nature  and  distribution  of  the 
scatterers  and  to  the  extent  that  the  scatterers  are  influenced 
by  a  fluid  process,  about  the  process  itself.  The  information, 
if  properly  extracted,  can  be  used  to  quantify  fishery  species 
of  commercial  interest  (Holliday,  1974),  to  examine  plank¬ 
tonic  communities  (Greenblatt,  1980),  to  monitor  pollution 
in  industrial  dumping  sites  ( Orr  and  Hess,  1 978 ),  to  identify 
ocean  bottom  types  (deMoustier,  1985),  or  for  remote  sens¬ 
ing  of  oceanic  fluid  processes  (Pinkel,  1981 ).  As  reverbera¬ 
tion  is  being  reevaluated  in  terms  of  its  information  content, 
so  must  the  signal  processing  techniques  used  vis-a-vis  rever¬ 
beration  be  reconsidered.  For  instance,  rather  than  indis¬ 
criminate  reverberation  suppression,  the  situation  may  call 
for  the  extraction,  from  a  composite  reverberation  return,  of 
the  component  created  by  the  class  of  scatterers  associated 
with  the  process  of  interest.  In  this  scenario,  “signal”  and 
“noise”  are  both  constituent  components  of  the  received  re¬ 
verberation  process.  Here,  we  offer  a  potential  solution  to 
this  signal  processing  problem  based  on  an  adaptive  beam¬ 
forming  technique. 

As  their  name  implies,  adaptive  filters  are  capable  of 
responding  to  changing  conditions  through  a  iudimentary 
“learning”  process.  The  filter  parameters  adjust  themselves 


under  the  guidance  of  an  appropriate  cost  function  and  track 
the  evolving  characteristics  of  the  input  signal(s).  These 
structures  grew  out  of  the  demand  for  systems  capable  of 
operating  in  uncertain,  time- varying  environments  and  we'  \ 
made  possible  by  the  increasing  availability  of  computa¬ 
tional  power..  They  can  be  viewed  as  realizable  approxima¬ 
tions  to  optimal  Wiener  filters  when  only  a  single  sample 
function  of  the  “signal”  and  “noise”  processes  is  available. 
In  this  case,  the  optimum  filtering  problem  may  be  ap¬ 
proached  from  a  statistical  viewpoint  by  invoking  stationar- 
ity  and  ergodicity  and  performing  time  averaging,  or 
through  a  deterministic  least-squares  formulation.  Both  so¬ 
lutions  can  be  implemented  either  in  a  block-processing  or  a 
time  recursive  mode  and  both  give  rise  to  efficient  lattice 
structures.  Nonstationary  processes  are  treated  as  locally 
stationary  by  restricting  the  time  interval  over  which  opti¬ 
mization  is  performed.  An  adaptive  filter  solves  the  Wiener 
problem,  or  its  deterministic  counterpart,  over  the  restricted 
optimization  interval  controlled  by  an  adaptation  coeffi¬ 
cient,  which  provides  for  a  “fading”  of  past  information  in 
favor  of  recent  values.  Adaptive  algorithms  can  be  used  to 
control  the  element  weights  of  an  array,  changing  its  beam 
pattern  in  an  optimum  manner  to  reject  interference.  This 
class  of  adaptive  filtering  is  known  as  adaptive  beamforming. 

I.  CONSTRAINED  ADAPTIVE  BEAMFORMING 

Conventional  arrays  £i  -  usually  designed  under  the  as¬ 
sumption  that  noise  is  spatially  disorganized.  These  fixed- 
weight  arrays  suffer  a  degradation  of  performance  in  the 
presence  of  dir  tional  interference  possessing  a  degree  of 
spatial  coherence.  Such  interference  may  be  due  to  natural 
sources  or  reverberation  returns.  Other  factors,  such  as  ar¬ 
ray  motion,  mul'ipaths,  and  constantly  changing  interfer¬ 
ence  characteristics  contribute  to  further  deterioration. 
Ad;  p  ..ve  arrays,  on  the  other  hand,  are  capable  of  reducing 
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or  eliminating  directional  noise  components  and  of  respond* 
ing  to  changing  conditions  by  adjusting  their  pattern  re¬ 
sponse  according  to  an  appropriately  chosen  error  criterion. 
In  general,  adaptive  arrays  are  in  the  form  of  a  space-time 
filter.  The  exact  manner  in  which  control  of  the  available 
spatial  and  temporal  degrees  of  freedom  is  relegated  to  the 
adaptive  processor  and  the  choice  of  a  performance  measure 
cepends  largely  on  the  objectives  and  priorities  of  the  specif¬ 
ic  application.  A  good  tutorial  introduction  to  the  adaptive 
array  problem  is  given  by  Gabriel  ( 1976 ) .  Our  objective  here 
is  to  utilize  the  general  adaptive  algorithms  derived  under 
the  minimum  mean-square  error  ( MMSE)  criterion,  in  con- 
juction  with  the  adaptive  noise  canceling  ( ANC)  principle, 
in  order  to  eliminate  reverberation  interference  entering 
through  the  sidelobes  and  the  mainlobe  edges  of  the  receiv¬ 
ing  array. 

In  an  adaptive  array  application,  it  is  important  to  uti¬ 
lize  all  available  a  priori  information  about  signal  and  inter¬ 
ference  in  order  to  ensure  signal  preservation  at  the  output. 
Constraints  may  be  incorporated  into  the  adaptive  algo¬ 
rithm  to  maintain  a  chosen  frequency  response  in  a  desired 
direction  (Frost,  1972).  “Pilot’-  signals,  simulating  actual 
signals  of  interest,  have  also  been  used  to  describe  a  desired 
“main”  look  direction  through  a  two-mode  adaptation  pro¬ 
cess  (Widrow  et  ai,  1967).  An  alternative  is  to  impose  pre¬ 
adaptation  constraints  in  the  time,  space,  and  frequency  do¬ 
mains  ( Applebaum  and  Chapman,  1976).  Restrictions  may 
be  applied  in  the  time  domain  by  allowing  the  array  to  adapt 
only  when  no  signal  is  present;  in  the  frequency  domain,  by 
allowing  adaptation  only  to  energy  received  outside  the  sig¬ 
nal  band;  or  preadaptation  spatial  filtering  may  be  used  to 
remove  the  signal  from  the  reference  channel,  thus  protect¬ 
ing  it  from  cancellation.  Such  prefiltering  may  range  from 
complete  conventional  beamforming  (Fig.  1(a)]  to  simple 
element-to-element  subtraction  [Fig.  1(b)].  The  former  can 
be  applied  only  when  the  spatial  characteristics  of  the  inter¬ 
ference,  as  well  as  the  signal,  are  known  and  requires  an 
external  steering  mechanism.  The  latter  is  more  appropriate 
when  the  directional  characteristics  of  the  interference  are 
unknown,  or  variable,  and  results  in  ( constrained )  spatially 
adaptive  cancellation  beams. 

Such  spatial  prefiltering  operations,  by  eliminating  the 
signal  of  interest  from  the  reference  channel (s),  in  effect 
place  constrained  adaptive  beamforming  in  the  context  of 
ANC.  This  approach  affords  us  the  flexibility  of  using  any 
efficient  solution  to  the  general  filtering  problem  as  the  adap¬ 
tive  processor  controlling  the  array  element  weights. 

A.  The  adaptive  noise  canceling  (ANC)  concept 

The  noise  canceling  solution  stems  from  the  following 
dual-channel  filtering  problem;  Given  two  discrete  time  sto¬ 
chastic  processes  x(rt),  d(n),  estimate  the  second  process 
( primary  signal )  by  operating  on  the  first  process  ( reference 
signal).  The  MMSE  solution  is  given  by 

R^(n,fc)  =  £  F^“(n)R«(M,/c),  n  - p<.k<n,  (1) 

i-O 

where  Fj,' ’  is  a  linear,  causal  filter  of  order  p  and  RXI,Rdx  are 
the  autocorrelation  function  of  the  reference  signal  and  the 


cross-correlation  function  between  the  primary  and  refer¬ 
ence  signals,  respectively. 

The  set  of  circumstances  giving  ris  t  to  the  noise  cancel¬ 
ing  concept  is  illustrated  in  Fig.  2(a).  Ine  primary  channel 
consists  of  the  signal  s  corrupted  by  a  form  of  additive  noise 
n0  and  the  reference  channel  consists  of  a  process  n ,  related 
in  some  unknown  way  to  the  primary  noise.  The  key  require¬ 
ment  is  that  the  signal  be  uncorrelated  with  both  the  primary 
noise  and  the  reference  process: 

0,  R„t  =0. 

We  then  have 

^■dx  ~  R-n^ 

which  leads  to  the  noise  canceling  solution  as  a  special  case 
of  the  general  filtering  problem  solution; 

R^,  («,*)  =  £  («)/!,,,,  («,*),  n-p<k<,n.  (2) 

i-O 

Solving  the  filtering  problem  in  this  setting  is  equivalent  to 
producing  the  best  MMSE  estimate  of  the  primary  noise  pro¬ 
cess.  The  ANC  output,  obtained  by  subtracting  this  estimate 
from  the  primary  input,  will  consist  of  the  signal  component 
s  plus  a  residual  error  n0  —  n0 . 

A  more  detailed  model  for  the  noise  canceling  structure 
is  shown  in  Fig.  2(b) .  The  “mismatch”  between  the  primary 
and  reference  inputs  is  represented  by  a  linear  transfer  func¬ 
tion  H(z).  Additional  uncorrelated  noise  components  m0 
and  m,  are  included  in  the  two  channels.  This  particular 
model,  which  is  representative  of  many  situations  of  practi¬ 
cal  interest,  has  been  studied  in  detail  by  Widrow  et  al. 
( 1975),  who  determined  that  the  uncorrelated  noise  compo¬ 
nents  have  a  deleterious  effect  on  cancellation.  Other  factors 
limiting  ANC  performance  are  the  presence  of  signal  com¬ 
ponents  in  the  reference  channel  and,  somewhat  surprising¬ 
ly,  a  high  signal-to-noise  ratio  in  the  primary  channel. 


B.  Adaptive  algorithms 

Rather  than  evolving  from  a  general  mathematical  pre¬ 
mise,  early  adaptive  filters  were  designed  to  solve  specific 
real-time  engineering  problems  (e.g.,  Glaser,  1961;  Davis¬ 
son,  1966;  Anderson,  1969).  The  work  of  Widrow  and  Hoff 
( 1960)  produced  the  first  generalized  adaptive  structure, 
based  on  the  method  of  steepest  descent,  known  as  the  least 
mean-square  (LMS)  algorithm.  The  LMS  has  been  used  as 
the  centrai  adaptive  processor  in  applications  such  as  adap¬ 
tive  array  processing  (Widrow  et  al,  1967;  Griffiths,  1969; 
Frost,  1972),  adaptive  modeling  (Schade,  1971),  adaptive 
channel  equalization  (Lucky,  1965),  and  adaptive  noise 
canceling  (Widrow  etal.,  1975).  Adaptive  lattice  structures 
were  shown  to  outperform  the  LMS  by  Satorius  and  Alex¬ 
ander  ( 1979)  and  Hodgkiss  and  Presley  ( 1981 )  in  channel 
equalization  and  frequency  tracking  applications,  respec¬ 
tively.  Lattice  joint-process  filters  suitable  for  noise  cancel¬ 
ing  have  been  suggested  (Griffiths,  1978),  but  few  applica¬ 
tions  in  cases  of  practical  interest  have  been  reported  (e.g., 
Alexandrou,  1985b). 
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FIG.  1.  Adaptive  beamforming.  (a) 
Fixed,  preformed  reference  beam.  An 
external  steering  mechanism  is  required 
to  track  the  changing  directional  proper¬ 
ties  of  the  interference.  (b)  Constrained 
adaptive  beamformer.  The  multichannel 
adaptive  filter  weights  the  constrained 
elements  to  form  a  spatially  adaptive  ref¬ 
erence  beam. 


II.  DETERMINISTIC  LEAST-SQUARES  ERROR  (LSE) 
SOLUTIONS 

The  development  of  the  LMS  and  its  successful  applica¬ 
tion  in  a  plethora  of  experimental  settings  have  demonstrat¬ 
ed  that  adaptive  structures  of  considerable  diversity  can  be 
realized  through  a  central  adaptive  estimator  which  is  a 
practical  implementation  of  a  Wiener  filter.  In  applications 
of  practical  interest,  where  real-time  processing  Is  of  the  es¬ 
sence,  only  a  single  sample  function  of  the  underlying  ran¬ 
dom  process  is  typically  available  and  time  averaging  must 
be  used.  When  wide  sense  stationarity  ( WSS)  holds  at  least 
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locally  and  ergodicity  can  be  safely  invoked  this  will  still  lead 
to  the  optimum  solution.  These  assumptions  form  the  basis 
for  the  development  of  stochastic  approximation  adaptive  al¬ 
gorithms  such  as  the  LMS  and  the  gradient  lattice  (GRL) 
filter  (Makhoul,  1978).  During  transition  periods,  when  the 
statistics  of  the  signal  change  drastically,  all  stochastic  ap¬ 
proximation  solutions  are  clearly  suboptimal.  In  addition, 
the  misadjustment  noise  of  the  gradient  methods  is  often  a 
limiting  factor  on  the  algorithm  convergence  speed. 

An  alternative  to  the  stochastic  approximation  methods 
is  the  deterministic  least-squares  approach.  When  WSS 
holds,  the  solutions  derived  through  this  approach  perform 
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and  with 

ep,./(n)  =  d(n)-d(n) 


we  seek  to  minimize 

HU.f)  =tT'2e/,il'/{n)eZJ/(n), 

i 

where  i,f  will  be  specified  later  (N,  </  <Nf,  i 

The  LSE  solution  to  the  filtering  problem  is  easily  ob¬ 
tained  through  the  projection  theorem  and  is  given  by 

Fj(/,/)Rv  (/,/)  =  R*,  (4  A  (3a) 

where 

R*.p(4/>  =  2  ^nn-prfnn-pt  =  Xp  (/  ,/)X"(4A 


1  REFERENCE'  i - 1  I  . 

CHANNEL! _ I 

FIG.  2.  The  adaptive  noise  ca  ;elmg  concept,  (a)  General  structure,  (b) 
Widrow  era/' s  ( 1975)  linear  model.  Here,  m0  and  m,  represent  uncorre- 
lated  additive  noise  components  in  the  two  channels  The  “mismatch" 
between  the  primary  and  reference  channels,  expressed  by  H(z),  must  be 
compensated  for  by  the  adaptive  filter  transfer  funct'on  A(z). 

identically  to  their  stochastic  (time  average)  counterparts. 
However,  they  enjoy  an  important  advantage  in  that  they 
continue  to  be  optimum  (in  a  deterministic  least-squares 
sense)  even  during  abrupt  changes  in  the  signal  statistics. 
This  contributes  to  transient  behavior  superior  to  the  gradi¬ 
ent  solutions  and  makes  the  least-squares  structures  excel¬ 
lent  candidates  for  our  choice  of  adaptive  algorithm.  For 
instance,  an  algorithm  of  this  class  was  shown  to  converge 
successfully  in  the  presence  of  amplitude  transients  that 
drove  the  LMS  and  the  GRL  to  instability  (Alexandrou, 
1985a).  This  is  relevant  here  because  the  onset  of  boundary 
reverberation  in  a  volume  reverberation  backround  can 
readily  amount  to  an  abrupt  intensity  “step.” 

The  deterministic  least-squares  problem  is  formulated 
as  follows:  Given  t\  o  discrete  time  series, 

j>(;i),d(n)},  N,<n<N/t 

of  dimension  m  and  /,  respectively,  we  form  the  estimate  of 
d(n ): 

d{n)  =  £  Fljf){n)x(n -j)  =  Fj(n)x(nn_p), 
where 

F TP  =  {F™,F"\...fP) 

I  x(n)  \ 
x(n-l)\ 


X{n-n  —  p) 


x,  (/,/)=  .  , 

\x(i—p)  ■■■  x{f-p)} 

R  dx.,,U,f)  =&lr>X?U,f)- 
The  minimum  least-squares  error  is 
Edp{i,f)  -  min{F) 

=  *„.oo (4/>  -FX UA  (3b) 

Similarly,  the  LSE  one-step  predictor  is  given  by 

A;(i,/)R  ,„(/,/>  =  [£J(4A  0 . 0],  (4) 

where 

aJ(/,/)  =  [imA  (Pp)a,/)], 
e;  (/,/) = rx_  oo  (/,  f)  -  a;q,.p  (4  A 

i 

Q*.,  =  X 

n  —  i 

Note  that  the  exact  form  of  Xp(i,/)  depends  on  the 
choice  of  i  and /.  The  following  choices  are  most  commonly 
made,  resulting  in  sample  correlatior  matrices  R  with  differ¬ 
ent  symmetry  and/or  shift  invariance  properties.  Such  prop¬ 
erties  are  crucial  to  the  development  of  algorithmic  imple¬ 
mentations  of  the  above  solutions. 

( 1 )  Prewindowed  case: 
i  =  N„  /=  Nf. 

(2)  Nonwindowed  case: 
i~N,+p,  f-  Nf. 

(3)  Prewindowed  and  postwindowed  cases: 
i  =  N„  f=N,+p. 

A.  The  “prewindowed”  least-squares  lattice  (LSL) 

In  this  case,  for  N,  =  0  and  Nf  =  N,  the  sample  correla¬ 
tion  matrix  R  of  order  p  is  of  the  form 

R p.n  =  xPi/vXpW, 
where 

/x(0)  •••  x(p)  ■■■  x(N)  \ 

°  . 

\  0  0  x(0)  •••  x(N-p)/ 

Note  that  Rp  N  is  not  Toeplitz,  but  it  consists  of  a  prod¬ 
uct  ot  two  (upper  and  lower  triangular)  Toeplitz  matrices 
and  might  be  thought  of  as  being  “near  Toeplitz.”  In  fact,  it 


1277 


J.  Acoust.  Soc.  Am..  Vol.  82,  No.  4,  October  1987 


Dimitri  Alexandrou:  Boundary  reverberation  rejection 


1277 


falls  within  the  class  of  the  so-called  a-stationary  matrices 
which  can  be  expressed  as 

R  =  T+fa,U,U»  a,=  ±  1,  (5; 

i=  i 

where  T  is  Toeplitz  and  the  {U, }  are  upper  triangular  Toe- 
plitz  matrices  (Friendlander  et  ah,  1979).  These  matrices, 
also  known  as  “low  shift  rank”  matrices,  have  certain  shift- 
invariance  properties  that  lead  to  efficient  recursive  solu¬ 
tions.  Specifically  in  the  prewindowed  case,  Rp  A  is  a -sta¬ 
tionary  with  a  =  2  and  satisfies  the  following  recursive 
identities: 


order  and  time  update 

t.oo  (^) 


(RxM(N)  Qxp  (TV)  \ 
VQ?,(AO  Rx.p_,(Ar-l)j’ 


order  update: 

R,„(A’) 


o  =  (R— 

\  V„(. 


_,(A0  V"  (A7) 


x.P^'A 

x.p  (N)  Rx.pp(N)/ 


time  update: 

Rx.pUV)  =  Rxp(N  —  1)  +  x(.v_P.Mx"WN, 
where 


(6a) 

(6b) 

(6c) 


yx.p(K)=  £  x{n  -p)x" 

R  =  0 

and  all  other  quantities  are  as  previously  defined.  These 
identities  form  the  basis  for  the  exact  order-  and  time-update 
recursions  for  the  prewindowed  prediction-error  LSL.  The 


joint-process  LSL  car.  be  obtained  through  an  embedding 
technique,  whereby  the  joint  sample  correlation  matrix 
RJp  (AO  is  defined  as 

Rz.p=Zp(A0Z»(A0, 

where 


(d{  0) 

d{p)  • 

d(N) 

„  ran  I  *<°> 

xip)  ■ 

xiN) 

Z,W=  o 

\  0 

o 

o 

■■  x(N-p), 

Then,  the  filtering  normal  equations  can  be  reexpressed  in 
terms  of  R,p  (AO  as  follows: 

FJpRz,p(A0  =  [.EP(A0,0,...,0],  (7) 

where 

F/p  (N)  —  [//,Fp  (AO  ]  • 

Here,  Rxp  has  shift  properties  milar  to  Rx  p .  The  two 
sets  of  recursions  can  be  combined  to  concurrently  solve  the 
forward  prediction,  backward  prediction,  and  filtering  prob¬ 
lems,  giving  rise  to  the  joint-process  LSL.  The  needed  alge¬ 
braic  manipulations  are  given  by  Lee  (1980).  A  schematic 
diagram  of  the  scalar  LSL  is  given  in  Fig.  3. 

1.  Exponential  weighting 

An  adaptive  version  of  the  LSL  filter  can  be  easily  de¬ 
rived  by  introducing  an  exponential  “window”  in  the  error 
i  orm.  Specifically,  one  may  choose  to  minimize 


(a)  r0(n)  r,(n)  r;(n)  v,{n) 


r„(n) 


e„(n) 


■d„(n) 


r„(n) 


e„(n) 


FIG.  3  The  scalar  LSL  filter,  (a)  Prediction 
error  (widening)  filter:  k‘  and  k’  are  the 
forward  ar  '  backward  reflection  (PARCOR) 
eoefficien  (n)  and  rjtt)  are  the  rth-nr- 
der  forward  and  backward  prediction  error 
sequences  (b)  Joint-process  filter;  the 
r,  (n  J  constitute  an  orthogonal  basis  for  the 
reference  process  x(it)  and  are  weighted  by 
the  cross-channel  coefficients  k  f  to  produce 
an  estimate  of  the  primary  process  d(n) 
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|;UV)=tr  ][(l-aUiL)'v  %■*(*)<*("), 

n  -  0 

where  is  a  constant>0,  so  that  past  errors  wili  have  a 
smaller  influence  on  the  estimate.  This  constant  “fade”  fac¬ 
tor  is  easily  incorporated  into  the  lattice  recursions  and  con¬ 
trols  the  effective  “averaging  interval,”  i.e.,  the  number  of 
samples  over  which  the  algorithm  solves  the  filtering  prob¬ 
lem  a?  any  given  instant.  This  "effective”  number  of  samples 
may  be  defined  in  terms  of  a  “time  constant”  as  follows: 

( 1  -  Olsl  )' *'■ - "  =  exp  -  [  ( N -  n )/NcS ] . 

Therefore, 

Neg  =  min{iV,  -  l/ln(  1  -  aUL )}. 

2.  Summary  of  the  LSL  algorithm 

The  complex,  multichannel,  joint-process,  exponential¬ 
ly  weighted,  prewindowed  least-squares  lattice  algorithm  is 
summarized  as  follows: 


initialization  (i  =  0, 


r,(  ■ 

-  1 )  =  0,  &p, 

(6a) 

e;< 

~l)  =  elm, 

(8b) 

A,( 

—  1 )  =  0,  /#0, 

(8c) 

r,- 

,  (  —  1 )  —  0,  i-PPt 

(8d) 

time  update  (n>0): 

c o(«)  = 

r0(n)  =x(/i), 

(8e) 

E5(it)  = 

=  E5(/t)  =  (1  -aLSLjEJfn 

-l)+x(«) 

xw(/t), 

(8f) 

7-  l(«)  : 

=  0, 

(8g) 

el  i(n) 

=  d(»i). 

(8h) 

order  update  (/  =  0,1 

A,  («)  = 

(1  -afLSL)A,(n  -  1) 

e,_i(n)r f_,(i  -  1) 

'V0, 

(8i) 

i  - 1) 

Kf(«)  = 

=  ▲frCn)Et- !',(«),  /#0, 

(8j) 

k;(«)  = 

A,(n)E,Vi  («  -  1),  'V0, 

(8k) 

«,(«)  = 

e,_  i  (n)  +  Kf (n )r, _  ,  (n  — 

l),  'VO, 

(81) 

r  ,(«)=: 

r,_  i  (n  —  1)  +  Kf(n)e,_ ,  (n),  jy 0, 

(8m) 

Ef(»t)  = 

:  Ef_,  (n)  —  A;(n)E,T'i  ("  ■ 

-1)A"(«), 

'V0, 

(8n) 

e;(«)  = 

E,r_ ,  (u  —  1 )  — 

i  (n)A,  («), 

'V0, 

(8o) 

Y.-iln) 

=  Z.-2(«)  +<f_,(n)E;_, 

(n)r ,_,(«). 

'V0, 

(8p) 

A  ?(/!)  = 

(1  -aL sL)Af(n  -  1) 

-  [ef_,  (n)rff(»)]/[l 

F.-i(«)]. 

(8q) 

K?(n)  = 

Af(/i)E  r'(n). 

(8r) 

ef(rt)  = 

ef_ ,  (n)  +  Kf(/i)r, (n). 

(8s) 

3.  On  the  choice  of  adaptation  coefficient 

If  WSS  hoids  for  all  time,  then  a ^  =  0  is  the  optimum 
choice  and  the  time-invariant  Wiener-Hopf  solution  is  ob¬ 
tained.  In  cases  of  practical  interest,  local  WSS  is  a  more 
realistic  assumption  and  the  effective  “averaging  interval”  is 
dictated  by  the  time  scale  of  variation  of  the  signal  statistics. 
Therefore,  the  lower  limit  of  the  adaptation  coefficient  aLSL 
( or  the  upper  limit  of  the  allowable  averaging  interval  >  is  set 
by  the  time  scale  over  which  local  stationary  can  be  as¬ 
sumed.  In  the  gradient  descent  algorithms  (LMS.GRL)  the 
upper  bound  of  the  adaptation  coefficient  is  imposed  by  nus- 
adjustmen:  noise,  caused  by  excessively  large  steps  m  the 
giadient  search  process.  In  the  case  of  the  exact  determinis¬ 
tic  least-squares  solution  (LSL),  there  is  no  misadjustment 
error.  The  upper  bound  of  is  strictly  a  question  of  statis¬ 

tical  reliability. 

Insight  may  be  gained  by  considering  the  variance  of  the 
following  cross-correlation  estimate  between  two  band¬ 
width-limited  white  noise  processes  x(n )  and y{n ) : 

=  Yx(i)yU  +  k),  (9) 

,to 

where 

Ncir  =  —  l/ln(l  —  a  lsl  )• 

The  variance  of  the  estimate  is  given  by  (see  Bendat  and 
Ptersol,  1971) 

var[RXJ,(fc)  1  **  (l/A^r/^  (0)^(0)  +  *;!,(*)]. 

(10) 

Equation  ( 10)  indicates  that  large,  spurious  cross-cor¬ 
relation  values  can  be  generated  between  two  uncorrelated 
processes  [Rxy(k)  =0]  for  sufficiently  small  Ntn  and/or 
large  autocorrelations  within  the  two  processes.  Although 
direct  comparison  between  this  simple  result  and  the  behav¬ 
ior  of  the  lattice  parameters  for  arbitrary  signals  is  not  war¬ 
ranted,  the  LSL  joint-process  filter  has  been  observed  to 
obey  this  general  rule.  Specifically,  the  filter  has  shown  a 
tendency  toward  cancellation  of  the  primary  channel  when 
the  primary  and  reference  channels  consisted  of  independent 
white  noise  sequences,  for  aUL  >0.02.  This  empirically  de¬ 
rived  value  will  be  used  as  the  upper  bound  for  a^i 
throughout  this  article. 

III.  REVERBERATION  IN  THE  CONTEXT  OF  ADAPTIVE 
BEAMFORMING 

Oceanic  reverberation  is  usually  classified  as  volume, 
surface,  and  bottom  reverberation.  The  scatterers  responsi¬ 
ble  for  volume  reverberation  are  mostly  bio’  gical  in  nature 
(Clay  and  Medwin,  1977).  Inorganic  partic.es  are  insignifi¬ 
cant  contributors  and  reflections  from  sound  velocity  micro- 
structure  are  effectively  marked  by  biological  scattering 
(Kaye,  1978).  Zooplankton  such  as  copepods  are  the  domi¬ 
nant  source  of  volume  scattering  in  the  near-surface  region. 
In  deeper  water,  biological  scatterers  are  often  distributed 
within  diffuse  deep  scattering  layers  (DSL)  consisting  of 
siphonophores,  copepods,  pteropods,  euphausiids,  and  me- 
sopelagic  fish.  Surface  reverberation  is  generated  by  the  en¬ 
tire  spectrum  of  the  rough  air/sea  interface  and  is  a  function 
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of  wind  speed  and  the  transmitted  frequency  (Urick  and 
Hoover,  1956;  Chapman  and  Harris,  1962).  Specular  reflec¬ 
tions  from  normally  inclined  wave  facets  and  scattering  from 
an  isotropic  layer  of  bubbles  in  the  near  surface  have  also 
been  suggested  as  reverberation  sources  at  high  and  low 
grazing  angles,  respectively  (Medwin,  1966;  Clay  and  Med- 
win,  1964).  Bottom  reverberation  is  an  extremely  complex 
phenomenon  owing  to  the  diversity  of  ocean  floor  types,  lat¬ 
eral  inhomogeneity,  and  potential  contribution  of  subbot¬ 
tom  layers.  It  appears  that  both  particle  size  and  bottom 
relief  are  important  factors  (McKinney  and  Anderson, 
1964;  Buckley  and  Urick,  1968). 

The  received  backscattered  return  is  typically  a  compos¬ 
ite  process  with  two  or  more  reverberation  types  contribut¬ 
ing  at  any  given  instant.  This  is  especially  true  in  shallow- 
water  environments.  Our  objective  here  is  to  achieve 
selective  reverberation  cancellation,  whereby  the  desired  re¬ 
verberation  component  is  preserved  and  the  rest  suppressed 
or  eliminated.  In  this  context,  the  spatial  correlation  charac¬ 
teristics  of  the  various  reverberation  components  are  criti¬ 
cally  important.  Certain  theoretical  predictions  of  the  point¬ 
scattering  reverberation  model  provide  valuable  insight  into 
the  spatial  correlation  properties  of  surface,  volume,  and 
bottom  reverberation,  and  point  to  an  adaptive  system  con¬ 
figuration  conducive  to  selective  reverberation  cancellation. 

A.  The  point-scattering  model  of  reverberation 

According  to  the  point-scattering  model  developed  by 
Faure  (1964),  Ol’shevskii  (1967),  and  Middleton  (1967a, 
b,  1972a,b),  reverberation  is  described  as  a  random  process 
constructed  by  a  linear  superposition  of  the  individual  ech¬ 
oes  emanating  from  a  large  number  of  point  reflectors  locat¬ 
ed  independently  in  a  homogeneous  medium.  For  a  trans¬ 
mitted  signal  j(t),  the  backscattered  signal  is  represented  by 

r{t)  =  jr  a1G’,(f)s[a,(/-r,)].  OD 

i=  i 

The  returning  echo  from  each  scatterer,  indexed  by  is  de¬ 
layed  by  t, ,  the  two-way  travel  time  between  the  source  and 
the  ith  scatterer.  The  a,  and  a,  are  random  variables  repre¬ 
senting  distributions  of  acoustic  cross  section  and  Doppler 
factor  of  the  scatterers,  respectively.  For  a  monostatic  sonar, 

G,(t)  =  gB?(r,)F(t), 

where  g  is  a  system  gain  factor,  B]  is  the  two-way  beam 
pattern  sensitivity  in  the  direction  of  the  /th  scatterer,  and 
F(t)  represents  the  two-way  propagation  loss. 

Assuming  that  for  a  sufficiently  large  scattering  region 
the  mean  scatterer  density  is  constant,  the  number  of  scat¬ 
terers  n  contributing  to  the  reverberation  return  at  any  given 
time  may  be  described  as  a  Poisson  random  process.  Based 
on  this  assumption,  the  second-order  statistics  of  reverbera¬ 
tion  can  be  calculated  for  general  experimental  geometries, 
beam  patterns,  and  transmitted  signals. 

B.  Spatial  correlation  properties  of  reverberation 

Spatial  correlation  expresses  mathematically  the  fact 
that  returns  arriving  from  widely  separated  angular  direc¬ 
tions  tend  to  add  out  of  phase  at  spatially  separated  receiv- 


(b)  x' 


FIG.  4.  Spatial  correlation  geometry,  (a)  Volume  reverberation,  (b) 
boundary  reverberation  (4>0  =  v/1  for  vertical  separation). 


ers.  The  geometry  for  the  volume  backscattered  return  is 
depicted  in  Fig.  4(a).  The  two  receivers  R ,  and  R  2  are  a 
distance  s  apart  and  vertically  aligned.  Volume  reverbera¬ 
tion  arrives  from  a  range  of  elevation  and  azimuthal  angles 
determined  by  the  directional  characteristics  of  the  trans¬ 
mitting  and  receiving  arrays.  For  omnidirectional  transmis¬ 
sion  and  reception  and  narrow-band  (quasiharmonic) 
transmitted  signals,  the  spatial  correlation  coefficient  of  vol¬ 
ume  reverberation  is  (Ol’shevskii,  1967) 

Rv(s)  =sin(ks)/ks,  (12a) 

with 

k  =  27 t/A.  =  2  irf /c, 

where  /  is  the  transmitted  frequency  and  c  is  the  speed  of 
sound.  Thus  the  spatial  correlation  coefficient  of  volume  re¬ 
verberation  is  a  decaying  function  and  will  tend  to  zero  as  the 
transmitted  frequency  and/or  receiver  separation  increase. 

The  geometry  for  the  boundary  return  is  described  in 
Fig.  4(b).  The  two  receivers  lie  in  the  xz  plane  and  the  line 
connecting  them  forms  an  angle  <b0  with  the  x  axis.  In  this 
case,  for  omnidirectional  transmission  and  reception  and  for 
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ranges  much  larger  ihan  the  receiver  separation,  the  spatial 
correlation  coefficient  is  given  by 

RB(s)  =J0(ks  cos  4>0),  02b) 

where  J0  is  the  modified  Bessel  function  of  the  first  kind.  For 
d>0  =  0  (horizontal  separation),  the  result  is  again  a  decay¬ 
ing,  oscillating  function.  However,  for  <b0  =  tt/2  (vertical 
separation),  we  have 

«fl(j)=/o(0)  =  l. 

Therefore,  no  loss  of  spatial  correlation  is  suffered  for  verti¬ 
cal  separation  of  the  receivers.  This  can  be  explained  intu¬ 
itively  in  terms  of  the  lack  of  vertical  extent  of  the  insonified 
surface  patch.  The  above  relations  are  similar  to  some  earlier 
results  on  the  correlation  properties  of  spatially  distributed 
ambient  noise  fields  (Cron  and  Sherman,  1962;  Jacobson, 
1962).  For  directional  transmission  and  reception,  Jackson 
and  Moravan  ( 1984)  show  that  the  vertical  correlation  of 
volume  reverberation  will  be  zero  for  receiver  separations 
larger  than  the  sum  of  the  dimensions  of  the  transmitting 
and  receiving  array.  On  the  other  hand,  surface  and  bottom 
reverberation  will  be  highly  correlated  between  the  same 
vertically  separated  receivers,  at  least  for  ranges  large  rela¬ 
tive  to  the  distance  of  receiver  separation. 

Experimental  evidence  exists  supporting  these  model 
predictions.  For  instance,  Urick  and  Lund  (1964)  have 
shown  varying  degrees  of  correlation  between  vertically  sep¬ 
arated  receivers  for  different  reverberation  types,  with 
boundary  returns  displaying  decidedly  higher  values  than 
the  volume  component.  The  same  investigators  (Urick  and 
Lund,  1970a,b)  have  shown  that  the  spatial  correlation  of 
composite  reverberation  in  a  shallow-water  environment  is 
substantially  higher  in  the  vertical  than  in  the  horizontal. 
This  was  confirmed  in  a  more  recent  experiment  involving 
pulsed,  high-frequency  sonar  and  directional  receivers  ( Wil¬ 
son  and  Frazer,  1983). 

The  predicted  disparities  in  vertical  correlation  form  a 
natural  premise  for  the  application  of  the  ANC  principle, 
whereby  the  correlated  component  between  the  primary  and 
reference  channels  is  canceled.  If  the  primary  and  reference 
channels  consist  of  the  outputs  of  two  vertically  separated 
receivers  in  a  composite  reverberation  field,  one  would  ex¬ 
pect  boundary  reverberation  (the  correlated  component)  to 
be  canceled  and  volume  reverberation  (the  uncorrelated 
component)  to  be  preserved  at  the  filter  output.  Therefore,  a 
vertical  array  would  be  conducive  to  removing  boundary 
reverberation  while  preserving  the  volume  component.  This 
prediction  is  supported  by  the  experimental  results  present¬ 
ed  in  Sec.  IV.  It  should  be  noted  that  if  suppression  of  the 
entire  reverberation  background  is  desired,  e.e.,  in  the  pres¬ 
ence  of  a  discrete  coherent  signal,  then  a  horizontal  array  is 
in  order. 

C.  Reverberation  simulation 

revgen,  (REverberation  GENerator)  (Goddard, 
1985)  is  a  software  implementation  of  the  point-scattering 
model.  Simulated  returns  from  a  large  number  of  discrete 
scatterers,  distributed  randomly  throughout  the  volume  and 
the  boundaries,  are  summed  coherently  at  each  receiver  to 
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obtain  a  reverberation  time  series.  The  revge:  t  output  is  a 
dual  digital  data  stream  representing  the  I  (>>i  phase)  and  Q 
(quadrature)  components  of  the  complex-basebanded  re¬ 
verberation  signal.  A  backscattering  coefficient  for  each  re¬ 
verberation  type  and  the  densities  of  their  Poisson  distribu¬ 
tions  are  specified  by  the  user.  Scattering  layers,  random 
scatterer  motion,  platform  trajectories,  attenuation  and  re¬ 
flection  losses,  arbitrary  (multiple)  transmitting  and  receiv¬ 
ing  beam  patterns,  and  several  transmitted  signal  types  may 
be  specified  through  appropriate  revgen  parameters  to  cre¬ 
ate  realistic  experimental  settings.  Here,  we  use  revgen 
simulations  in  parallel  with  real  data  to  investigate  the  can¬ 
cellation  performance  of  fhe  adamive  beamforming 

scheme. 

IV.  BOUNDARY  REVERBERATION  REJECTION  IN  A 
SHALLOW-WATER  ENVIRONMENT 

In  this  section,  the  adaptive  beamforming  concept  is  ap¬ 
plied  to  real  reverberation  data  from  a  shallow-water  active 
sonar  experiment.  The  shallow-water  environment  creates 
some  unique  difficulties.  As  opposed  to  deep-water  settings 
where  the  acoustic  paths  are  few  and  identifiable,  shallow- 
water  reverberation  returns  arrive  in  rapid  succession,  im¬ 
posing  stringent  adaptation  requirements  on  the  algorithm 
and  complicating  the  performance  evaluation  process.  Vol¬ 
ume  reverberation  is  assumed  to  be  the  component  of  inter¬ 
est,  with  boundary  reverberation  the  interference.  This  sce¬ 
nario  may  arise,  for  instance,  in  a  fisheries  sonar  application. 
The  performance  of  a  horizontally  directed  sonar  deployed 
from  a  surface  vessel  is  often  limited  by  surface  reverbera¬ 
tion.  The  general  method  used  is  adaptive  beamforming  with 
spatial  constraints  and  the  LSL  algorithm  is  used  exclusively 
as  the  central  adaptive  processor.  Given  the  distinctly  differ¬ 
ent  scattering  mechanisms  responsible  for  the  three  rever¬ 
beration  types,  the  independence  assumptions  needed  for 
ANC  are  readily  satisfied.  No  further  assumptions  need  to 
be  made  about  the  relative  intensity  and  duration  of  signal 
and  interference  and  no  modifications  of  the  algorithm  arc 
necessary  in  order  to  impose  the  constraints. 

A.  Experiment  description 

The  experiment  was  conducted  in  the  deepest  part  of 
Dabob  Bay,  Washington  (maximum  depth  —200  m).  The 
physical,  biological,  and  acoustical  properties  of  Dabob  Bay 
have  been  studied  intensively  (Helton,  1976).  Dabob  Bay  is 
characterized  by  large  fluctuations  of  temperature  and  salin¬ 
ity.  In  addition  to  the  usual  seasonal  vanations  in  the  top 

—  30  m,  cellular  inhomogeneities  exist  in  the  deeper  regions, 
caused  by  intrusions  of  Pacific  Ocean  water.  These  transient 
water  masses  cause  wide  fluctuations  in  the  sound  velocity 
structure  and  create  the  need  for  up-to-date  sound  velocity 
profiling  during  acoustic  experiments.  Coring  device  sam¬ 
pling  (Bums,  1962)  has  shown  that  the  deep  regions  of  Da¬ 
bob  Bay  are  covered  by  silt  and  mud  with  some  sand  and  fine 
gravel.  A  maximum  bottom  backscattering  coefficient  of 

—  20  dB  has  been  measured.  Correspondingly,  the  surface 
backscattering  coefficient  was  found  to  be  in  the  range  of 

—  35  to  —  25  dB.  Intensive  volume  sampling  and  volume 

Dimitn  Alexandrov.  Boundary  reverberation  rejection  1 281 


p 


FIG.  5.  Shallow-water  experimental  geometry'  The  surface,  volume,  and  bottom  contribute  to  the  backscattered  return. 


scattering  measurements  over  a  period  of  2  years  have  re¬ 
vealed  multiple  scattering  layers  with  substantial  seasonal 
variations  (Anderson,  1981 ).  The  maximum  measured  vol¬ 
ume  scattering  coefficient  is - 60  dB. 

The  overall  experimental  geometry  is  depicted  in  Fig.  5. 
The  sonar  system  moved  with  constant  velocity  along 
straight  line  paths.  It  carried  a  transducer  array  whose  ele¬ 
ments  were  combined  by  a  programmable  conventional 
beamformer  capable  of  producing  several  transmit/receive 
beam  sets.  During  one  of  the  experimental  runs,  a  “wide” 
(60°X60°)  transmitting  beam  was  used  and  9  row  sums  of 
the  receiving  elements  were  individually  available,  thus 
forming  a  9-element  vertical  array.  Acc <■  ding  to  the  results 
of  Sec.  Ill  this  arrangement  has  good  potential  for  selective 
reverberation  cancellation.  We,  therefore,  focus  our  atten¬ 
tion  to  this  particular  data  set.  The  sound  velocity  profile 
effective  for  this  run  and  the  corresponding  ray  tracing  dia¬ 
gram  are  shown  in  Fig.  6.  The  sonar  system  transmitted  1 80- 
ms  pulses  from  a  depth  of  30  m  and  recorded  digitally  the 
complex-basebanded  reverberation  return.  In  addition,  a 
narrow-band  calibration  signal  is  present,  arriving  broadside 
to  the  array,  which  is  valuable  in  assessing  the  cancellation 
results.  A  revgen  simulation  preserving  the  experimental 
geometry  and  all  known  characteristics  of  the  sonar  system 
was  carried  out.  The  receiving  ( row )  elements  were  approxi¬ 
mated  by  rectangular  pistons  at  A  /2  spacing  and  a  constant 
sound  velocity  profile  was  assumed. 

B.  Data  description 

The  primary  beam  was  constructed  by  adding  the  out¬ 
puts  of  the  nine  row  elements.  The  simulated  primary  beam 
is  shown  in  Fig.  7(a).  In  addition,  eight  constrained  refer- 
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ence  beams  were  created  by  pairwise  subtraction  among  the 
nine  rows,  in  the  manner  depicted  in  Fig.  1(b).  The  beam  of 
a  single  constrained  reference  beam,  characterized  by  a  cen¬ 
tral  null,  is  the  dotted  curve  in  Fig.  7(b) .  The  real  and  simu¬ 
lated  (composite)  reverberation  signals  entering  the  pri¬ 
mary  beam  are  displayed  in  parallel  in  Fig.  8  The  first  plot  is 
a  range  Doppler  map  (RDM)  display  intended  to  show  the 
time-varying  character  of  the  data  in  the  frequency  domain. 
It  was  constructed  by  taking  successive  128-point  FFTs  of 
windowed  (Kaiser-Bessel,  aKB  =  2.5)  data  segments  with 
50%  overlap.  The  dB  magnitude  of  each  transform  corre¬ 
sponds  to  a  single  line  in  the  pseud i  3-D  plot.  The  main 
reverberation  “ridge”  occupies  the  near-dc  region,  following 
complex  basebanding.  The  second  plot,  an  RDM  with 
~90%  overlap,  provides  a  more  detailed  view  of  the  near 
ranges  where  boundary  reverberation  enters  through  the  si- 
delobes.  The  data  have  been  frequency  shifted  to  compensate 
for  the  Doppler  shift  due  to  sonar  motion.  Clearly  visible  is 
the  typical  “hook”  pattern  of  negative  Doppler  features  in 
the  initial  —0.5  s  of  the  ping.  This  pattern  is  entirely  due  to 
high-angle  boundary  returns  which  are  imparted  a  lesser  net 
Doppler  shift  by  the  velocity  of  the  sonar.  With  increasing 
range,  as  the  lower  sidelobes  and  eventually  the  mainlobe 
become  the  boundary  reverberation  outlets,  the  pattern 
merges  into  the  main  ridge.  This  artifact  is  useful  for  identi¬ 
fying  sidelobe  returns  and  establishing  their  cancellation  fol¬ 
lowing  adaptive  processing. 

Although  the  real  and  simulated  data  have  similar  over¬ 
all  structure,  several  differences  are  apparent.  The  real  data 
are  characterized  by  a  somewhat  wider  hook  which  may  be 
an  indication  that  the  real  beam  pattern  has  an  appreciable 
sidelobe  response  at  a  higher  angle  than  the  simulated  beam 
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FIG.  6.  Velocity  profile  and  ray  trace  diagram 
(a)  Real  data,  (b)  revgen  simulation.  Note 
that,  in  the  real  data  case,  a  shadow  zone  is 
formed  in  the  near  surface  and  no  surface  rever¬ 
beration  is  produced  at  the  far  ranges.  In  the 
simulated  case,  on  the  other  hand,  surface  re¬ 
verberation  continues  to  be  produced  at  the  far 
ranges. 
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FIG.  7.  Beam  patterns,  (a)  Primary  beam  pattern,  produced  by  adding  all 
nine  rows;  (b)  single  vertical  CT  bearing  cuts  of  the  primary  beam  (solid 
line)  and  of  a  single  constrained  beam  pattern  produced  by  subtracting  two 
rows  (dotted  line) 


or  that  the  actual  velocity  of  the  sonar  was  underestimated  in 
the  simulation.  In  addition,  the  simulated  reverberation  is 
better  sustained  at  long  ranges.  Two  reasons  may  be  ci  ed  for 
this.  First,  no  decaying  function  of  grazing  angle  was  intro¬ 
duced  in  the  simulation  for  boundary  reverberation.  This 
decay  with  grazing  angle  is  a  well-documented  fact  (Urick, 
1975 ) .  The  second  reason  may  be  inferred  from  the  ray  trac¬ 
ing  diagram  (Fig.  6)  and  the  time  of  arrival  versus  angle  of 
arrival  plot  (Fig.  9).  Because  of  the  well-developed  thermo- 
cline,  a  shadow  zone  is  formed  in  the  near  surface  past  a 
range  of  —500  m.  The  surface  contribution  disappears  ab¬ 
ruptly  at  that  point.  In  the  simulated  case,  an  isovelocity 
profile  is  used  and  the  surface  continues  to  contribute  at  long 
ranges.  Finally,  a  high-Doppler  calibration  signal  can  be  ob¬ 
served  in  the  real  data  RDM  at  f  — 1.5  s.  It  arrives  broadstue 
to  the  array.  This  signal  was  not  included  in  the  simulation. 

C.  ANC  with  a  single  constrained  reference  element 

A  single  constrained  reference  element,  created  by  sub¬ 
tracting  one  row  clement  from  another,  has  a  beam  pattern 
[  Fig.  7  ( b )  ]  which  at  first  glance  appears  to  be  well  suited  for 
our  purposes.  Excluding  its  broadside  null,  it  is  nearly  omni¬ 
directional  in  the  vertical,  thus  covering  all  possible  direc¬ 
tions  of  arrival  for  boundary  reverberation.  One  may  expect 
it  to  be  effective  in  canceling  the  sidelobe-bom  boundary 
returns  and  to  some  extent  mainlobe  interference,  while  pro¬ 
tecting  much  of  the  volume  signal  in  its  null. 

The  data  were  processed  in  this  manner,  using  the  scalar 
LSL  with p  —  3  and  =  0.02.  The  single  reference  chan- 
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FIG.  8.  Shallow-water  data  Real  (left-hand  side)  and  revgen  (nght-hand  side),  (a)  Standard  RDM,  (b)  expanded  near-range  plot  (note  the  typical  hook 
pattern  caused  by  the  high-angle  boundary  reverberation  returns),  and  (c)  time  series  plot. 


nel  consisted  of  rows  5-6.  The  results  for  both  the  real  and 
simulated  data  are  displayed  in  Fig.  10.  In  addition  to  the 
RDMs  and  the  time  series  plots,  a  cancellation  curve  is  in¬ 
cluded  which  is  simply  the  dB  magnitude  of  the  ratio  of 
primary  signal  power  over  filter  output  power.  Cancellation 
is  generally  low  in  the  near  ranges.  In  addition,  the  simulated 
case  displays  a  significant  increase  in  background  noise  level 


for  f<0.5  s.  Moderate  cancellation  was  achieved  in  both 
cases  for  intermediate  ranges  ( 1 .0<  ?<  2.0 ) ,  while  the  behav¬ 
ior  in  the  far  ranges  is  distinctly  different  for  real  and  simu¬ 
lated  data.  For  th  latter,  cancellation  continues  to  improve 
with  range,  while  little  or  no  cancellation  is  observed  for  the 
former. 

The  low  overall  cancellation  in  the  near  ranges  is  a  man- 
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FIG.  9  Angle  versus  time  of  arrival  (a)  Real,  (b)  revgen  data  These 
plots  were  intended  to  illustrate  the  origin  of  the  "uncorrelated  noise”  ef¬ 
fect.  In  the  early  ranges,  the  constrained  reference  beam  (dotted  line)  re¬ 
ceives  reverberation  from  angular  directions  where  the  primary  beam  ( solid 
line)  has  a  null  response. 


ifestation  of  the  ANC  theoretical  prediction  regarding  the 
effect  of  uncorrelated  noise  components  in  either  the  pri¬ 
mary  or  the  reference  channel  ( Widrow  et  al.,  1975).  This 
effect  comes  into  play  here  as  follows:  The  transmit  beam 
insonifies  the  entire  angular  region  from  —  60°  to  +  60°. 
The  reference  beam,  being  nearly  omnidirectional  in  the  ver¬ 
tical,  receives  reverberation  from  the  entire  range  with  the 
exception  of  the  near  broadside.  The  main  beam,  on  the  oth¬ 
er  hand,  has  a  definite  sidelobe  response  which  operates  on 
the  high-angle  returns.  Each  time  boundary  reverberation 
arrives  through  the  nulls  between  the  sidelobes  it  is  excluded 
from  the  primary  channel.  The  same  components  enter  the 
reference  channels  without  being  attenuated  and  cause  the 
uncorrelated-to-correlated  noise  power  ratio  to  increase,  re¬ 
sulting  in  reduced  signal-to-noise  ratio  at  the  filter  output. 
This  explains  the  observed  rise  of  the  noise  “floor”  [Fig. 
10(a)  ].  The  somewhat  better  cancellation  apparent  for  the 
high-angle  (  +  60°)  returns  in  the  real  versus  the  simulated 
data  may  be  an  artifact  of  the  RDM  displays  caused  by  dy¬ 
namic  range  differences,  unavoidable  in  view  of  the  necessar¬ 
ily  incomplete  model  of  reality  which  is  used  in  the  simula¬ 
tion.  The  particular  sidelobe  structure  of  the  real  main  beam, 
which  may  be  substantially  different  from  the  simulated 
beam,  could  also  explain  this  difference. 

The  difference  in  far-range  cancellation  performance  is 
due  to  the  grazing  angle  effect.  In  the  real  case,  there  is  sim- 
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ply  little  boundary  revet  beration  left  in  the  far  ranges  and 
therefore  the  cancellation  potential  is  low.  This  effect  was 
not  included  in  the  simulation  and  as  a  result  significant 
cancellation  was  achieved  there.  Interestingly,  there  is  a 
clear  trend  of  improving  cancellation  performance  with 
range,  which  is  consistent  with  another  ANC  theoretical 
prediction.  Specifically,  the  output  noise  spectrum  was 
found  to  be  proportional  to  the  signal-to-correlated  noise 
ratio  at  the  primary  input  (Widrow  et  al.,  1975).  In  our 
experimental  geometry,  this  ratio  shou'd  decrease  with 
range  as  the  boundary  returns  begin  to  enter  through  the 
mainlobe  instead  of  the  sidelobes.  Therefore,  one  may  expect 
the  output  noise  power  to  decrease  and  cancellation  to  im¬ 
prove  with  range. 

Two  additional  observations  can  be  made  concerning 
the  real  data  results.  First,  the  high-Doppler  “control”  sig¬ 
nal  was  virtually  unaffected  by  this  processing.  This -is  an 
indication  that  the  constraining  scheme  is  effective.  Second, 
a  large  zero-Doppler  feature  was  revealed  [marked  VI  m 
Fig.  10(c)  ]  at  t  ~  1.0  s.  This  is  probably  a  volume  reverbera¬ 
tion  feature,  the  origin  of  which  will  be  discussed  later. 

1.  Remarks 

(i)  Processing  with  i  single  constrained  reference  beam 
would  be  more  effective  in  the  sidelobe  cancellation  mode  if 
the  transmit  beam  was  identical  to  the  main  receive  beam.  In 
that  case,  the  “uncorrelated  noise”  effect  would  not  exist  and 
near-range  cancellation  would  improve. 

(ii)  Regarding  the  constrained  reference  cancellation 
beam,  a  choice  must  be  made  about  the  width  of  the  broad¬ 
side  null.  It  would  be  desirable  to  have  a  wide  null  initially, 
when  interference  arrives  through  the  sidelobes  and  a  pro¬ 
gressively  narrower  null  as  the  main  beam  begins  to  intersect 
the  boundaries. 

(iii)  An  alternative  explanation  for  the  improved  can¬ 
cellation  observed  in  the  far  ranges  is  in  terms  of  the  in¬ 
creased  spatial  coherence  of  boundary  reverberation  pre¬ 
dicted  by  the  point-scattering  model. 

D.  An  application  of  constrained  adaptive  beamforming 

The  multichannel  joint-process  filter  can  compensate 
for  noise  components  not  present  in  the  primary  channel, 
which  are  correlated  between  the  reference  channels.  There¬ 
fore,  if  the  multichannel  LSL  filter  were  to  be  used  in  the 
present  case,  with  the  eight  constrained  reference  elements 
comprising  the  vector  reference  input,  one  may  expect  it  to 
be  more  effective  than  the  single-channel  filter  by  alleviating 
the  near-range  “uncorrelated  noise”  effect.  Equivalently, 
the  multichannel  adaptive  algorithm  should,  m  theory,  form 
an  adaptive  reference  (cancellation)  beam,  with  the  tenden¬ 
cy  to  emulate  the  sidelobe  features  of  the  main  beam.  This 
adaptive  reference  beam  will  still  have  the  central  null  ac¬ 
cording  to  the  beam  pattern  product  theorem  (Urick,  1975 ) 
and,  therefore,  the  basic  constraining  scheme  will  still  be  in 
effect. 

The  results  for  the  same  ping  processed  through  the  8- 
channel  LSL  filter  are  displayed  in  Fig.  1 1 .  Overall,  cancella¬ 
tion  is  substantially  improved  over  the  single  reference  ele- 
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FIG  10.  ANC  with  single  constrained  reference  element  Real  (left-hand  side)  and  REVGEN  simulation  (nght-handside).  (a)  Standard  RDM,  (b)  «tpand- 
ed  RDM.  Substantial  reverberation  cancellation  is  seen  for  both  real  and  simulated  data  in  the  far  ranges,  but  near-range  performance  is  marginal.  Note  the 
increase  in  noise  level  m  the  simulated  results  (c)  Output  time  senes  and  (d)  dB  cancellation  curve.  A  volume  “signal”  (VI )  is  exposed  in  the  real  data 
output.  Cancellation  is  seen  to  improve  with  range  for  simulated  data  and  to  level  out  for  real  data. 


ment  case.  As  evidenced  by  the  virtual  elimination  of  the 
hook  pattern  from  the  RDM  plots,  the  sidelobe  boundary 
interference  has  been  successfully  removed.  A  dramatic  im¬ 
provement  is  observed  in  far-range  cancellation.  As  a  result. 


in  the  real  data  case,  volume  feature  VI  is  even  more  clearly 
in  evidence  at  r  —  1 .0  s  and  an  additional  feature  (V2)  was 
exposed  at  t~  1.6  s. 

In  order  to  examine  more  closely  the  cancellation  per- 
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FIG.  1 1.  Adaptive  beamforming  with  eight  constrained  reference  elements.  Real  (left-hand  side)  and  revgen  simulation  ( right-hand  side)  (a)  Standard 
RDM;  (b)  expanded  RDM.  Cancellation  improves  dramatically  for  both  real  and  simulated  data  in  all  ranges,  (c)  Output  time  senes  and  dB  cancellation 
curve.  An  additional  volume  “signal”  (V2)  is  revealed. 


formance  of  this  arrangement,  a  separate  simulation  run  was 
performed  with  only  volume  scatterers  present.  Thus  the 
desired  “signal”  was  isolated.  Figure  12(a)  shows  that  the 
multichannel  constrained  adaptive  beamformer  is  generally 
very  successful  toward  achieving  our  objective,  especially  for 
t>  1-0  s.  Insight  may  be  gained  on  the  adaptive  operation 
performed  by  this  structure  by  considering  the  adaptive  can- 
ce’lation  beam  it  creates.  Vertical  “cuts”  of  this  beam,  at 
zero  bearing,  were  calculated  from  each  set  of  eight  cross¬ 


channel  reflection  coefficients  at  ~0.024-s  intervals  and  are 
displayed  in  a  “waterfall”  plot  [Fig.  12(b)  ].  Note  that  the 
adaptive  cancellation  beam  still  has  the  central  adaptive 
null,  as  expected.  Moreover,  it  strives  to  approximate  the 
sidelobe  stiucture  of  the  main  beam  in  the  early  ranges  of 
high-angle  boundary  returns.  At  longer  ranges  ((>1.0  s), 
when  the  sidelobe  activity  subsides,  the  adaptive  cancella¬ 
tion  beam  concentrates  on  the  mainlobe  boundary  returns  to 
the  extent  allowed  by  the  broadside  null. 
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FIG  12.  Selective  reverberation  cancellation  performance,  (a)  Composite  reverberation  (top  solid  line),  pure  volume  component  (bottom  solid  line),  and 
output  of  8-channel  adaptive  beamformer  (dotted  line)  With  the  exception  of  a  near-range  region  (0.5-1  0  s),  selective  reverberation  cancellation  is 
achieved  (b)  “Waterfall”  plot  of  the  adaptive  reference  (cancellation)  beam  produced  by  the  algorithm.  The  adaptive  beam  strives  to  emulate  the  sidelobe 
structure  of  the  primary  beam  in  the  near  ranges 
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FIG.  13.  revcen  symmetric  scenario  (a)  Cancellation  performance  improves  in  the  near  ranges  and  (b)  a  symmetric  adaptive  cancellation  beam  is 
produced  Vertical  "cuts"  of the  adaptive  beam  pattern  at  (c)  r  =  0.5s,  (d)  r=  1.0s, and  (e)  r  =  2.0s  show  three  stages  of the  adaptation  process.  Initially, 
the  reference  beam  (dotted  line)  emulates  the  sidelobe  structure  of  the  primary  beam  (solid  line),  thus  eliminating  the  near-range  “uncorrelated  noise” 
effect.  As  the  high-angle  sidelobe  returns  diminish  and  boundary  reverberation  begins  to  enter  through  the  mainlobe,  the  reference  beam  progressively 
narrows  its  center  null  to  compensate  for  the  near-grazing  boundary  returns. 
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FIG  14.  Volume  “feature”  interpretation  The  two  volume  "signals"  ex¬ 
posed  in  the  data  by  the  adaptive  beamforming  operation  (a)  correspond 
well  in  range  with  two  caustics  resulting  from  thesharp  velocity  gradient  m 
the  near  surface.  The  ray  trace  diagram  (o)  covers  an  angular  range  of 
±  i'.  which  is  covered  by  the  null  of  the  reference  channel  Therefore,  the 
strong  volume  reverberation  returns  originating  at  the  caustic  reg-ons  are 
protected  from  cancellation  and  preserved  at  the  filter  output 


1.  Simulation  with  symmetric  experimental  geometry 

In  order  to  further  examine  the  nature  of  the  adaptive 
cancellation  beam,  the  same  processing  was  earned  out  in  a 
more  benign  symmetric  expenmental  geometry.  In  this 
revgen  simulation,  the  overall  depth  was  set  to  300  m  and 
the  sonar  was  positioned  at  150  m.  The  surface  and  bottom 
backscattering  coefficients  were  both  set  equal  to  —  30  dB. 
In  the  sequence  of  plots  presented  in  Fig.  13,  it  is  apparent 
that  the  symmetrical  arrangement  is  beneficial  to  the  cancel¬ 
lation  performance  in  the  near  ranges.  The  adaptive  cancel¬ 
lation  beam  is  itself  symmetnc  and  its  behavior  easier  to 
interpret.  In  its  sidelobe  cancellation  mode  (r  =  0.5)  it  suc¬ 
cessfully  emulates  the  sidelobe  structure  of  the  main  beam. 
As  the  mainlobe  begins  to  intersect  with  the  bo"ndaries,  the 
adaptive  filter  is  seen  to  alter  adaptively  the  effective  width 
of  the  broadside  null,  making  it  progressively  narrower  with 
increasing  range  [plots  for  t  =  1.0  and  t  —  2.0  in  Fig.  13(d) 
and  (e)  ].  This  is  intuitively  satisfying  and  explains  the  dra¬ 
matic  improvement  in  far-ranse  cancellation  performance 
achieved  by  the  adaptive  beamformer  over  single-channel 
(i.e.,  fixed  null  shape)  ANC. 
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E.  Volume  “feature  ’  interpretation 

A  number  of  hypotheses  were  considered  regarding  the 
origin  of  the  volume  features  (VI  and  V2)  exposed  by  the 
adaptive  beamforming  operation  on  the  real  data.  First,  the 
possibility  that  they  are  acoustic  echoes  from  schools  of  fish 
or  cellular  water  masses  within  the  angular  range  of  the  null 
was  debated.  This  hypothesis  was  rejected  based  on  the  pro¬ 
cessing  results  for  another  ping  preceding  our  data  by  ~  104 
s.  During  this  time,  assuming  they  were  moving  slowly  rela¬ 
tive  to  the  sonar,  one  would  expect  the  echoes  to  be  at  dis- 
cemibly  different  langes  for  the  two  pings.  This,  m  fact,  was 
not  the  case;  the  features  appeared  at  precisely  the  same 
ranges. 

Alternatively,  the  features  may  be  explained  by  the  pres¬ 
ence  of  a  strong,  well-defined  scattering  layer  at  a  shallow . 
depth.  Such  a  layer,  consisting  mainly  of  prespawning  popu¬ 
lations  of  Pacific  herring,  has  been  consistently  observed  in 
the  spring  months  (Friedl,  1970;  Anderson,  1981 ).  If,  due  to 
refraction  effects,  the  near  broadside  rays  residing  in  the  null 
region  intersect  this  layer  at  appropriate  ranges,  they  could 
possibly  account  for  the  features.  The  high-resolution  ray 
tracing  plot  created  to  provide  support  for  this  hypothesis 
exposed  additional  evidence  (Fig.  14).  Specifically,  two 
caustics  are  clearly  visible  at  ranges  which  correspond  w ell 
with  the  feature  position.  The  intense  scattering  caused  by 
'.he  presence  of  the  caustics  in  the  strong  shallow  scattering 
layer  can  probably  account  for  the  two  volume  features. 


V.  SUMMARY 

An  adaptive  beamforming  technique  with  simple  spatial 
constraints  was  considered.  The  theoretical  foundation  of 
the  joint-process  least-squares  lattice  algorithm  was  ex¬ 
posed.  The  predictions  of  the  point-scattering  model  regard¬ 
ing  the  spatial  correlation  properties  of  oceanic  reverbera¬ 
tion  were  found  to  be  valuable  in  determining  optimum 
adaptive  array  configurations.  Data  from  a  shallow-water 
active  sonar  experiment  were  used  to  demonstrate  the  poten¬ 
tial  of  removing  boundary  reverberation  while  preserving 
the  volume  component.  Several  predictions  of  Widrow  et 
al.'s  ( 1975)  ANC  model  were  found  to  have  direct  bearing 
on  our  results.  Constrained  adaptive  beamforming  imple¬ 
mented  through  the  multichannel  joint-process  least- 
squares  lattice  was  shown  to  offer  superior  performance  over 
single-channel  ANC.,  revgen  simulations  of  the  experi¬ 
ment  proved  to  be  very  useful  in  assessing  algorithm  perfor¬ 
mance  and  interpreting  the  cancellation  results. 
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Abstract — A  difital  array  of  120  acoustic  channels  and  900  m  ip  length  sour  is  are  present  [2],  [3].  This  component  may  include  re¬ 


in  bet*  contracted  to  stud;  low-frequency  (20-200  Hz)  ambient  noise 
ia  the  ocean.  The  array  nay  be  deployed  vertically  or  horizontally  from 
tie  research  platform  FLIP  and  the  array  elements  are  localized  with  a 
Ugh-freqnency  acontte  transponder  network.  This  paper  describes  the 
Snstruaeatation,  telemetry,  and  navigation  systems  of  the  array  daring  a 
vertical  deployment  I*  the  northeast  Pacific.  Preliminary  ambient  noise 
spectra  are  presented  for  various  array  depths  and  local  wind  speeds. 
Ambient  noise  in  the  frequency  band  above  100  Hz  or  below  25  Hz  in- 
creates  with  local  wind  speed.  However,  in  the  frequency  band  25-100 
Hz  ambient  noise  is  independent  of  wind  speed  and  may  be  dominated 
by  shipping  sources. 

Keywords — acoustic  array,  ambient  noise,  acoustic  navigation,  low  fre¬ 
quency. 

I.  Introduction 

A  DIGITAL  ARRAY  of  120  acoustic  channels  and  900  m 
in  length  has  been  constructed  for  the  study  of  low- 
frequency  (20-200  Hz)  ambient  noise.  A  large  aperture  ar¬ 
ray  is  required  for  high-resolution  directional  information  at 
low  frequencies.  A  well-filled  array  is  required  to  provide  low 
side-lobe  levels  for  the  study  of  ambient  noise.  This  paper  de¬ 
scribes  the  characteristics  of  a  large  aperture  linear  array  of 
hydrophones  which  may  be  deployed  vertically  or  hoitr  ‘ally 
from  the  research  platform  FLIP.  Preliminary  observations 
of  low-frequency  ambient  noise  from  a  vertical  deployment  of 
this  array  in  the  northeast  Pacific  during  September  1987  are 
described. 

Oceanic  ambient  noise  is  the  prevailing  sustained  back¬ 
ground  of  sound  in  the  ocean.  These  noise  levels  place  con¬ 
straints  on  the  operation  of  acoustic  sensors  in  the  ocean.  For 
this  reason,  it  is  valuable  to  understand  the  sources  generat¬ 
ing  the  sound,  the  absolute  levels,  and  the  spectral  shapes  of 
ambient  noise.  Low-frequency  noise  is  particularly  important 
because  of  its  low  attenuation  and  therefore  its  ability  to  prop¬ 
agate  over  long  distances.  Only  recently  has  it  been  practical 
to  investigate  the  directionality  of  ambient  noise  sources  at 
low  frequency  (less  than  100  Hz)  due  to  power,  size,  and  cost 
constraints. 

The  following  have  been  identified  as  sources  of  low- 
frequency  noise  in  the  ocean  [1].  shipping,  wind  and  waves, 
seismic  disturbances,  and  nonlinear  ocean  wave  interactions. 
In  the  frequency  band  between  20  and  100  Hz,  shipping 
is  thought  to  be  the  dominant  noise  source  where  shipping 
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verberative  paths,  perhaps  related  to  prominent  bottom  topo¬ 
graphic  features,  as  well  as  forward  scattering  and  channeling 
of  the  shipping  sources.  The  variation  of  shipping  noise  may 
depend  on  whether  the  sources  are  of  single-ship  or  multiple- 
ship  origin  and  whether  they  are  local  or  distant. 

At  frequencies  between  100  and  200  Hz,  sea  surface  noise 
generated  by  wind  and  waves  may  be  dominant.  The  reported 
contribution  by  local  wind  sources  in  this  band  is  quite  vari 
able,  ranging  from  a  difference  in  spectral  level  of  approx¬ 
imately  +  18  to  -4  dB/jiPa/  VHz  [1],  [4],  The  sources  re¬ 
sponsible  for  this  variation  may  be  identified  by  examining  the 
vertical  and  horizontal  directionality.  Distant  storms  and  noise 
generated  at  the  edges  of  the  ocean  due  to  waves  breaking  on 
cliffs,  rocks,  or  beaches  as  well  as  the  nonlinear  generation 
due  to  interference  of  incoming  and  reflected  coastal  swell 
may  produce  an  azimuthally  nonuniform  contribution;  local 
wind  should  induce  a  vertically  variable  contribution. 

The  unique  capabilities  of  the  array  described  in  this  pa¬ 
per  allow  its  deployment  as  a  high-performance  vertical  or 
horizontal  array.  Arrays  previously  used  to  measure  low- 
frequency  vertical  directionality  are  listed  in  Table  I  to  fa¬ 
cilitate  comparison.  The  number  of  elements  (120)  and  large 
aperture  (900  m)  of  our  array  is  substantially  greater  tnan 
previously  reported  vertical  arrays.  This  large  aperture  will 
allow  higher  resolution  vertical  directionality  than  was  previ¬ 
ously  available.  Horizontal  directionality  is  usually  measured 
using  large  aperture  towed  arrays.  Our  array  can  be  moored 
horizontally  because  of  its  neutrally  buoyant  design.  When  op¬ 
erated  in  the  horizontal  configuration,  the  flow  noise  affecting 
our  array  is  significantly  less  than  for  towed  arrays,  leading 
to  improved  array  performance. 

n.  Array  Description 

This  section  describes  the  array  electrical  and  mechanical 
design.  The  array  has  a  modular  design,  which  ricilitates  as- 
sembly  and  transportation  and  allows  for  a  variable  aperture. 
It  is  separable  into  identical  hose  sections  of  ran  elementsieach, 
joined  together  by  in-iine  interchangeable  pressure  cases.  rach 
of  the  ten  elements  consists  of  two  hydrophones,  a  preamp  a 
filter,  and  a  line  drive  submerged  in  insulating  no'-'T'na  oil. 
The  oil-filled  hoses  are  nenrraily  buoyant  in  seawater,  nec¬ 
essary  for  horizontal  deployments.  The  interelement  spacing 
is  ';.5  m  and  the  elements  are  secured  within  the  2.54-  cm 
diameter  uretfaene  hose  by  a  kevlar  line  which  is  terminated 
near  each  end  of  die  hose  subsection.  In-line  pressure  cases 
are  located  between  each  hose  section  of  ten  hydrophones. 
The  pressure  cases  are  45  cm  long  with  a  7.6  on  outside  di- 
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TABLE  I 

COMPARISON  OF  VERTICAL  ARRAYS 


Number  of 
Elements 

Aperture 

meters 

Hydrophone 

Spacing 

Frequency 

Range 

Hz 

Deployment 
Depth  (actual) 
meters 

Deployment 

Platform 

Reference 

120 

903 

uniform 

20-200 

400-3100 

FLIP* 

this  paper  1988 

11 

34 

nested  subarrays 

62.5-1K 

200 

free  floating 

(14]  Buckingham  and  Jones  1987 

48 

115 

uniform 

0-450 

sound  channel 

FUP 

[15]  Hodgkiss  and  Fisher  1986 

27 

93 

uniform 

0-600 

sound  channel 

FUP 

[IS]  Hodgkiss  and  Fisher  1986 

31 

180 

nested  subarrays 

20-800 

300 

free  floating 

(3]  Burgess  and  Kewley  198. 

31 

310 

uniform 

45-100 

1500 

surface  ship 

[16]  Wales  and  Diachok  1982 

12 

237 

logarithmic 

<200 

300-3100 

free  floating 

[17]  Browning, et  at  1982 

20 

<560 

(variable) 

uniform 

5-400 

700-4800 

FUP 

[18]  Tyce  1982 

[19]  Anderson  1979 

40 

97  _J 

geometric 

112-1414 

4400 

anchored  to  bottom 

[20]  Axelrod, «  at  1965 

*  FLIP  (Floating  Instrument  Platform)  is  a  manned  109  meter  spar  buoy  stable  platform  operated  by  the  University  of  California,  San  Diego, 
Marine  Physical  Laboratory. 


ameter  providing  a  low-profile  cross  section.  These  pressure 
cases  house  nonpressure  tolerant  electronics  for  processing 
and  telemetering  the  hydrophone  signals.  Hydrophone  data 
are  transmitted  asynchronously  along  the  ray  to  a  telemetry 
module  near  the  FLIP  end  of  the  array.  This  module  buffers 
the  data  and  synchronously  transmits  it  through  a  double-lay 
armored  electrical  cable  to  the  surface  where  it  is  recorded  by 
the  data  acquisition  computer  (LSI- 11).  The  tension  carried 
by  the  electrical  support  cable  is  transferred  at  the  teleme¬ 
try  module  to  a  1 .5  cm  diameter  kevlar  line.  During  vertical 
deployment,  the  hoses  and  pressure  cases  are  attached  to  the 
kevhr  line,  which  has  1500  kg  at  its  bottom  to  maintain  veri- 
cality.  The  array  is  deployed  from  the  research  platfc  rm  FLIP, 
which  maintains  station  by  a  multipoint  moor.  The  array  is 
suspended  from  a  hyarographic  winch  which  allows  it  to  be 
lowered  to  a  specified  depth  below  FLIP.  Fig.  1  shows  a 
schematic  of  the  array  configuration  during  vertical  deploy¬ 
ment. 

The  coaxial  armored  uplink  electrical  cable  carries 
frequency-multiplexed  data  in  three  bands:  uplink  data,  down¬ 
link  commands,  and  dc  power.  The  spectrum  allocated  for 
the  uplink  data  is  from  100  kHz  up  to  approximately  1.5 
MHz.  The  uplink  data  rate  is  1  Mbit  per  second  encoded 
using  a  Miller  code  to  reduce  the  bandwidth  required  to  ap¬ 
proximately  500  kHz.  The  downlink  spectrum  is  from  100 
kHz  down  to  1.5  kHz.  The  downlink  data  is  encoded  on  a  20 
kHz  carrier  which  is  used  to  synchronize  the  data  sampling 
clocks,  and  a  command  synchronization  bit  sequence  is  trans¬ 
mitted  every  2  ms.  The  cable  simultaneously  carries  the  dc 
power  for  the  array.  Each  section  of  the  array  has  dc-to-dc 
converters  that  produce  5  V  at  600  mA  and  ±  15  V  at  150  mA 
to  power  the  electronics  and  hydrophone  elements.  The 
sections  are  in  series  so  that  they  use  the  same  current;  the 
voltage  necessary  for  the  complete  array  is  5.1  V  times  the 
number  of  sections.  The  power  loss  in  the  armored  uplink 
cable  is  proportional  to  the  current  and  does  not  change  as 
the  number  of  array  sections  changes. 


VERTICAL  ARRAY 
CONFIGURATION 
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A.  Uplink  Data  Stream  Fig.  1.  Vertical  deployment  configuration  of  the  low-frequency  digital 

_.  ,  ...  ,  ,  .  ,  acoustic  array  from  the  research  platform  FLIP. 

This  section  describes  the  uplink  of  acoustic  data  from  the 
array  hydrophones  to  the  topside  electronics.  The  data  stream 
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Fig  2.  Amy  Section  Block  Diagram.  The  hydrophone  signals  are  ampli¬ 
fied  and  filtered  at  the  element  except  for  hydrophones  4  and  S  which  are 
f'ltered  within  the  processor  pressure  case  to  allow  detection  of  12  kHz 
navigauon  signals.  The  low-frequency  acoustic  signals  are  multiplexed, 
amplified,  captured  by  the  sample  and  hold,  and  converted  to  digital  for¬ 
mat  before  being  transmitted  to  the  surface. 


originates  at  the  hydrophones.  It  is  amplified,  filtered,  con¬ 
verted  to  a  digital  signal,  reformatted,  and  finally  transmit¬ 
ted  to  the  surface,  as  shown  schematically  in  Fig.  2.  The  hy¬ 
drophones  are  an  Aquadyne  AQ- 1  with  sens  tivity  of  -204  dB 
re  1  \lfi Pa  and  a  capacitance  of  12  nf.  AQ-1  hydrophones 
have  been  calibrated  with  respect  to  pressure,  temperature, 
and  frequency  and  exhibit  a  well-behaved  response  over  the 
range  of  operating  conditions  (5).  There  is  a  2  to  3  dB  re  1 
Vlfi Pa  sensitivity  increase  from  low  pressure  (near  surface)  to 
high  pressure  (6000  m),  a  0.2  to  0.3  dB  re  1  V/jiPa  sensitivity 
increase  from  0°  to  22®C,  and  a  ±  0.2  dB  re  1  V/jiPa  sensitiv¬ 
ity  variation  across  a  frequency  band  of  10  to  1000  Hz.  There 
are  two  hydrophones  per  array  element  wired  in  series  for  an 
element  sensitivity  of  -198  dB  re  1  V/jzPa.  Tested  at  a  con¬ 
stant  depth  (1830  m),  the  relative  phase  of  the  hydrophones 
is  within  0.2*C;  changing  the  depth  from  just  below  the  sur¬ 
face  to  full  operating  depth  (3000  m)  induced  a  relative  phase 


change  between  hydrophones  which  was  less  than  0.6°C.  The 
hydrophone  output  is  applied  to  a  very  low  noise  FET  pream¬ 
plifier  with  40  dB  of  gain.  The  minimum  expected  acoustic 
noise  level  is  approximately  45  dB  re  1  jxPa  at  100  Hz  [6], 
and  the  electrical  noise  in  the  preamp  is  approximately  15  dB 
below  this  level.  The  output  is  filtered  by  a  six  pole  low-pass 
phase-matched  filter  with  a  comer  frequency  at  220  Hz  and 
whose  in-band  gain  is  1 .  The  preamp  has  a  low-frequency  cut¬ 
off  below  10  Hz.  A  differential  line  driver  is  used  to  transmit 
the  signal  a  distance  of  up  to  33.75  m  to  a  processor  pres¬ 
sure  case.  There  are  ten  elements  per  75  m  section  with  each 
processor  receiving  five  inputs  from  the  hose  on  either  side. 
The  two  elements  immediately  adjacent  to  a  processor  pres¬ 
sure  case  are  filtered  within  the  pressure  case  rather  than  at 
the  element  to  provide  12  kHz  acoustic  information  required 
for  navigating  the  array. 

In  the  processor  pressure  case  the  hydrophone  inputs  are 
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20  AhmAt  ScCTOnS 


Fig.  3.  Array  System  Diagram.  The  data  path  from  the  hydrophones 
through  the  telemetry  module  to  the  magnetic  tape  drive  is  shown.  Rirther 
analysis  of  the  data  is  accomplished  by  a  Sun  Microsystems  computer. 


selected  sequentially  by  a  differential  multiplexer,  converted 
from  differential  to  single-ended  signals,  amplified  by  a  pro¬ 
grammable  variable  gain  amplifier,  captured  by  a  sample- 
and-hold  circuit,  and  converted  to  12-bit  digital  form.  The 
programmable  analog  to  the  digital  converter  (A/D)  clock  is 
synchronous  with  the  20  kHz  downlink  carrier  and  the  A/D 
outputs  are  stored  in  a  \6-word  first-in,  first-out  (FIFO)  reg¬ 
ister  before  the  processor  is  interrupted  to  take  the  ten  data 
words.  Of  the  16  bits  per  word  stored  in  the  FIFO,  12  bits  are 
data,  2  bits  are  from  the  navigation  detectors,  and  2  bits  are 
hardware  status  flags.  Prior  to  low-pass  filtering,  the  signals 
from  the  hydrophones  adjacent  to  the  pressure  case  are  routed 
to  a  12  kHz  detection  circuit.  This  circuit  compares  the  signal 
level  in  a  narrow  band  receiver  at  12  kHz  to  the  level  of  a 
broadband  receiver  to  determine  the  presence  or  absence  of  a 
12  kHz  acoustic  transponder  signal.  The  1-bit  detection  from 
each  circuit  is  multiplexed  in  with  every  12-bit  hydrophone 
word  as  it  is  stored  in  the  FIFO.  The  CPU  is  a  Motorola  68000 
operating  at  4  MHz  with  4K  ROM  and  2K  RAM  available. 
The  4  MHz  clock  is  derived  from  a  16  MHz  crystal  which 
is  phase-locked  to  a  1  MHz  clock  signal  from  the  teleme¬ 
try  module.  The  software  is  interrupt  driven  with  the  highest 
priority  interrupt  responding  to  the  A/D.  The  second  highest 
priority  interrupt  transmits  the  processed  data,  and  the  lowest 
priority  services  the  downlink  synchronization  sequence.  The 
signal  data  are  loaded  from  the  FIFO  into  a  large  RAM  buffer 
and  the  CPU  then  processes  the  data  before  presenting  it  ,o  a 
Manchester  code  repeatsr/encoder  for  transmission. 

The  array  processors  are  synchronized  by  a  downlink  pat¬ 
tern  at  a  500  Hz  rate.  At  synchronization,  the  farthest  pro¬ 
cessor  from  FLIP  transmits  a  synchronization  word  and  the 
processed  data.  Each  sequential  orocessor  repeats  the  data  bit 


stream  from  the  processor  before  and  inserts  its  data,  fol¬ 
lowed  by  a  unique  ID  within  a  specified  time  window.  At  the 
FLIP  end  of  the  array  is  a  telemetry  module  which  contains  a 
Manchester  decoder,  a  FIFO,  and  a  synchronous  data  trans¬ 
mitter.  The  telemetry  module  converts  the  asynchronous  data 
along  the  array  to  a  continuous  data  stream  up  the  cable  to 
FLIP.  It  adds  frame  sync  words  and  null  data  when  the  FIFO 
is  empty  or  busy.  The  data  are  driven  up  the  cable  at  a  1  MHz 
rate  by  a  high-power  amplifier  (10  W)  to  overcome  the  cable 
attenuation  of  approximately  45  dB  re  1  V.  A  bit  stream  con¬ 
taining  data  from  all  the  array  sections  is  available  every  2  ms 
at  FLIP. 

B.  Topside  Electronics 

This  section  describes  the  topside  data  acquisition  system 
(Fig.  3).  The  uplink  data  are  presented  to  an  array  inter¬ 
face  module  where  they  are  amplified  and  filtered  to  re¬ 
move  cable  phase  distortion.  A  Decom  Systems  bit  synchro¬ 
nizer/decommutator  locks  to  the  frame  sync  word  in  the  data 
bit  stream,  and  decodes  the  data  to  present  it  to  the  data  acqui¬ 
sition  computer  (LSI-1 1).  The  decommutator  hardware  trans¬ 
fers  data  by  direct  memory  access  (DMA)  to  a  ring  buffer 
initialized  by  the  LSI- 11  data  handler  software.  An  8-word 
header  containing  a  buffer  ID  and  the  timing  sequence  infor¬ 
mation  are  prepended  to  the  ring  buffer  prior  to  the  transfer 
to  magnetic  tape.  The  tape  transfer  is  a  DMA  directly  from 
the  ring  buffer,  thus  avoiding  the  overhead  in  an  intermediate- 
user  buffer  transfer,  A  buffer  counter  tracks  the  number  of 
transfers  and  when  the  tape  is  full,  the  handler  automatically 
begins  accessing  a  second  tape  drive  without  ioss  of  data. 
Confirmation  of  the  buffer  ID,  buffer  counter,  frame  sync 
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word,  and  processor  ID  during  data  processing  verifies  the 
integrity  of  the  telemetry  system. 

C.  Downlink  Command  Stream 

This  section  describes  the  downlink  of  commands  from  the 
topside  data  acquisition  computer  to  the  array  processors.  The 
commands  are  entered  at  the  operator  terminal  or  from  a  set 
of  switches  on  the  array  interface  box  and  may  be  sent  at  any 
time  while  the  array  is  operating.  The  maximum  downlink 
command  data  rate  is  limited  to  550  baud  to  ensure  accuracy 
of  transmission.  Downlink  commands  may  be  specified  to  se¬ 
lected  processors  or  broadcast  to  all  processors.  There  are 
three  categories  of  commands  as  shown  in  the  Appendix:  di¬ 
agnostic,  control/initialization,  and  data  format.  The  diagnos¬ 
tic  commands  assist  in  localizing  errors  in  individual  pro- 
.  cessors.  Allowable  functions  are  to  test  memory,  read  spe¬ 
cific  memory  locations,  enter  and  execute  additional  machine 
codes,  alter  the  processor  position  in  the  uplink  bit  stream, 
turn  off  the  transmitter  of  the  addressed  processor  so  the  bit 
stream  is  passed  around  it,  and  full  or  partial  resets.  The 
control/initialization  commands  modify  array  operational  pa¬ 
rameters.  The  select  command  determines  the  order  in  which 
hydrophones  or  other  sensors  are  digitized,  allowing  selection 
of  other  sensors  such  as  depth  gauges.  The  scan  off  command 
causes  only  one  hydrophone  to  be  digitized  per  section.  The 
A/D  rate  command  selects  the  rate  at  which  the  A/D  will  dig¬ 
itize  the  incoming  data.  The  variable  gain  command  selects 
specific  gain  outputs  of  a  two-stage  amplifier.  Gains  between 
5  and  1000  are  obtained  by  selecting  a  gain  of  1,  5,  10,  or  20 
from  the  first  stage,  and  5,  8,  12,  or  20  from  the  second  stage. 
The  data  format  command  determines  the  data  format  to  be 
transmitted.  The  formats  available  transmit  some  combination 
of  test  data,  hydrophone  data,  navigation  data,  and  a  processor 
ID.  The  navigation  receiver  commands  select  which  naviga¬ 
tion  receiver  signal  is  digitized  when  the  choice  of  dzta  format 
restricts  the  number  of  navigation  bits  transmitted. 

D.  Navigation 

Array  acoustic  navigation  is  accomplished  by  detecting  the 
return  signals  of  near-bottom  acoustic  transponders  [7].  A 
minimum  set  of  three  ransponders  are  interrogated  from 
FLIP  at  unique  frequencies  and  their  replies  are  detected  by 
the  array.  The  time  of  arrival  of  each  reply  corresponds  io  the 
range  between  the  transponder  and  the  array  element.  A  set 
of  ranges  are  determined  for  each  array  navigation  element 
and  input  to  a  program  which  calculates  the  element  position 
[81- 

The  transponders  employed  were  developed  for  navigation 
of  the  MPL  deeptow  fish  (9).  For  navigation  of  the  array  these 
transponders  are  deployed  in  roughly  a  one  nautical  mile  equi¬ 
lateral  triangle  about  FLIP.  They  are  anchored  above  the  sea¬ 
floor  by  a  100  m  length  of  line.  To  increase  the  reliability  of 
detection  in  a  noisy  environment,  the  receive  circuitry  in  the 
transponders  (as  well  as  in  the  array)  compares  the  energy  in 
a  narrow  band  (200  Hz  effective  bandwidth)  about  the  inter¬ 
rogation  frequency  to  the  total  energy  received  in  a  passband 
of  approximately  1.5  kHz.  By  adjusting  the  bandwidth  and 


Q  of  ine  recognition  circuitry,  a  short  recognition  time  ( >  1 
ms}  and  high-noise  rejection  are  achieved.  The  interrogation 
signal  for  each  transponder  is  unique  and  upon  detection  the 
transponder  replies  with  a  3  ms  12  kHz  pulse  which  is  re¬ 
ceived  and  recorded  by  FLIP  and  the  array. 

The  transponder  locations  must  be  surveyed  to  determine 
their  relative  positions  before  using  them  to  navigate  the  ar¬ 
ray.  The  transponder  locations  are  adjusted  relative  to  an  ar- 
bitary  origin  by  an  iterative  approach  which  alternately  deter¬ 
mines  ship  positions  from  known  transponder  positions  and 
transponder  positions  from  known  ship  positions  using  a  least 
squares  approach.  The  initial  transponder  positions  are  de¬ 
rived  from  the  GPS  satellite  navigation  positions  of  the  surface 
ship  as  each  transponder  is  deployed.  A  data  set  containing 
slant  ranges,  depths,  and  initial  positions  is  input  to  the  nav¬ 
igation  program.  The  difference  between  the  calculated  and 
measured  positions  defines  the  error. 

The  horizontal  projection  (Fig.  4)  of  the  slant  range  between 
a  given  transponder  and  a  given  surface  ship  position  (or  FLIP 
or  an  array  element)  is 

Hproj(ntr,  npos)  =  \ /S2  -  D2 

where  ntr  indicates  a  particular  transponder,  npos  indicates 
a  particular  fix  or  position  of  the  surface  ship,  5  is  the  slant 
range  from  the  source  to  that  transponder,  and  D  is  the 
transponder  depth  minus  the  source  depth.  The  mean  squared 
error  is 

Emst  =  l/NY(rngxy  -  Hprojf 

where  rngxy  =  \t(xa  ~xTl)2  V (yss  -yr\ )2- 

The  horizontal  range  rngxy  is  determined  by  the  initial  X-Y 
positions  where  x5S  and  yss  represent  the  position  of  the  sur¬ 
face  ship,  and  N  is  the  number  of  transponders.  If  ihe  RMS 
errors  are  large  the  position  is  adjusted  and  the  process  iter¬ 
ates.  The  adjustment  is  calculated  using  the  method  of  steepest 
descent  to  follow  the  mean  squared  error  gradient  to  a  mini¬ 
mum.  For  known  transponder  positions,  the  perturbed  ship 
position  in  the  x  direction  is 

=  ^ss  “b  ,  yss) 

Xss  =  +  h  *  (xss  -  XT}) 

*  ( rngxy  -  Hproj  )/(rngxy  *  N) 

where  xss  is  the  surface  ship  position,  h  is  the  step  size,  and 
E'm„  is  the  derivative  of  the  error  function  with  respect  to  *. 
The  y-direction  adjustment  is  calculated  similarly.  When  the 
RMS  error  becomes  small  the  current  position  is  saved. 

Assuming  the  ship  positions  are  known  parameters,  the 
transponder  positions  are  adjusted  using  the  same  technique. 
Upon  completion  of  the  adjustment  loops,  the  RMS  error  for 
all  the  transponders  is  evaluated.  If  this  error  is  not  acceptable 
then  the  RMS  errors  associated  with  each  ship  position  are 
examined,  and  any  position  with  an  error  greater  than  a  speci¬ 
fied  value  is  deleted  and  the  process  begins  again.  This  method 
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Fig  4.  Navigation  Overview.  The  horizontal  projection  is  estimated  first 
by  using  the  slam  range  and  depths  (Hproj)  and  then  with  the  initial  X-Y 
positions  (rngxy). 


gives  the  best  transponder  locations  in  a  least  squares  sense 
which  are  consistent  with  the  available  data. 

The  array  and  FLIP  are  navigated  similarly.  A  hydrophone 
located  at  the  bottom  of  FLIP  transmits  a  series  of  transpon¬ 
der  interrogation  sequences  (TIS).  A  TIS  consists  of  four  in¬ 
terrogation  pulses  at  10  second  intervals  followed  by  a  silent 
interval.  The  first  three  pulses  are  at  the  interrogate  frequen¬ 
cies  of  the  bottom  transponders  (10,  10.5,  and  1 1  kHz).  Upon 
receiving  an  interrogation  pulse  a  bottom  transponder  replies 
with  a  3  ms  pulse  at  12  kHz.  The  fourth  TIS  pulse  is  at  12 
kHz  and  is  received  by  the  array  navigation  elements  to  in¬ 
dicate  array  depth.  The  array  therefore  receives  four  consec¬ 
utive  12  kHz  pulses,  whose  timing  indicates  the  transponder 
ranges  and  depth  beneath  FLIP.  The  array  samples  the  12 
kHz  pulses  at  an  operator  selected  rate  (typically  0.4  ms).  The 
CPU  decimates  the  data,  if  necessary,  to  provide  a  contin¬ 
uous  time  series  consistent  with  the  number  of  bits  allowed 
for  navigation.  The  interrogation  sequence  is  synchronized 
with  the  data-sampling  timebase  in  the  array  and  the  initia¬ 


tion  time  of  the  sequence  is  recorded.  The  navigation  time 
series  establishes  the  range  from  each  navigation  element  to 
each  transponder  after  removal  of  the  FLIP  to  transponder 
ranges  [10].  The  ranges  are  corrected  for  a  varying  sound 
speed  profile  by  integrating  over  the  ray  path. 

Calculation  of  the  hydrophone  position  begins  by  determin¬ 
ing  the  two  intersections  of  the  horizontal  projection  arcs.  The 
third  transponder  range  determines  which  intersection  is  used 
as  an  initial  position.  The  data  are  iterated  as  described  above 
to  reduce  the  RMS  errors.  Array  element  relative-location 
accuracies  of  a  few  meters  may  be  achieved  by  this  method. 
Examples  of  array  element  positions  from  the  September  1987 
experiment  are  displayed  in  Fig.  5,  demonstrating  the  relative 
motion.  The  symbols  represent  the  array  at  4-hour  incre¬ 
ments  spanning  a  24-hour  period.  There  is  less  than  1°  tilt 
from  vertical  across  the  900  m  array  aperture.  The  motion 
of  FLIP  (Fig.  5(a))  appears  to  be  driven  by  the  increasing 
northerly  wind  over  the  time  analyzed  (see  Fig.  7(h)).  The 
north-south  motion  of  the  array  (Fig.  5(b))  is  affected  by  the 
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Fig.  5.  Navigation  Results  A  time  senes  of  the  array  position  sampled 
approximately  every  four  hours.  The  three  plots  represent  (a)  a  plan  view 
in  X  (increasing  toward  the  east)  and  Y  (increasing  toward  the  north), 
(b)  X  versus  depth,  and  (c)  Y  versus  depth.  The  scale  in  (b)  and  (c)  is 
enhanced  in  the  vemcal  direction  by  8:1. 


motion  of  FLIP  as  well  as  by  tidal  motion.  The  east-west 
array  motion  (Fig.  5(c))  is  primarily  tidal  driven,  with  a  semi¬ 
diurnal  period. 

HI.  Experimental  Results 

Ambient  noise  data  are  presented  from  an  experiment  con¬ 
ducted  in  the  northeast  Pacific  during  the  month  of  Septem¬ 
ber  1987.  Data  were  collected  for  20  days,  approximately 
400  km  southwest  of  Monterey  in  4700  m  of  water.  The  pro¬ 
grammable  data  sampling  rate  was  selected  at  500  Hz.  The 
array  was  deployed  at  various  water  depths  spanning  the  wa¬ 
ter  column  from  400  to  3100  m.  Fig.  6  shows  the  array  depth 
superimposed  on  the  sound  speed  profile. 

Representative  hydrophone  spectra  illustrate  the  ambient 
noise  variation  with  depth  and  wind  speed  (Fig.  7).  The.  local 
wind  speed  varied  from  4  to  28  knots  during  the  experiment, 
as  shown  in  Fig.  7(h).  The  spectra  were  obtained  by  inco¬ 


herently  averaging  1 1  8192-point  FFT’s  with  50  percent  over¬ 
lap  of  a  Kaiser-Besse)  (a  =  2.5)  windowed  time  series.  The 
nearly  uniform  spectral  level  with  depth  is  consistent  with 
previous  investigation  [1],  [11].  There  is  a  distinct  difference 
in  the  spectral  level  as  wind  speed  increases  (Fig.  7(a)).  The 
spectra  at  a  nominal  depth  of  1300  m  at  low  wind  speed  (Fig. 
7(c))  and  high  wind  speed  (Fig.  7(d))  were  bin-averaged  and 
subtracted  to  provide  a  measure  of  wind  speed  dependent  vari¬ 
ations.  Significant  variation  in  ambient  noise  was  observed  for 
frequencies  above  100  Hz  and  below  25  Hz.  Above  100  Hz  at 
all  depths,  an  increasing  amplitude  and  distinct  whitening  of 
the  spectra  are  observed  at  high  wind  speeds.  The  observed 
spectral  difference  is  consistent  with  previous  observations 
which  are  diverse  m  this  frequency  band  (see  Table  II).  An¬ 
other  effect  of  the  wind  is  the  level  of  mechanical  vibration 
or  st'u.-n  which  increases  with  wind  speed.  This  source  is  im¬ 
portant  to  the  spectral  shape  below  25  Hz  at  all  wind  speeds 
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Fig.  6.  Array  deployment  depths  for  the  September  198?  sea  test.  The  three 
curves  represent  (a)  the  historical  sound  speed  profile  (dotted)  obtained 
from  the  National  Oceanographic  Data  Center  [12]  for  the  experimental 
area  and  time  of  year  (designated  area  25C),  (b)  the  sound  speed  (dashed) 
obtained  from  a  conductivity,  temperature,  and  depth  (CTD)  cast  deployed 
from  a  surface  ship  approximately  40  km  northwest  of  the  array,  and  (c) 
the  sound  speed  (solid)  calculated  with  the  Clay  and  Med  win  [13]  equation 
using  the  historical  salinity  data  and  an  expandable  bathymetric  thermistor 
(XBT)  deployed  from  FLIP. 


but  may  contribute  energy  between  25  and  50  Hz  at  higher 
wind  speeds.  The  bandwidth  of  the  strum-contaminated  noise 
increases  with  wind  speed  and  decreases  with  depth.  Prelim¬ 
inary  analysis  indicates  that  frequency-wavenumber  filtering 
is  effective  in  removing  this  noise  source  as  the  vibrational 
modes  travel  at  velocities  other  than  that  of  acoustic  energy. 
In  the  band  between  20  and  100  Hz  the  spectral  amplitude 
was  independent  of  wind  speed.  Shipping  may  be  the  dom¬ 
inant  noise  source  for  this  band  as  the  experiment  was  in 
an  area  of  high  shipping  density.  The  characteristk  spectral 
“hump”  due  to  shipping  noise  was  observed  throughout  the 
experiment.  Narrow  band  60  Hz  harmonics  are  seen  in  the 
spectra  and  are  generated  either  mechanically  or  electrically 
by  FLIP.  There  was  no  evidence  of  60  Hz  lines  during  lab¬ 
oratory  studies  of  the  array  self-noise.  The  60  Hz  harmonics 
are  narrow  band  and  do  not  degrade  the  broadband  signal 
analysis. 

A  graphic  display  of  the  output  of  120  acoustic  channels 
recording  an  air  gun  source  at  a  range  of  500  km  is  shown 
in  Fig.  8.  In  this  plot,  a  compressive  pressure  field  or  pos¬ 
itive  voltage  excursion  is  represented  as  a  filled  line,  and  a 
rarefactive  pressure  field  or  negative  voltage  is  represented 
as  an  unfilled  line.  The  air  gun  is  seen  as  a  series  of  im¬ 


pulsive  arrivals  that  appear  as  both  downward-propagating 
and  upward-propagating  wavefronts  across  the  array.  In  ad¬ 
dition,  two  modes  of  mechanical  vibration  (strum)  are  iden¬ 
tified  in  Fig.  8.  The  first  mode  is  a  longitudinal  vibration 
of  the  kevlar  support  cable,  probably  due  to  vertical  motions 
of  FLIP  pulling  on  the  support  cable.  This  mode  propagates 
with  a  phase  velocity  of  1800  m/s  and  appears  as  a  downward- 
propagating  transient  in  Fig.  8.  The  second  mode  is  a  trans¬ 
verse  vibration  of  the  kevlar  support  cable,  probably  excited 
by  water  currents.  This  mode  propagates  with  a  phase  veloc¬ 
ity  of  40  m/s  and  appears  as  a  series  of  shallow-dipping  lines 
with  an  interference  pattern  every  75  m,  corresponding  to  the 
spacing  of  the  inline  pressure  cases. 

The  air-gun  array  is  seen  as  a  contaminant  in  the  spectra  of 
the  ambient  noise  data  in  Fig.  9(a).  Since  the  operating  area 
of  the  air  run  array  was  shallow,  the  ocean  bottom  altered  the 
arrival,  part  of  which  is  coupled  into  the  deep-sound  channel 
by  down-slope  conversion  and  arrives  at  the  array  in  a  multi- 
path  structure  at  a  variety  of  angles.  The  effect  of  the  profiler 
is  significant  in  the  125  to  250  Hz  band  where  it  dominates 
the  spectra,  clearly  distorting  the  10  dB/octave  roll-off  of  the 
ambient  noise.  The  profiler  was  extracted  (Fig.  9(b))  by  re¬ 
moving  the  visible  signal  from  the  array  time  series,  padding 
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Fig.  7  Ambient  Noise  variation  with  respect  to  wind  and  depth.  Depths 
indicated  are  nominal  for  the  top  of  the  array,  the  hydrophone  shown  is 
675  m  deeper.  Data  were  extracted  on  different  days  to  obtain  similar 
environmental  conditions.  The  narrow  band  signals  are  interference  from 
equipment  aboard  FLIP,  (a)  Wind:  4  kts.  Depth:  400  m,  Jday:  257.  (b) 
Wind:  28  kts.  Depth:  400  m,  Jday:  259  (c)  Wind:  6  kts.  Depth:  1300 
m.  Jday  267  (d)  Wind:  26  kts.  Depth-  1300  m,  Jday:  270.  (e)  Wind: 
4  kts,  Depth:  2300  m.  Jday  267.  (0  Wind  22  kts.  Depth:  2300  M.  Jday 
269  (g)  Representative  spectral  difference  with  wind  (h)  Wind  speed 
(kts)  versus  Julian  Day  (Jday) 
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TABLE  II 

AMBIENT  NOISE  SPECTRAL  VARIATION  WITH  WIND  SPEED 


Change  in 

Change  in 

Frequency 

Reference 

Spectral  Level 

Wind 

dB//uPa/VHz 

kts 

Hz 

+18 

4-6  to  28-33 

200 

[1]  Uncle  1983 

+9 

1040 

200 

[3]  Burgess  and  Kewley  1983 

+7 

4-6  to  22-28 

200 

this  paper  1988 

+7 

10-30 

150 

[11]  Moms  1978 

+5 

18-28 

177 

[21]  Shooter  and  Gentry  1981 

44 

5-28 

177 

[22]  Penone  1969 

0 

6-22 

200 

[15]  Hodgkiss  and  Fisher  1986 

-4 

low-high 

165 

[4]  Wilson  1983 

LONGITUDINAL  AIR 

STRUM  GUN 


TIME,  seconds 


Fig.  8.  Time  series  display  of  acoustic  channels.  A  graphic  display  of 
the  acoustic  channels  shows  an  air-gun  source  at  a  range  of  300  miles  as 
well  as  longitudinal  and  transverse  (shear)  strum  modes.  The  ume  record 
displayed  is  3.S  s.  The  array  hydrophone  depths  span  400  to  1300  m. 


the  remaining  data  to  8192  points,  and  producing  the  spec¬ 
trum  as  described  above  with  no  overlapping.  The  profiler 
spectra  in  Fig.  9(c)  was  calculated  from  800  data  points  and 
a  rectangular  window  and  shows  the  ambient  noise  spectra 
uncontaminated  by  the  profiler. 

IV.  Summary 

This  paper  describes  the  low-frequency  digital  acoustic  ar¬ 
ray  designed  and  built  at  the  Marine  Physical  Laboratory.  A 
sea  test  has  been  conducted  to  verify  the  vertical  deployment, 
telemetry,  acoustic  navigation,  and  ambient  noise  measure¬ 
ment  capabilities  of  the  array.  Navigation  of  the  array  was 
conducted  within  an  acoustic  transponder  net.  Ambient  noise 
spectra  from  single  array  elements  were  consistent  with  pre¬ 
vious  observations  of  shipping  noise  in  the  frequency  band 
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Fig.  9.  Seismic  profiler  contamination.  The  effect  of  the  air  gun  signal  is 
seen  to  be  clearly  contaminating  the  ambient  noise  spectra,  (a)  Original 
data  acquired  at  array  channel  24.  (b)  Same  data  as  (a)  with  profder  signal 
removed,  (c)  Profiler  signature  as  seen  at  the  array. 


between  20  and  100  Hz.  At  frequencies  above  100  Hz  sea- 
surface  noise  generated  by  local  wind  and  waves  was  ob¬ 
served.  Current  efforts  are  underway  to  investigate  the  ver¬ 
tical  directional  spectra  of  low-frequency  noise  measured  by 
the  array. 


Appenddc 

Available  array  commands  are  listed  below.  The  left-most 
hex  number  is  the  command  byte  s  t  by  the  data  acquisition 
computer  to  the  array.  The  last  column  indicates  the  type  of 
command:  Diagnostic,  Initialization,  or  data  Format. 
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APPENDIX 

AVAILABLE  ARRAY  COMMANDS 


41 

DUG  2  ON 

memory  i m  torn  full  uui 

D 

42 

SHOW  HWID 

hwid  pot  ond  eaec 

D 

43 

DUG  1  ON 

analog  dan  test 

D 

44 

DUG  1  OFF 

D 

45 

SCAN  ON 

•can*  through  phenes 

I 

46 

SCAN  OFF 

mUcu  only  am  phene 

I 

47 

SOFT  INTTIAl  ITT 

D 

41 

OFF 

I 

49 

4A 

AD  RATE 

5K 

10* 

I 

1 

4B 

20K 

I 

4C 

NAVRCVR  0 

I 

4D 

NAVRCVR1 

I 

-i£-_ 

ENABLE  RAM  PROGRAM 

p 

JO-SF 

SELHYDN 

N-Othrouih  F 

I 

60 

I 

luge  1 

I 

61 

5 

I 

62 

10 

I 

63 

20 

I 

64 

VAR  GAIN 

5 

tug*  2 

1 

65 

t 

1 

66 

12 

I 

67 

50 

I 

6$ 

1000 

1 

69 

6 

•awl  5’«  (navO)  or  A'i  (n«vl) 

F 

6A 

5 

1 0  phone*  ♦  n»v 

F 

6B 

4 

2  nav  words  ♦  S  phonas 

F 

6C 

DATA  FORMAT 

3 

10  phone* 

F 

60 

2 

tSMtap* 

F 

6E 

1 

Utflt 

F 

iSF 

OFF 

F 

-ZfcTE-i 

NEW  POSITION- 

N 

N^Jthrcuih  F 

leC 

to 

OFF 

g«  data  'ram  marwry  location 

D 

tl 

2000 

(aded  mu  cocc)  addmsa.  oomanu 

D 

12 

2200 

D 

K 

2400 

D 

M 

2S20 

D 

15 

2SD4 

D 

16 

2900 

D 

17 

DUG  3 

2940 

D 

St 

29SO 

D 

19 

29C0 

D 

SA 

— 

D 

SB 

— 

D 

SC 

2FX 

D 

SO 

2FC0 

D 

*e 

— 

D 

IF 

INCREMENT  ADD*  IN  DUO  3  ABOVE 

p 

90 

UPLINK  OFF 

D 

91 

UPLINK  ON 

D 

92 

STORED  SINUSOID  ON 

D 

93 

STORED  SINUSOID  OFF 

D 

99 

FULL  RESET 

D 
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Acoustic  navigation  of  a  large-aperture  array 
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Acoustic  travel  time  measurements  were  used  to  navigate  the  elements  of  a  large-aperture  (900 
m)  acoustic  array.  Array  navigation  system  performance  was  evaluated  during  a  vertical 
deployment  in  the  northeast  Pacific  from  the  Research  Platform  FLIP.  A  network  of  bottom- 
moored  acoustic  transponders  were  interrogated  from  FLIP  and  their  12-kHz  replies  were 
detected  by  receivers  at  75-m  intervals  along  the  array.  A  nonlinear  least-squares  algorithm 
was  used  to  estimate  FLI  P  and  array  element  positions  from  the  travel  time  measurements. 

The  FLIP  positions  derived  from  this  procedure  agreed  with  positions  obtained  from  global 
positioning  system  (GPS)  satellite  navigation  to  within  a  10-m  rms  error.  Navigated  positions 
for  FLIP  were  internally  consistent  wiih  a  0.5-m  mean  rms  error  and  standard  deviation  of  1.1 
m,  and,  for  an  array  element,  were  consistent  with  a  2.8-m  mean  rms  error  and  standard 
deviation  of  0.8  m.  The  resulting  time  series  of  array  and  FLIP  motions  were  analyzed  with 
respect  to  wind,  tidal,  and  internal  wave  forcing  functions.  Wind  and  tidal  forcing  had  the 
greatest  influence  on  FLIP  motion,  whereas  array  motion  was  governed  by  FLIP  movement, 
tides,  and  higher  frequency  sources.  Low-frequency  array  motion,  with  periods  on  the  order  of 
hours,  was  a  result  of  FLIP  towing  the  array  over  a  horizontal  range  of  300  m  in  response  to 
the  wind  and  the  semidiurnal  tidal  oscillations;  the  array  remained  within  a  30-m  horizontal 
range  of  FLIP’S  position.  Higher  frequency  array  motion  had  apparent  internal  wave  and 
surface-coupled  components.  The  array  shape  was  primarily  straight  and  nearly  vertical,  to 
within  approximately  a  2°  tilt,  responding  as  a  simple  pendulum  with  small  displacements. 

PACS  numbers:  43.30.Tg,  43.30.Pc 


INTRODUCTION 

Large-aperture  arrays  produce  high-resolution  direc¬ 
tional  information  by  coherently  combining  signals  from  in¬ 
dividual  array  elements.  Estimates  of  array  shape  and  ele¬ 
ment  position  are  required  because  of  the  dependence  of 
signal  phase  on  element  location.  Navigation  of  such  an  ar¬ 
ray  in  the  ocean  requires  a  known  reference  system  in  a  spa¬ 
tially  and  temporally  varvmg  environment. 

A  900-m,  120-element  low-frequency  acoustic  array  has 
been  developed  by  the  Marine  Physical  Laboratory.1  The 
array  is  capable  of  being  deployed  either  horizontally  or  ver¬ 
tically  up  to  an  ocean  depth  of  6000  m.  Navigation  of  the 
array  is  an  important  part  of  its  operational  requirements, 
and  the  subsystem  implemented  to  perform  the  navigation 
must  meet  design  specifications.  The  position  accuracy  re¬ 
quired  is  on  the  order  of  a  a  few  meters,  the  duration  of 
deployment  is  on  the  order  of  a  few  months,  an  operating 
range  of  up  to  10  km  must  be  viable,  the  size  and  cost  of  the 
system  must  be  manageable,  and  the  sampling  rate  must  be 
sufficient  to  resolve  array  dynamics  due  to  water  currents, 
tides,  wind,  and  internal  waves. 

Techniques  commonly  used  in  ocean  positioning  sys¬ 
tems  are  based  on  mechanical,  electromagnetic,  or  acoustic 
signals.2  Mechanical  positioning  devices  capable  of  the  accu¬ 
racy  required  include  inertial  systems,  pressure  sensors,  tilt 
meters,  measuring  rods,  and  sound-velocity  meters.  The  cost 
and  size  requirements  of  instrumenting  a  900-m  array  with 
such  sensors  is  prohibitive,  although  pressure  sensors  would 
be  valuable  for  depth  positioning  in  combination  with  an 


alternative  method  for  lateral  positioning.  Electromagnetic 
systems  are  either  surface  systems  whose  signals  are  highly 
attenuated  in  water,  but  are  useful  in  determining  the  posi¬ 
tion  of  a  surface  vessel  [e.g.,  global  positioning  system 
(GPS)  ],  or  pulsed  laser  systems,  which  are  currently  under 
development  and  may  eventually  be  effective  within  a  1-km 
range.  The  technique  most  easily  adapted  to  a  large-aperture 
array  deployed  in  the  water  column  is  acoustic  navigation. 
Such  systems  measure  the  travel  time  of  acoustic  energy  be¬ 
tween  two  points  in  the  ocean,  are  capable  of  resolving  a  few 
meters  over  a  10-km  range,  and  are  cost  effective  and  physi¬ 
cally  manageable.  Based  on  these  considerations,  an  acoustic 
navigation  subsystem  was  implemented  as  an  integral  part  of 
the  array. 

This  paper  describes  the  array  acoustic  navigation  sys¬ 
tem  and  its  operation  during  a  vertical  deployment  from  re¬ 
search  platform  FLIP  in  the  northeast  Pacific  during  Sep¬ 
tember  1987.  We  describe  the  hardware  system  and 
acquisition  of  travel  time  measurements  and  associated  er¬ 
rors,  the  nonlinear  least-squares  estin.ation  of  array  spatial 
positions,  and  the  analysis  of  FLIP  and  array  time-varying 
positions  during  the  deployment. 

I.  TRAVEL  TIME  MEASUREMENTS 

Acoustic  travel  time  measurements  provide  range  infor¬ 
mation  between  the  array  and  a  set  of  fixed  transponders. 
Ranges  to  three  or  more  transponders  are  required  to  calcu¬ 
late  FLIP  position  and  array  position  and  shape.  The  naviga¬ 
tion  hardware  consisted  of  three  bottom-moored  acoustic 
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transponders,  two  transmit/receive  timing  units  mounted 
on  FLIP  (one  used  for  FLIP  navigation  and  the  other  used 
for  array  navigation),  an  interrogator/receiver  located  at 
the  bottom  of  FLIP,  and  12-kHz  acoustic  detectors  distrib¬ 
uted  along  the  array.  At  a  known  time,  an  interrogation  sig¬ 
nal  is  transmitted  from  FLIP.  This  signal  triggers  a  tran¬ 
sponder  reply,  which  is  then  detected  by  the  array  or  the 
FLIP  receiver  systems.  The  travel  times  between  the  tran¬ 
sponders,  FLIP,  and  the  array  detectors  determine  array 
position  and  shape. 

During  the  September  1987  experiment,  three  naviga¬ 
tion  transponders  were  moored  in  4700  m  of  water  distribut¬ 
ed  as  a  1.4-nmi.  equilateral  triangle.  FLIP’S  position  was 
maintained  in  the  middle  of  the  transponder  network  by  a 
three-point  moor  (Fig.  1 ).  The  transponders  were  deployed 
by  a  surface  support  ship  operating  a  GPS  satellite  naviga¬ 
tion  receiver  for  absolute  positional  information.  The  tran¬ 
sponders  .were  equipped  with  anchors  ( 175  kg)  to  secure 
them  to  the  ocean  floor  and  to  constrain  the  horizontal 
movement  during  descent  (sink  rate  of  approximately  1  m/ 
s).  Once  anchored,  they  were  buoyed  up  above  the  anchor 
on  100-m  lines  with  130  kg  of  hard  floats  to  limit  horizontal 
movement  and  to  insure  a  well-separated  bottom  bounce. 
They  operate  by  detecting  a  unique  interrogation  frequency 
and  answer  with  a  3-ms  reply  at  12  kHz  with  a  source  level  of 
190  dB  re:  1  /^Pa  at  1  m.  The  recognition  delay  time  is  0.5-3 
ms,  depending  upon  the  signal-to-noise  ratio  of  the  interro¬ 
gation  pulse.3  Due  to  the  proximity  of  FLIP,  the  absence  of 
other  local  high-frequency  ambient  noise  sources,  and  the 


manual  control  of  the  signal  level  on  the  interrogator,  a  con¬ 
stant  delay  time  of  less  than  1  ms  was  assumed. 

The  interrogation  pulse  was  sent  from  a  transceiver 
mounted  on  the  bottom  of  FLIP  ( 90  m  in  depth ) .  The  pulse 
level  was  adjusted  manually  above  the  ambient  noise  for  con¬ 
sistent  transponder  replies.  The  pulse  was  triggered  by  either 
a  chart  recorder  to  navigate  FLIP  or  by  an  interrogation 
circuit  to  navigate  the  array..  The  chart  recorder  trace,  set  on 
a  1-s  sweep  rate,  displayed  the  12-kHz  transponder  reply 
amplitudes  in  a  500-Hz  band  with  a  transmit/receive  delay 
time  of  0.1  ms.  Each  transponder  was  interrogated  individu¬ 
ally  to  facilitate  identification  and  transmitted  a  reply  pulse 
once  per  second  for  45  s.  Since  the  FLIP  transducer  was  also 
used  as  an  interrogator  for  the  array  navigation  system,  only 
one  system  was  operable  at  a  time;  consequently,  the  FLIP 
navigation  data  were  recorded  only  once  an  hour  during  the 
sea  test.  Round-trip  travel  times  for  navigating  FLIP  were 
measured  by  hand  from  the  chart  recorder  output  with  an 
estimated  rms  error  of  2-3  ms. 

The  900-m  array  was  deployed  vertically  from  FLIP 
with  a  420-kg  weight  attached  to  maintain  a  nearly  vertical 
orientation.  The  array  navigation  subsystem  was  composed 
of  12  hydrophones  with  75-m  spacing  monitored  by  12  indi¬ 
vidually  programmable  processors.  Each  of  the  12  proces¬ 
sors  contained  tuned  navigation  receivers  designed  to  detect 
the  12-kHz  transponder  reply  in  a  200-Hz  band  with  a  —  6- 
dB  minimum  signal-to-noise  ratio.  Once  per  minute  the 
transceiver  on  FLIP  sent  out  four  continuous  wave  (cw) 
pulses:  three  transponder  interrogation  pulses  at  unique  fre- 
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FIG.  1 .  Navigation  overview.  The  navigation  sys¬ 
tem  consists  of  three  bottom-moored  transpon¬ 
ders  in  an  equilateral  triangle  about  FLIP  and  the 
vertical  array.  The  transponders  are  designated 
by  color:  red,  green,  and  blue.  The  initial  xyz  co¬ 
ordinates  of  the  transponders  and  of  the  FLIP  in¬ 
terrogator  are  indicated. 
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quencies,  followed  by  a  12-kHz  pulse  for  depth  ranging.  The 
12-kHz  transponder  replies  (3  ms)  and  the  12-kHz  depth 
ranging  pulse  from  FLIP  ( 10  ms)  were  recorded  in  series  as 
a  1-bit  output,  every  0.4  ms,  by  each  of  the  array  navigation 
receivers.4  The  interrogation  sequence  clock  was  synchro¬ 
nized  with  the  array  time  base  and  recorded  every  12s  ms. 

To  calculate  the  travel  time,  the  time  of  the  transponder 
return  must  be  determined  relative  to  the  time  of  the  interro¬ 
gation.  The  signal  that  initiated  the  interrogation  sequence 
was  found  by  linear  interpolation  of  the  interrogation  se¬ 
quence  clock.  The  leading  edge  of  the  return  pulse  was  found 
by  recreating  each  of  the  navigation  receiver  time  series,  and 
correlating  these  binary  time  series  with  a  pulse  replica  that 
depended  upon  array  receiver  noise  and  level  of  detection  A 
1-s  time  interval  of  the  normalized  correlation  output  from 
each  array  section  is  shown  for  each  transponder  in  Fig.  2. 
The  correlated  arrival  across  the  array  of  each  transponder 
reply  is  evident,  interference  and  noise  are  also  evident  To 
reduce  the  chances  of  detecting  a  noise  spike  or  multipath 
return  rather  than  the  direct  transponder  reply,  a  window 
was  invoked,  for  each  expected  return,  during  which  time 
the  data  were  considered  valid,  as  shown  by  the  small  X ’s  in 
Fig.  2.  The  window  was  allowed  to  track  the  replies  as  the 
arrival  time  shifted  with  array  movement.  The  travel  time 
from  the  transponder  to  the  array  is  the  difference  between 
the  departure  time  of  the  interrogation  pulse  and  the  arrival 
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FIG  2  Travel  time  acquisition  The  12  array  navigation  receivers  each  de¬ 
tect  four  pulses  during  a  65  536-s  navigation  cycle  The  normalized  corre¬ 
lated  receiver  outputs  are  shown  during  a  1-s  interval  bracketing  each  detec¬ 
tion  The  top  trace  represents  the  output  of  the  top  receiver,  followed  in 
sequence  by  the  deeper  receivers  The  first  three  panels,  (a)-(c).  represent 
the  detected  direct  replies  of  the  three  transponders  (shown  as  R,„Gd  and 
Bd ),  arriving  at  the  deepest  receiver  first  and  traveling  up  through  the  w  ater 
column  to  the  shallower  receivers  The  last  panel  (d)  shows  the  detection  of 
the  direct  arrival  Fd  of  the  12-kHz  pulse  transmitted  from  FLIP  traveling 
down  the  array,  with  a  surface  reflection  F„  arriving  about  1 20  ms  after  the 
direct  pulse.  The  small  X ’s  represent  the  valid  data  window  The  start  time 
is  indicated  at  the  top  of  each  panel 


time  of  the  transponder  reply  at  the  leading  edge  of  the  corre¬ 
lated  signal,  minus  the  travel  time  from  FLIP  to  the  tran¬ 
sponder  and  the  constant  system  delay  times. 

The  accuracy  of  the  travel  time  measurements  from  the 
bottom  transponders  to  the  array  was  dependent  upon  fac¬ 
tors  which  affected  all  receivers  similarly  as  well  as  factors 
which  varied  with  each  array  receiver.  The  factors  with  simi¬ 
lar  effect  were:  (a)  the  ambient  ocean  noise  in  the  12-kHz 
band,  which  has  a  typical  level  of  30-40  dB///r  Pa  A  Hz ,  and 
(b)  the  transponder  performance,  which  introduced  timing 
errors.4  The  factors  that  varied  with  receiver  were  (a)  the 
drift  in  the  array  clock  phase  lock  loop  .circuitry,1  which 
introduces  small  sampling  time  errors,  (b)  the  FLIP-to- 
transponder  travel  time  measurement  error  (the  resultant 
error  in  array-to-transponder  range  increases  as  depth  de¬ 
creases)  of  order  1.5  ms,  and  (c)  the  noise  and  detection 
threshold  of  the  array  receivers.  The  variation  in  receiver 
error4  is  confined  for  the  most  part  to  a  mismatch  in  detec¬ 
tion  threshold  due  to  the  temperature  sensitivity  of  the  capa¬ 
citors  and  high  failure  rate  of  the  inductors  in  the  tuned 
navigation  filters.  This  variation  was  estimated  at  2  ms  dur¬ 
ing  laboratory  tests  of  the  frequency  sensitivity  of  the  tuned 
filters  and  by  examining  the  length  of  the  transponder  re¬ 
plies.  Estimated  accuracy  for  the  array  travel  time  measure¬ 
ments  depends  upon  the  receiver  and  the  depth  of  the  array, 
but  is  typically  2-3  ms.  The  accuracy  of  the  depth  measure¬ 
ments  is  affected  by  the  ambient  ocean  noise,  the  transmitter 
and  detector  noise,  the  detection  threshold,  and  the  strength 
of  the  transmission.  Measurement  error  in  the  travel  times 
corresponding  to  array  depth  is  estimated,  by  calculating  the 
mean  standard  deviation  in  ten  sample  segments,  to  be  <1 
ms. 


II.  POSITION  ESTIMATION 

In  this  section  we  determine  where  the  three  bottom 
transponders,  the  FLIP  transceiver,  and  the  array  elements 
are  in  three-dimensional  space  using  the  travel  times  de¬ 
scribed  in  the  previous  section.  The  general  problem  of  fit¬ 
ting  model  functions  and  parameters  (spatial  positions)  to 
data  (travel  times)  appears  in  many  areas  of  applied  science: 
as  parameter  estimation,  regression  analysis,  inverse  prob¬ 
lems,  filtering,  process  identification,  and  as  an  optimization 
problem  in  numerical  analysis.  Noise  m  the  travel  time  mea¬ 
surements  demands  that  sufficient  data  be  collected  to  over¬ 
determine  the  solution,  requiring  some  kind  of  approxima¬ 
tion  method. 

To  approximate  the  solution,  an  optimization  prob’em 
is  defined  to  minimize  differences  between  model  predic¬ 
tions  and  actual  measurements.  The  function  relating  pre¬ 
dictions  and  measurements  is  determined  by  the  problem 
structure  and  the  error  distribution.  One  approach  is  to 
minimize  the  diffeiences  with  respect  to  a  matrix  norm. 
Analytically,  the  most  viable  choice  is  the  U  norm  resulting 
in  the  familiar  least-squares  solution.  If  tu-  errors  are  inde¬ 
pendent  and  distributed  normally  with  a  constant  standard 
deviation,  then  the  l2  norm  coincides  with  the  maximum- 
likelihood  solution.  The  least-squares  solution  magnifies  the 
effect  of  large  errors,  and  if  the  error  distribution  contains 
many  outliers,  the  /,  norm  solution  as  described  by  Gill  et 
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a/.5  and  imp!eis“ited  by  Duckworth6  is  an  attractive  alter¬ 
native.  Although  the  /,  solution  is  more  robust  m.  terms  of 
outliers,  these  errors  continue  to  affect  the  solution,  and  if 
possible,  the  source  of  the  large  errors  should  be  identified 
and  removed. 

'  lie  navigation  problem  discussed  here  has  been  imple¬ 
mented  in  terms  of  the  least-squares  method.  The  transpon¬ 
der  net  was  localized  separately  to  allow  accurate  real-time 
array  navigation.  The  FLIP  travel  time  measurements  used 
for  the  transponder  localization  are  very  consistent  and  pos¬ 
sess  no  large  errors,  having  been  selected  and  verified  by 
hand.  A  priori  knowledge  of  array  measurement  errors  did 
not  justify  deviation  from  the  assumed  normal  distribution. 
In  the  remainder  of  this  section,  a  general  descriptiomof  the 
least-squares  formulation  is  followed  by  the  specific  details 
of  implementation  for  the  array  navigation  problem.  The 
effect  of  specific  types  of  errors  are  simulated,  and  results  are 
presented. 

A,  Least-squares  formulation 

The  FLIP,  array,  and  transponder  positions  were  esti¬ 
mated  by  a  nonlinear  least-squares  method.  The  problem  is 
to  find  the  set  of  model  parameters  ( spatial  positions ) ,  given 
a  set  of  measurements  (travel  times)  plus  noise.  The  least- 
squares  estimator  identifies  a  realistic  set  of  parameters  by 
minimizing  the  sum  of  the  squared  errors,  with  an  error  de¬ 
fined  by  the  difference  between  the  model  estimate  and  the 
measurement.  For  the  set  of  measured  data  vk  (travel 
times),  and  the  set  of  parameters  a:  (spatial  positions)  that 
determine  the  model  estimate  vk  {a/J  —  l, AO  ( referred  to  as 
v  or  vk ),  the  least-squares  formulation  is  written  as 

mini takefm.)  -  T  lUtZ*!? 

«=*"  V  *ri  <r 

=  (v-v)rR-'(v  — v)),  (I) 

where  the  measurement  errors  are  uncorrelated  with  identi¬ 
cal  variances  cF  and  R  is  (Fl.  error  function  F(v;a) 
(designated  as  F),  is  minimized  when  the  gradient  g  of  F 
with  respect  to  the  model  parameters  is  zero: 

g-~=  -2(>— v)rR-1-g-  =  °.  (2) 

da  da 

In  general,  a  closed-fonned  (one-step )  least-squares  so¬ 
lution  exists,  provided  the  set  of  equations  is  linear.  For  the 
array  navigation  problem,  however,  the  nonlinear  mapping 
from  measurement  space  to  parameter  space  necessitates  an 
iterative  implementation.  One  such  implementation  is  to  ap¬ 
proximate  the  smooth  nonlinear  functions  described  by  F 
with  a  linear  form  of  the  Taylor’s  series  expansion  about  an 
initial  position  a,.  In  the  immediate  neighborhood  of  a7, 

F(v;  a  +  Ap)  =  F(v;a,)  +  Y  b.p  ,  (3) 

j~ i  daJ  ,, 

where  A  is  a  scalar,  and  p  is  a  direction  vector  guiding  the 
search  for  a  minimum.  In  terms  of  the  iteration  step  the 
notation  may  be  written 

Fn  +  i  =F„ +hrglp„,  (4) 


where  g„  is  the  gradient  ofF  at  v„ .  An  acceptable  minimum 
for  F  is  obtained  by  calculating  theestnnate  v  from  an  initial 
“best  guess”  for  the  parameters  a7,  then  developing  a  pre¬ 
scription  which  improves  this  initial  guess,  updating  the  pa¬ 
rameters  until  specified  convergence  criteria  are  satisfied.  A 
typical  prescription  to  search  for  this  minimum  is  to  calcu¬ 
late  v  from  a, ,  F„  from  v  and  v„ ,  and  complete  the  following 
steps: 

( 1 )  Compute  a  vector  p„ ,  which  is  the  search  direction. 

(2)  Compute  a  scalar  step  length h„ . 

(3)  Update  the  estimate  of  the  parameters.  Set 
aB^,—  a„  +  hnpn,  followed  by  n—n  +  1,  such  that  the  cur¬ 
rent  error  function  is  now  Fn . 

(4)  Test  for  convergence.  Calculate  the  forward  prob¬ 
lem,  which  determines  current  estimates  v„  of  the  measure¬ 
ments  ▼  from  current  estimates  of  the  parameters  a„ .  Evalu¬ 
ate  F„  and  Fn  _ ,  with  respect  to  the  convergence  criteria.  If 
the  conditions  for  convergence  are  satisfied,  the  current  val- 

*  ue  of  a„  is  considered  the  solution;  if  not,  iterate  by  returning 
to  step  (1). 

To  implement  this  prescription,  the  choice  of  conver- 
gence  criteria,  search  direction,  and  step  length  must  becon- 
sidered;  particular  choices  relating  to  array  navigation  will 
be  specified  m  the  following  section.  In  an  optimization 
problem,  the  value  at  theminimumus  generally  not  known, 
and  convergence  criteria  must  be  developed  to  bound  the 
number  of  iterations  while  insuring  that  the  resulting  solu¬ 
tion  is  sufficiently  close  to  thereal  minimum.  Specific  crite¬ 
ria  are  defined  by  evaluating  the-magnitude  of  the  gradient, 
the  magnitude  of  the  squared  error,  and  therelativedecrease 
in  squared  error  with  respect  to  acceptable  tolerances  <5,  as 
follows: 

lift.  I!  <5i. 

Fn  «S2.  (5) 

~F„  <<53, 

number  of  iterations  >  l/<54. 

Combinations  of  the  above  criteria  determine  convergence. 
The  rate  of  convergence  depends  on  the  specific  minimiza¬ 
tion  algorithm,  the  initial  estimates,  and  the  error  surface. 
The  resulting  minimum  point  could  be  either  global  or  local. 
In  general,  there  areno  sure  methods  for  finding  globai  mini- 
ma;  so  the  accuracy  of  the  initial  guess  is  important. 

Considerthe  iterativestep  in  the  direction  of  search  p„ . 
In  the  univariate  case,  the  only  possible  directions  are  posi¬ 
tive  or  negative.  As  the  number  of  dimensions  increases, 
even  to  two,  the  number  of  possible  directions  is  theoretical¬ 
ly  infinite.  Search  direction  procedures  are  typically  direct 
search  methods  or  gradient  search  methods  with  the  con¬ 
straint  that  p„  bea  descent  direction.  Iterative  search  proce¬ 
dures  having  the  general  form  of  Eq.  (4)  lend  themselves 
easily  to  gradient  search  methods.  To  define  the  gradient 
search,  ifg„  is  nonzero,  there  must  exist  a  vector  p„  such  that 
gnPn  <  0.  If  p„  is  defined  as  —  g„  and  the  gradient  is  non¬ 
zero,  this  direction  is  clearly  a  descent  direction;  the  algo¬ 
rithm  is  recognized  as  the  steepest  descent  method.  It  fol¬ 
lows  from  Eq.  (4)  that,  for  any  small  positive  step  size  h, 
Fn+l  <F„.  Since  p„  is  in-the  direction  of  decreasing  F,  con- 
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vergencc  is  guaranteed  provided  hn  is  chosen  so  that  F  is 
“sufficiently  decreased”  at  each  iteration.  To  be  appropriate, 
h„  must  be  small  enough  to  locate  the  minimum  precisely 
and  satisfy  the  immediate  neighborhood  constraint  of  Eq. 
(4),  and  at  the  same  time  large  enough  to  search  efficiently. 
The  existence  of  guaranteed  convergence  of  the  algorithm 
does  not  imply  that  convergence  will  be  achieved  in  ar  ac¬ 
ceptable  number  of  iterations  A  common  method  used  to 
define  h  is  to  follow  the  steepest  descent  gradient  for  the 
current  iteration  to  a  minimum.  This  dictates  that  the  m-xi 
search  dnection  be  nearly  orthogonal  to  the  current  direc¬ 
tion  gjjp„  ^ ,  —0.  It  has  been  shown  that  the  directions  gener¬ 
ated  by  this  method  asymptotically  converge  to  only  two 
directions  for  many  problems,7  increasing  the  number  of 
iterations  immensely.  This  situation  may  be  avoided  by  con¬ 
straining  the  step  size  to  be  less  than  that  required  to  reach 
the  minimum  m  the  current  search  direction.  Although  the 
rate  of  convergence  of  the  steepest  descent  algorithm  is 
known  to  be  less  than  other  methods,5  the  method  is  numeri¬ 
cally  stable  and  produces  efficient  results  for  the  nav,bciion 
problem  considered  here. 

B.  Implementation 

We  disc"ss  implementation  ct  the  general  nonlinear 
least-squares  method  described  above  for  the  transponder 
and  array  navigation.  The  ingredients  required  for  a  least- 
squares  problem  are-  ( 1 )  the  measurements  v  <  travel  time 
detected  at  the  array  and  at  FLIP),  (2)  the  initial  estimates 
of  the  model  parameters  a  represented  by  FLIP,  array,  and 
transponder  xyz  positions;  (31  the  mapping  function  from 
measurement  space  to  parameter  space  which  calculates  the 
estimates  v  from  the  parameters  a;  and  (4)  the  measurement 
errors  a.  Although  the  general  prescription  is  the  same  for 
the  transponder  localization  and  for  the  FLIP/array  local¬ 
ization,  there  are  differences  in  the  implementation.  In  the 
followmg,  the  explicit  inputs  and  implementation  of  the 
transponder  localization  are  specified  and  the  differences  re¬ 
lating  to  FLIP  and  the  array  are  discussed  separately. 

1.  Transponder  localization 

The  two-way  travel  times  between  FLIP  and  the  tran¬ 
sponders  were  used  toi  transponder  xy  position  estimation 
Normally,  data  from  an  extensive  surface  ship  survey  would 
be  available  for  localizing  the  transponders.8  g  Although 
such  a  survey  was  conducted  during  the  experiment,  the  data 
were  not  complete  due  to  extreme  noise  levels  of  the  surface 
ship;  instead  data  recorded  by  the  FLIP  navigation  system 
wcre  used.  The  FLIP  data  set  defines  a  survey  configuration 
that  is  not  optimal;  however,  the  extended  time  series  and 
low  errors  of  the  FLIP  measurements  increase  their  reliabil¬ 
ity  over  the  other  available  data  sets  The  unconventional 
geometry  of  these  measurements  imparts  a  n_ rational  sym¬ 
metry  to  the  problem,  as  FLIP  was  moored  in  the  center  of  a 
roughly  equilateral  triangle  defined  by  the  transponder  posi¬ 
tions.  This  physical  constraint  limits  the  amount  of  indepen¬ 
dent  xy  information  which  is  contained  in  the  measurements 
and  tends  to  destabil  ze  the  problem.  Our  resu’ts  were 
-chieved  by  choosing  well-constrained  initial  climates  of 
the  xyz  positions  described  below,  by  fixing  the  -  component 
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during  the  xy  iteration,  and  by  perturbing  the  FLIP  posi¬ 
tions  independently  of  the  transponder  positions.  Although 
the  discussion  of  the  general  least-squares  method  was  de¬ 
scribed  in  terms  of  measured  and  estimated  travel  times,  the 
actual  implementation  described  here  was  in  terms  of  hori¬ 
zontal  ranges.  Describing  the  implementation  in  terms  of  the 
variables  intr  iduced  in  Eqs.  ( 1 )  and  (2),  the  model  param¬ 
eters  aJtj  =  l  ,N  are  the  x  and  y  positions  of  the  FLIP  trans¬ 
ceiver  and  the  three  bottom  transponders,  these  x  and  y  posi¬ 
tions  define  the  horizontal  projection  cf  the  slant  range  from 
a  specified  transponder  to  FLIP  which  is  designated  i\, 
k—\,M .  The  measured  travel  times  are  converted  to  hori¬ 
zontal  projections  designated  vk ,  k  =  \,M,  using  the  sound- 
speed  profile  and  the  fixed  z  components  as  shown  below 
The  number  of  horizontal  projections,  M,  is  3  for  the  FLIP 
perturbations  and  424  for  the  transponder  perturbations 
The  conversion  from  the  time  domain  to  the  space  do¬ 
main  is  dependent  on  the  sound-speed  profile  between  the 
specified  depths.  Assuming  a  constant  sound  speed  and 
therefore  a  straight-line  path,  and  expanding  the  vk's  for 
illustrative  purposes,  the  measured  travel  time  is  related  to 
the  xyz  position  as 

measured  travel  time  =  t  (slant  range), 

slant  range  =  [(x,  -  xr)2  +  u>,  -y.):  +  (z,  -zr)2]'/:, 

(6) 

horizontal  projections  =  [  (slant  range)2  —  (&)2] 1/2 

=  [(*«  -*,)*+ O', 

^  =  h , 

where  c  is  the  sound  speed,  x,,  y„  z,  iepresent  the  transmit¬ 
ter  position,  and  xr,yr,zr  represent  the  receiver  position. 
Sound  speed  is  not  constant  m  the  ocean  environment,  and  a 
varying  sound-speed  profile  refracts  acoustic  eneigy;  how¬ 
ever,  under  certain  conditions  the  harmonic  mean,  defined 
below,  may  be  used  with  minimal  error  To  estimate  the  trav  - 
el  time  differences  due  to  refraction  effects,  the  configuration 
parameters  of  the  September  experiment  were  specified  and 
employed  as  inputs  to  the  generic  sonar  model  (GSM).10 
The  GSM  used  generalized  ray  theory  to  calculate  the  ray 
path  and  travel  time  between  a  transmitter-receiver  pair  by 
assuming  a  sound-speed  profile  with  layers  of  constant  gra¬ 
dient.  The  simulation  results  are  presented  in  Table  I,  which 
shows  the  travel  time  between  an  estimated  transponder  po¬ 
sition  and  the  array  for  a  12-kHz  signal  as  a  function  of 
range,  depth,  and  sound-speed  profilt  The  first  column 
( right  of  the  double  vertical  bar)  represents  the  results  gen- 
e  ated  by  assuming  a  constant  sound  sneed  of  1500  m/s 
throughout  the  water  column.  T  ie  second  column  was  cal¬ 
culated  using  F  q.  ( 8 )  shown  below,  with  historical  values  of 
temperature  and  salinity  versus  depth  for  the  area  covering 
the  experiment  site  and  at  the  same  time  of  year. 1 1  The  third 
column  used  the  harmonic  mean,  calculated  by  numerical 
integration  of  the  layered  sound-speed  profile  over  depth  to 
estimate  travel  time  which  is  then  converted  to  speed 


The  C(z)  profile  used  in  the  equation  above  is  shown  in  Fig. 
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TABLE  I.  Sound-speed  impact  or,  travel  time  predictions 


Travel  time  versus  sound  speed,  depth,  and  range:  4573-m  source  depth,  12  kHz 
GSW  travel  time  predictions  using 


Horizontal 

range 

(m) 

Receiver 

depth 

(m) 

Constant 

profile  ( 1500  m/s) 

(s) 

Histone 

profile 

(s) 

Harmonic 

mean 

(s) 

Deep  ctd 
(s) 

2000 

90 

3  2724 

3.2740 

3  2733 

3.2733 

2600 

3.4546 

3  4563 

3  4556 

3.4555 

2000 

400 

3.0848 

3.0839 

3  0834 

3  0833 

2600 

3  2774 

3  2765 

3  2759 

3  2758 

3000 

3.4243 

3.4241 

5000 

4  3387 

4.3382 

10000 

7.2181 

7.2152 

2600 

1600 

2.6325 

2.6184 

2.6175 

2.6174 

2600 

2800 

1.0972 

1.0726 

1.0717 

10716 

3  and  was  calculated  using  the  Unesco  equations12  relating 
conductivity  to  practical  salinity,  sound  speed  to  salinity, 
temperature  and  pressure,  and  pressure  to  depth.  The  con¬ 
ductivity,  temperature,  and  pressure  measurements  were  ob¬ 
tained  from  a  surface  ship  on  Julian  day  267  at  2243  GMT 
within  1  nmi.  of  the  FLIP  location.  Sound-speed  profiles  in 
the  ocean  are  typically  very  stable  below  the  thermocline 


VELOCITY,  meiers/second 


FIG.  3.  Array  navigation  receivers  and  sound-speed  profiles.  The  array  was 
deployed  at  three  nominal  depths  during  the  September  1987  experiment 
The  sound-speed  profile  is  used  in  the  conversion  from  travel  times  to  spa¬ 
tial  position.  This  profile  was  calculated  from  CTD  data  (0-4900-m  depth) 
and  is  shown  with  three  profiles  calculated  from  XBT  data  (0-700-m 
depth)  chosen  to  illustrate  the  variation  in  sound  speed  in  the  upper  ocean 
layer.  The  navigation  data  from  one  receiver  per  array  section  (H5)  ware 
recorded  during  the  experiment. 
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over  the  time  period  of  the  sea  test.  The  diversity  in  tempera¬ 
ture  at  the  surface,  however,  imparts  temporal  variations  in 
the  sound-speed  profile.  Temperature  measurements  were 
acquired  from  expendable  bathythermograph  (XBT)  casts 
deployed  from  FLIP  twice  a  day  throughout  the  experiment. 
Temperature  data  were  converted  to  sound  speed  and  deci¬ 
mated  to  maintain  compatibility  with  program  inputs.  This 
conversion  utilized  historical  salinity  data1  ’  and  was  imple¬ 
mented  by  means  of  the  following  equation13: 

c  =  1449.2  -f  4.67-  0.05572  +  0.000  29  73 

+  (1.34-0.107)  (S-  35)  +0.016Z.  (8) 

where  c  —  sound  speed  (m/s),  7=  temperature  (°C),  S 
=  salinity (ppt),  and  z  —  depth(m).  The  maximum  differ¬ 
ence  in  the  harmonic  mean  due  to  time  variation  in  surface 
temperatures  was  about  1 . 5  m/s  ( G.  1  %  of  the  mean  speed ) . 
This  affects  the  predicted  travel  time  from  a  transponder  to  a 
receiver  at  90  m  by  about  3  ms  ( 5  m) .  An  error  in  the  sound- 
speed  correction  of  the  measurements  translates  into  a  larger 
array  element  rms  error,  but  does  not  substantially  affect  the 
estimated  position  due  to  the  nearly  equilateral  triangle  sys¬ 
tem  configuration.  Due  to  possible  inaccuracies  in  the  abso¬ 
lute  XBT  measurements,  and  the  robust  nature  of  the  system 
configuration  to  enors  in  sound  speed,  the  profile  was  con¬ 
sidered  constant  over  time  and  equal  to  the  CTD  calcula¬ 
tions.  The  fourth  column  in  Table  I  incorporates  refraction 
effects  due  to  changes  in  the  sound-speed  profile  C(z)  de¬ 
scribed  above.  The  predictions  generated  by  the  GSM  indi¬ 
cate  that,  for  the  geometrical  configuration  of  the  September 
experiment  shown  in  Fig.  4,  the  error  in  travel  time  by  as¬ 
suming  a  straight-line  path  with  a  constant  average  sound 
speed  is  0.1  ms;  therefore,  using  die  harmonic  mean  of  the 
sound-speed  profile  introduces  a  negligible  error  (0. 1 5  m )  in 
the  array  positions. 

GPS  satellite  positions  at  FLIP  and  the  surface  ship  dur¬ 
ing  transponder  deployment  provided  the  initial  estimates  of 
the  transponder  xy  positions  to  within  tens  of  meters.  The 
clement  weather  that  prevailed  during  the  deployment  and 
the  negative  buoyancy  of  the  transponders  (sink  rate  of  cs  1 
m/s)  constrained  the  final  transponder  positions.  Initial 
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transponder  depths  were  estimated  by  subtracting  the  tran¬ 
sponder  anchor  line  length  from  the  echo-sounding  depth 
measured  at  FLIP.  The  assumption  that  the  seafloor  is  flat  in 
the  immediate  area  of  FLIP  was  corroborated  by  bottom- 
mounted  Swallow  floats  deployed  approximately  2  nmi  from 
FLIP.  Preliminary  Swallow  float  depth  estimates  were  with¬ 
in  15  ni  of  the  depth  measured  at  FLIP.14 

Once  the  ingredients  for  the  least-squares  method  have 
been  accumulated,  the  xy  positions  are  adjusted  until  the 
root-mean-squared  (rms)  error  satisfies  the  convergence 
criteria  described  below.  This  is  accomplished  in  several 
steps  by  first  maintaining  constant  transponder  positions 
and  perturbing  the  FLIP  positions  to  minimize  errors,  then 
holding  the  current  FLIP  positions  constant  while  perturb- 
in;  the  transponder  positions  to  minimize  errors,  then  exam- 
in:  lg  the  mean-sq  jared  error  of  each  FLIP  position  to  deter- 
mi  “  (via  a  user  parameter)  whether  it  should  be  preserved 
as  .able  contributor,  and,  finally ,  repeating  the  entire  pro- 
cedu  e  until  the  convergence  criteria  are  satisfied.  The  per- 
turbai'on  adjustment  is  calculated  with  the  search  direction 
p  defined  as  the  negative  gradient  The  step  size  is  a  constant 
(1.5  mj  onless  the  rms  error  is  less  than  1  m,  at  which  time 
the  step  size  begins  to  decrease,  as  the  minimum  is  ap¬ 
proached,  the  stcy  ■  _  ,  ^aicuiaied  as  a  function  of  the  per¬ 
cent  change  m  berated  rms  error  The  convergence  criteria 
are  defined  as  follow  s.  an  absolute  rms  error  less  than  0  15m. 
a  0.15%  change  in  the  iterated  rms  error,  or  a  maximum  of 
30  iterations  A  typical  number  of  loop  iterations  required 
for  the  convergence  of  the  September  data  set  were  3  for  each 
inner  loop  and  5  for  the  entire  procedure,  the  repetition  fac¬ 
tor  depending  on  the  user  parameters  A  detailed  description 
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FiU  4  Spatial  configuration  for  navigation  during  ihc  September  1987  sea 
test  The  estimated  positions  of  FLIP  and  tne  three  botloir  transponders 
(red,  green,  and  blje.  measured  houriv  oser  18  da\s  arc  snovsn  1.1  plan 
View  The  FLIP  mooring  lines  arc  represented  b>  arrows,  and  the  slant 
range  (j)  and  horizontal  projection  of  the  slant  range  ( /i }  arc  indicated  in 
meters  from  an  arbitrar.  FLIP  x\  position  1 2400.  3200)  Transponder  ba¬ 
seline  distances  arc  indicated  in  meters 
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of  the  software  implementation  and  its  usage  is  docume  ,t- 
ed.4 

The  transponder  depths  were  iterated  manually  by  min¬ 
imizing  the  magnitude  and  dispeision  of  rms  error  across  the 
array.  Increasing  rms  errors  across  the  array  was  shown  by 
the  simulations  in  the  follow  ing  section  to  be  an  indication  of 
potential  transponder  depth  error.  The  FLIP  database  used 
to  navigate  the  transponders  was  sampled  during  tidal  peaks 
to  discern  whether  there  was  any  transponder  movement 
due  to  tidal  forces.  No  appreciable  difference  was  observed, 
and  transponder  movement  is  assumed  to  be  negligible 

2.  FLIP  and  array  element  localization 

The  localization  procedure  for  FLIP  and  the  array  is 
essentially  the  same  as  just  described  for  the  transponders. 
With  the  transponder  positions  known,  the  FLIP  and  array 
element  navigation  may  be  calculated  m  real  time 

The  initial  xy  positions  for  FLIP  and  the  array  elements 
were  estimated  from  GPS  positions  af.er  FLIP  was  moored 
(Julian  day  255/.  Because  FLIP'S  position  was  maintained 
by  a  three-pomi  moor,  the  same  initial  xy  position  was  used 
for  each  FLIP  position  throughout  the  test.  If  the  initial  xy 
position  must  be  estimated  vv  ithout  GPS.  or  if  the  receiv  er  is 
not  moored,  a  geometric  estimate  may  be  obtained  from  the 
transponder  slant  ranges.  The  depths  of  the  array  elements 
were  measured  to  within  1.5  m  by  the  12-kHz  transmitter 
mounted  on  the  bottom  of  FLIP.  The  12-kHz  signal  was 
received  by  each  navigation  receiver  without  exception, 
however,  there  were  indications  of  slight  interference  in  the 
v  ariabiluy  of  the  arrival  time  of  the  resulting  correlated  out¬ 
put.  The  array  was  deploy  ed  under  420  kg  of  tension  to 
maintain  verticality  and  constrain  the  depth  parameter.  A 
1.5-m  error  in  array  depth,  because  this  input  ts  a  fixed  pa¬ 
rameter  in  an  equilateral  configuration,  translates  into  less 
than  0.5  m  of  error  in  the  xy  positions  of  the  array,  but 
changes  the  rms  error  by  as  much  as  2.5  m. 

The  noise  in  the  travel  times  measured  at  the  array  is 
greater  than  the  noise  in  the  hand-picked  data  measured  at 
FLIP.  To  improve  the  data  quality,  various  averaging, 
thresholding,  and  interpolation  schemes  were  incorporated 
in  the  processing  software.  Should  a  receiver  not  detect  a 
return,  or  if  the  detected  return  is  not  within  specified 
thresholds  for  range  and  depth  variability,  various  interpola¬ 
tion/ext  rap  olation  software  options  may  be  specified.  Prior 
to  the  least-squares  iterations,  the  data  may  be  smoothed 
with  a  running  average  filter.  The  travel  time  measured  at 
the  array  represents  the  time  from  FLIP  to  the  transponder 
to  the  array.  Therefore,  to  acquire  .ne  transponder-to-array 
travel  time  required  by  the  least-squares  procedure,  the 
FLIP-to-transponder  travel  time  (measured  at  FLIP  once 
ar.  hour)  is  interpolated  and  subtracted  from  the  travel  time 
measured  at  the  array  (approximately  once  per  minute). 

C.  Simulations 

Simulations  were  conducted  to  examine  the  array  posi¬ 
tion  sensitivity  to  errors  m  the  transponder  and  FLIP  posi¬ 
tions.  The  spatial  configuration  showm  m  Fig.  4  simulates  the 
September  sea  test;  FLIP  is  centrally  located  in  the  tran- 
spondei  network  (identified  by  colors),  and  the  array  is  di- 
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rectly  under  FLIP.  The  transponders,  FLIP,  and  the  array 
are  it  dually  assigned  to  known  positions  with  slant  ranges  as 
shown.  Individual  transponder  positions  were  perturbed, 
and  the  resulting  array  shape  and  relative  errors  were  exam¬ 
ined  These  tests  were  conducted  first  wi:h  the  array  in  a 
vertical  orientation;  that  is,  the  array  element  xy  positions 
were  identical  to  those  of  FLIP,  and  then  with  the  array 
straight  but  tilted  2°  vertically  at  225°  azimuth.  Zero  mean 
Gaussian  random  errors  (a2  —  3  m)  were  added  to  the  array 
travel  times  with  no  appreciable  changes.  The  final  simula¬ 
tion  configuration  perturbed  the  FLIP  position.  Because  the 
array  measures  travel  times  of  rays  which  originate  at  FLIP, 
errors  in  the  FLIP  parameters  map  into  array  positional  er¬ 
rors. 

The  simulations  provide  an  understanding  not  only  of 
the  error  magnitudes  but  also  of  their  dispersion  across  the 
array.  With  a  +  25-m  perturbation  in  the  x  ory  transponder 
positions,  the  effect  on  the  initially  vertical  array  is  a  hori¬ 
zontal  translation  in  the  direction  of  the  transponder  with  no 
relative  arcov  positional  errors.  If  an  error  is  introduced  in 
the  z  transponder  position,  the  array  is  translated  and  tilted. 
Figure  5(a)  shows  the  effect  on  the  array  as  the  blue  '.ran¬ 
sponder  is  perturbed  by  100-m  errors  individually  in  x,  y, 
and  z,  and  finally  simultaneous  100-m  errors  ir.  x,  y,  and  z. 
Figure  5(b)  shows  the  effect  of  transponder  depth  errors. 
Each  transponder  depth  was  increased  and  decreased  indi¬ 
vidually  by  10  and  25  m.  The  dot  in  the  center  of  the  figure  is 
the  navigated  FLIP  and  array  positions  with  no  input  errors; 
with  a  decrease  in  transponder  depth,  the  input  slant  range 
appears  to  be  too  long,  and  the  array  is  translated  in  a  direc¬ 
tion  away  from  the  transponder  (notated  by  O’s  on  the  plot) 
and  tilted  slightly.  An  increase  in  transponder  depth  oper¬ 
ates  similarly  except  in  a  direction  toward  the  transponder 
(notated  by  X 's  on  the  plot) .  The  results  of  the  simulations 
for  the  positive  perturbations  with  an  initially  straight  verti¬ 
cal  array  are  compiled  in  Table  II.  The  columns  indicate 
which  parameter  has  been  perturbed  and  the  amount  of  the 
perturbation,  the  translation  and  percent  change  in  range 
calculated  for  the  FLIP  position,  system  tilt,  and  position 
rms  [-/F  in  Eq.  (1)].  The  rms  error  is  shown  as  a  single 
number  if  the  difference  in  erroi  across  the  array  was  less 
than  0.1  m;  when  two  numbers  are  shown,  they  indicate  the 
error  at  FLIP  and  the  bottom  element  of  the  array.  By  exam¬ 
ining  the  FLIP  translation  and  rms  error  for  the  x  and  y 
perturbations  in  Table  II,  it  is  clear  that  the  x  position  is 
controlled  by  the  red  and  blue  transponders,  and  they  posi¬ 
tion  is  controlled  by  the  green  and  blue  transponders  (Fig. 
4).  The  difference  in  rms  error  between  the  three  transpon¬ 
ders  reflects  the  percent  deviation  from  original  slant  range. 
The  dispersion  in  rms  error  was  used  to  iterate  the  transpon¬ 
der  depth.  This  can  be  seen,  for  example,  in  the  results  of  the 
red  depth  perturbation  which  yields  an  rms  error  at  FLIP  of 
1 8. 1 3  m,  monotonically  decreasing  to  1 3.48  m  at  the  bottom 
of  the  array.  The  array  and  FLIP  positions  are  more  sensi¬ 
tive  to  errors  in  transponder  depth  (60%-70%  of  perturba¬ 
tion)  than  to  errors  in  transponder  horizontal  (xy)  position 
(  <35%  of  perturbation ).  To  demonstrate  this  sensitivity 
dichotomy,  with  FLIP  at  a  known  position  x,y  from  a  tran¬ 
sponder,  the  FLIP  translation  T.  due  to  an  error  <5.  in  tran- 


2350  2400  2450 

meters 

FIG.  5.  Error  simulation  results.  The  effect  of  enors  in  the  transponder 
positions  on  the  shape  and  position  of  FLIP  and  the  array  are  shown  in  plan 
view,  (a)  Effect  on  the  FLIP  and  the  array  as  the  blue  transponder  is  per¬ 
turbed  inx.y,  and  z.  (b)  FLIPandarraypositionsasafunction  of  transpon¬ 
der  depth  error. 


sponder  depth,  and  the  translation  Txy  due  to  errors  8X  and 
8y  in  transponder  xy  position  are 

T.=h—ssin{cos~'[(d  +  8!)/s]},  (9) 

Txy  =  h—  [(x  —  8x)2  +  (y  —  Sy)2]’/2,  (10) 

where  h  and  s  are  the  horizontal  and  slant  ranges,  respective¬ 
ly  (Fig.  4),  d  is  the  difference  in  depth  between  the  base  of 
FLIP  and  the  transponder,  x  is  the  component  of  h  in  the  x 
direction,  and  y  is  the  component  of  h  in  the  y  direction. 
FLIP  is  constrained  by  the  experimental  configuration  to 
move  approximately  along  the  line  defined  by  the  horizontal 
projection  of  the  slant  range  between  the  perturbed  tran¬ 
sponder  position  and  FLIP  (Fig.  5).  To  quantify  these  trans¬ 
lational  errors,  a  8X  of  +  25  m  in  the  red  transponder  depth 
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TABLE  II  Navigation  error  simulations 


Perturbation 

(m) 

%  change 
in  slant 
range 

%  change 
m  horizontal 
range  projection 

Array  tilt 
(deg) 

Translation 
from  initial 
position 
(m) 

rms 

error 

(m) 

<+ of 
original 
perturbation 

FLIP  and  array  position  with  errors  in 

transponder  positions 

Red  X  +25 

0.1323 

0.7780 

0.000 

15.802 

7  11 

28.4 

Green 

-00693 

-  0.3095 

Blue 

-  0  0724 

-0  2940 

Red 

00110 

0  0640 

0.000 

1  498 

0.78 

3  1 

Green  X  +  25 

-  0.0124 

-0  0550 

Blue 

00012 

0  0049 

Red 

0  0576 

0  3333 

0000 

14.161 

6  45 

25.8 

Green 

0.0636 

0  2832 

Blue  X+  25 

-0  1360 

-  0.5463 

Red  T+25 

-  0  0638 

-  0.3663 

0.000 

7.540 

3.83 

153 

Green 

0.0383 

0  1708 

Blue 

0.0296 

0  1198 

Red 

-  0.0669 

-  0.3880 

0.000 

15  849 

8.73 

34.9 

Green  T+25 

0.1468 

0.6636 

Blue 

-  0.0879 

-  0.3567 

Red 

0.0454 

0.2627 

0.000 

11  149 

5.01 

20.0 

Green 

0  0499 

0.2223 

Blue  T+  25 

-  0  1069 

-0  4304 

Red  2+25 

-  0  3083 

-  1.8163 

0.482 

37.199 

18.13,  13.48 

72  5,  53  9 

Green 

0  1778 

0  7907 

Blue 

0  1624 

0.6567 

Red 

0.1405 

08112 

0  375 

29.320 

16.63,  12.39 

66.5,  49.6 

Green  Z  +  25 

-  0.2689 

-  1.2156 

Blue 

0  5489 

0  6022 

Red 

0  1225 

0.7076 

0  382 

30  631 

14.46,  10.76 

57.8, 43  0 

Green 

0  1419 

0.6311 

Blue  Z  +  25 

-0  2901 

-  1  1915 

FLIP  and  array  positions  with  errors  m  FLIP  slant  ranges 

Red  +25 

0  3561 

2  0453 

0.424 

41.936 

1904.  15  73 

76  2,  62.9 

Green 

-0.1865 

-  0  8346 

Blue 

-0.1855 

-  0.7543 

Red 

-0.1454 

-  0.8450 

0.315 

33.396 

18.28,  15.20 

73.1,60.8 

Green  +  25 

0.3122 

1.3848 

3lue 

-0.1780 

-  0.7236 

Red 

-0  1279 

-  0  7432 

0315 

34.034 

16  08,  13.51 

64.3,  54.0 

Green 

-0  1577 

-  0  7054 

Blue  +  25 

0.3284 

I  3250 

of  Eq.  (9)  resulted  in  a  58-m  FLIP  xy  position  translation, 
and  a  6X  of  +  25  m  in  the  red  transponder*  position  of  Eq. 
L 10 J  resulted  in  a  FLIP*;  position  translation  of  only  22  m, 
illustrating  that  FLIP  position  is  more  sensitive  to  errors  in 
transponder  depth  than  to  errors  in  transponder  xy  position. 

An  error  in  transponder  position,  when  the  arra>  is  tilt¬ 
ed,  will  distort  the  vertical  and  azimuthal  array  orientation 
as  well  as  impart  a  horizontal  translation.  A  tilted  arra>  will 
show  a  horizontal  translation  comparable  to  the  initially 
straight  array,  but  with  az:muthal  and  vertical  rotations. 
The  vertical  angular  change  is  similar  to  that  of  the  straight 
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array  for  depth  perturbations  and  small  but  "oticeable  for 
horizontal  perturbations.  The  azimuthal  angular  change  is 
larger  for  both  horizontal  perturbations  (0.2°  tilt)  and  depth 
perturbations  (14° tilt). 

An  error  in  a  FLIP  slant  range  also  translates  into  an 
array  positional  error.  An  error  •vhic'i  increases  the  slant 
range  moves  FLIP  away  from  the  transponder  and  the  array 
closer.  It  also  induces  a  tilt  in  the  array  making  the  bottom 
closer  to  the  no-error  position  than  the  top,  with  a  corre¬ 
sponding  decrease  in  rms  error  from  the  top  to  bottom  of  the 
array.  Table  ll  shows  the  results  of  the  simulation  for  +  25- 
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m  errors  in  the  three  transponder  slant  ranges  to  FLIP.  A 
transponder  slant  range  error  <5,  before  iteration  creates  a 
translation  error: 

Ts  =  h  —  yjhi  —  2Ss  s  —  S\,  (11) 

where  h  and  s  are  denned  above.  An  error  <5,  of  +  25  m  m 
the  red  transponder  slant  range  of  Eq.  ( 1 1 )  yields  a  transla¬ 
tion  error  T,  of  62.3  m  prior  to  iteration. 

These  simulations  quantify  the  mapping  of  errors  m 
slant  ranges  to  errors  in  xy  positions.  The  adjustment  direc¬ 
tion  is  dictated  by  the  measurement  direction  (horizontal 
range  direction)  not  the  model  parameters  (xy  directions) 
or  the  error  direction.  The  positions  of  FLIP  and  the  array 
are  more  sensitive  to  errors  in  transponder  depth  and  FLIP 
slant  range  than  to  errors  in  transponder  horizontal  position. 
This  sensitivity  is  attributed  to  the  constraint  of  fixed  array 
depth  (negligible  measurement  errors)  necessitating  an  as¬ 
signment  of  the  error  to  the  horizontal  projection.  These 
errors  are  magnified  and  dispersed  across  the  array.  It  was 
this  dispersion  that  was  used  as  a  criteria  for  convergence 
during  the  transponder  depth  iteration. 

III.  DATA  ANALYSIS 

The  motions  of  FLIP  and  the  array  are  analyzed  to  de¬ 
termine  the  navigation  accuracy,  the  extent  of  movement, 
and  possib'e  forcing  functions  driving  the  system.  The  data 
consisted  of  the  GPS  satellite  positions  acquired  at  FLIP, 
estimated  FLIP  positions  spanning  a  period  of  18  days  with 
1-h  samples,  and  estimated  array  element  positions  at  a  nom¬ 
inal  depth  of  850  m  spanning  a  24-h  period.  The  GPS  posi¬ 
tions  acquired  aboard  FLIP  provide  an  independent  verifi¬ 
cation  of  the  estimated  FLIP  positions  derived  using  the 
least-squares  method  ( Fig.  6 ) .  The  GPS  positions  in  latitude 
and  longitude  were  converted  to  meters  from  the  least- 
squares  origin  by  calculating  the  local  radius  of  the  earth  as 
follows13: 

*gps  =  'pm  <5(lat),  yGPS  =  r ^  cos(lat)  <5(long), 

(12) 

where 

^pos  ='■«,[!  — /sin:  (lat) ],  /= 

rc0 

and  where  xGPS  is  the  E-W  distance,  yGPS  is  the  N-S  dis¬ 
tance,  is  the  earth’s  radius  at  the  position  latitude,  S (lat) 

is  the  difference  in  latitude  between  the  position  and  the  ori¬ 
gin,  5(long)  is  the  difference  m  longitude,  iat  is  the  latitude 
of  the  position,  is  the  earth’s  equatorial  radius,  and 
is  the  earth’s  polar  radius.  The  observed  10-m  rms  deviation 
between  the  GPS  and  navigated  FLIP  positions  is  reassur¬ 
ing.  FLIP  positional  errors  are  independent  and  distributed 
normally  at  a  x2  0-05  level  of  significance  with  zero  mean 
and  a  standard  deviation  of  1.1  m.  Array  positional  errors 
may  deviate  from  a  normal  distribution,  depending  on  the 
nature  of  the  transponder  and  receiver  errors;  however,  con¬ 
sidering  a  receiver  with  low  errors,  a  Gaussian  distribution 
of  2.9-m  mean  and  0.8 1-m  standard  deviation  may  be  com¬ 
pared  at  a  0.05  level  of  significance.  Receiver  detection 
threshold  mismatch,  tuned  filter  drift,  sound-speed  errors, 
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FIG  6.  Comparison  of  GPS  satellite  and  acoustic  navigation  FLIP  GPS 
positions  are  compared  to  those  estimated  by  the  least-squares  navigation. 
The  large  symbols  represent  the  GPS  positions,  while  the  corresponding 
small  symbols  represent  the  estimated  FLIP  positions.  A  bias  correction, 
attributed  to  drift  of  the  transponders  during  deployment,  of  24  5  m  was 
added  to  the  N-S  estimated  positions  and  1  2  m  to  the  E-W  estimated  posi- 
uons. 


FIG  7  Errordistnbutions.Theenordistnbutionsshownfor  (a)  FLIPar.d 

v  b )  array  navigation  receiver  No  8  are  independent  as  illustrated  by  a  runs 
test  and  Gaussian  at  a  0.05  level  of  significance  The  high  transponder  errors 
shown  in  the  6-m  bin  of  (b)  are  due  to  transponder  detection  of  multipath 
propagation 
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the  FLIP  position,  as  described  earlier,  contribute  to  the 
mean.  FLIP  and  array  error  distributions  are  shown  in 
Fig.  7. 

The  motion  of  FLIP  (Fig.  8)  is  primarily  controlled  by 
the  local  wind  field,  the  tides,  and  the  tension  on  the  mooring 
lines.  The  low-frequency  movement  of  FLIP  is  driven  by  the 
local  wind,  primarily  from  the  north  during  the  entire  exper¬ 
iment.  The  wind  speed  has  a  high  negative  correlation  with 
FLIP  N-S  positions  and,  to  a  lesser  extent,  with  FLIP  E-W 
positions.  This  is  expected  as  a  northerly  w  ind  w  ould  tend  to 
push  FLIP  toward  the  south.  At  35°  N  latitude  the  tidal  cycle 
is  semidiurnal  and  appears  as  a  prominent  component  in  the 
FLIP  time  series.  The  tidal  component  in  the  E-W  displace¬ 
ment  lags  2  h  behind  the  N-S  component,  producing  a  clock¬ 
wise  rotation.  This  rotation  is  also  seen  in  the  calculated  tidal 
velocity  ellipse  (Fig.  9),  which  lags  the  displacement  by  90°. 
The  tidal  velocity  ellipse,  which  includes  contributions  from 
both  the  barotropic  and  baroclinic  modes,  is  comparable  to 
simultaneous  estimates  of  the  current  field  recorded  by  cur¬ 
rent  meters  deployed  at  50  and  100  m  below  the  water  sur¬ 
face  from  FLIP.  Although  the  current  meter  measurements 
are  relative  to  FLIP,  the  tidal  component  is  four  times  as 
strong  as  the  constrained  FLIP  component  and  leads  in 
phase  due  to  FLIP’S  massive  structure  and  mooring  line  ef¬ 
fects. 

The  motion  of  the  array  (Fig  10)  is  controlled  by  the 
tides,  the  motion  of  FLIP,  and  possibly  by  internal  waves  or 
higher  frequency  surface  motion  that  is  unresolved  in  the 
FLIP  data  set  The  longest  penod  movement  of  the  array 
was  on  the  order  of  days  and  was  driven  by  FLIP'S  response 


to  the  wind.  This  is  inferred  by  the  fact  that  the  top  array 
position  was  normally  within  a  30-m  horizontal  slant  range 
of  FLIP  which  traveled  over  300  m  during  the  1 8  days  (with¬ 
in  20  m  in  Fig.  11).  The  tidal  cy  cle  is  ev  ident  as  a  1 2-h  oscilla¬ 
tion  in  the  24-h  time  series  and  as  the  0.083-cy  cles/h  peak  in 
the  spectrum  (Fig  12). 

Higher  frequency  array  oscillations  appear  in  the  ex¬ 
panded  time  series  in  Fig  10(c)and  (d),  but  are  not  station¬ 
ary  over  the  24  h,  and  ail  but  disappear  from  the  spectrum 
during  the  averaging  process  Although  the  sample-to-sam- 
ple  motion  of  the  array  is  on  the  seme  order  as  the  errors,  the 
errors  are  independent  and  Gaussian,  therefore,  we  believe 
the  fluctuations  which  occur  with  periods  from  1 5  or  20  min 
to  2  h  are  due  to  array  motion  Modeling  the  array  as  a 
simple  pendulum  where  the  center  of  mass  is  derived  from 
four  component  parts  (uplink  wire,  array,  kevlar,  and 
weights),  the  natural  period  of  oscillation  is  just  over  1  min 
The  damped  natural  period  for  a  velocity  of  5  cm /s  is  7  min 
with  a  damping  coefficient  of  slightly  less  than  critical 
These  periods  do  not  coincide  with  the  0.3-2-h  periods  seen 
in  the  array  data,  therefore,  the  high-frequency  motions 
must  be  due  to  environmental  forcing  functions  We  believe 
it  to  be  a  combination  of  two  processes,  the  coupling  of  the 
array  to  the  surface  and  the  array  response  to  low -order 
modes  of  the  internal  wave  field.  The  motion  whose  periods 
are  from  tens  of  minutes  to  many  hours  ( Vaisala  frequency 
in  the  N  Pacific  is  about  20  min  near  the  surface)  and  whose 
horizontal  velocities  of  1-2  cm  A  are  reasonable  for  a  con¬ 
strained  internal  wave  response1*  where  typically  horizontal 
currents  are  about  5  cmA  and  horizontal  current  spectra  is 
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FIG.  6.  FLIP  posil.on  time  senes.  The  FLIP  positions  are  estimated  approximate]}'  ever}  hour  and  plotted  over  the  18-day  experiment  The  (a )  .r  positions, 
increasing  toward  the  cast,  and  the  (bj  y  positions,  increasing  toward  the  north,  are  plotted  on  the  same  rcbti'  a  vertuw..  ,,calc,  but  with  different  on6.n 
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FIG.  9  Tidal  current  ellipse  (a)  The  cur¬ 
rent  ellipse  estimated  from  FLIP  positions 
was  calculated  by  finite  differencing  the 
18-day  time  senes,  low-pass  filtenng  the 
resulting  hourly  estimates,  and  then  aver¬ 
aging  12-h  segments  The  tips  of  the  cur¬ 
rent  vectors  were  plotted  from  a  common 
origin  every  hour.  The  0  GMT  vector  is 
marked  in  cm/s.  (b)  Current  meter  data 
collected  from  FLIP  represent  !2-h  seg¬ 
ment  averages  (O)  over  13  days  at  43  m 
and  ( * )  over  20  days  at  90  m.  Mean  veloc¬ 
ities  of  7.5  cm/s  E,  0.64  cm/s  N  ct  the  43- 
m  site,  and  7.2  cm/s  E.  -  0.9 1  cm/s  N  at 
the  90-m  site  were  removed. 
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dominated  by  the  lower-order  modes.17  The  motion  spec¬ 
trum  does  not  show  a  peak  at  the  local  inertial  frequency,18 
nor  an  energy  dependency  between  the  inertial  frequen¬ 
cy  (0.048  cycles  A)  and  the  Vaisala  frequency  common  in 
internal  wave  energy  spectra.19  This  characteristic  shape  is 
apparently  filled  in  by  the  coupled  response  to  surface-gener¬ 
ated  motions. 

The  array  shape  during  the  experiment  is  modeled  to 
first  order  as  linear,  but  slightly  tilted  from  vertical.  The 
array  shape  plotted  approximately  every  3  min  over  the  first 
hour  of  data  shown  in  Fig.  10(c)  and  (d)  is  illustrated  in 
Fig.  1 1,  where  the  horizontal  deformation  has  been  exagger¬ 
ated  for  clarity,  extending  50  m  in  the  horizontal  compared 
to  1 600  m  in  the  v  ertical.  These  data  show  a  linear  array  with 
receiver  position  deformations  within  the  rms  error  estimat¬ 
ed  by  the  least-squares  procedure4;  and  with  a  deviation  of 
less  than  5  m  across  the  array  aperture,  the  array  maintains  a 
vertical  orientation  to  within  0.5°.  Array  shape  snapshots 
every  5  h  over  a  longer  period  of  time  (20  h)  show  that  the 
array  follows  FLIP  to  within  a  30-m  horizontal  range  (Fig. 
13).  The  change  in  array  depth  due  to  such  a  30-m  deviation , 
assuming  a  linear  cable,  is  only  1  m,  which  is  not  significant. 
The  array  tilt  from  vertical  is  withir  about  T,  considering 
the  maximum  deviation  (15  m)  shown  across  the  aperture. 
The  linear  pendulumlike  response  is  particularly  apparent  in 
the  2-min  snapshots  as  the  array  responds  to  the  motion  of 
FLIP  and  the  internal  wave  field.  Although  the  time  se¬ 
quence  of  the  individual  estimates  is  not  legible,  the  crossing 
uplink  lines  indicate  small-scale  array  oscillations,  through 
which  the  linear  shape  is  maintained.  It  is  evident  that  the 
array  exhibits  a  more  active  response  (presumably  due  to 
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FIG  1 1  Array  shape  over  1  h  The  positions 
of  the  ten  deepest  navigation  receivers  are 
shown  over  a  1-h  period  The  top  two  receiver 
positions  were  deleted  due  to  their  large  errors 
(Ref  4)  The  two  views  are  (a)  East- Wes! 
versus  depth,  and  (b)  North-South  versus 
depth  The  snapshots  are  spaced  at  approxi¬ 
mate  3-mm  intervals  and  represent  a  subset  of 
the  data  shown  in  Fig  10(c)  and  (d) 
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2390  2400  3180 


3.a0  3200 

meters 


interna]  waves)  than  FLIP  by  examining  the  extent  of  FLIP 
motion  compared  to  that  of  the  array  during  the  same  time 
period.  The  array  moves  over  four  times  the  E-W  distance 
covered  by  FLIP  during  the  hour  plotted  m  Fig  1 1 


I  (a)  Receiver  8  \  -  portion 
I 


-80 - * - - - * 

0  2  4  6  8  10  12  14  16  18  20 

CYCLES  /HOUR 

FIG  12  Spectrum  of  array  motion  For  receiver  (a)  x  position  and  (by  > 
position,  the  smv.au  curve  is  an  averaged  1024-pi  spectrum  of  15,  50G 
overlapped,  Kaiser -Bessel  windowed  (o  —  2  5)  12S-point  time  series  plot¬ 
ted  on  a  relative  scale  A  single  spectrum  was  computed  similarly  over  the 
entire  time  series  and  plotted  on  the  same  relativ  e  scale  with  0  039-cyclcs/h 
resolution  for  more  aveurat  identification  of  the  narrow-band  spectral 
components 


IV.  CONCLUSION 

The  least-squares  method  was  implemented  to  convert 
measured  trav  el  time  data  to  spatial  positions  for  FLIP  and  a 
900-m  vertical  array  by  assuming  a  harmonic  sound  speed 
and  normally  distributed  errors.  Navigated  FLIP  positions 
(mean  =  0.5  m,  a  =  1.1  m)  agree  with  positions  obtained 
from  GPS  satellite  navigation  to  within  an  rms  error  of  10  m. 
A  mean  error  as  small  as  2.8  m  with  a  standard  deviation  of 
0.8  m  was  ob  ained  for  array  element  position;  however, 
each  receiver  displayed  a  unique  error  distribution  due  to 
variations  in  the  array  detectors.  Array  positions  were  found 
to  be  particularly  sensitive  to  errors  in  FLIP  position  and 
transponder  depth.  This  was  demonstrated  by  simulation  to 
be  a  result  of  the  optimization  constraints.  Iteration  of  array 
depths  and  incorporation  of  a  more  accurate  sound-speed 
correction  has  shown  promising  preliminary  results  by  re¬ 
ducing  the  mean  rms  array  element  error  to  less  than  1  m. 
while'  leaving  the  estimated  positions  within  0.2  m  of  the 
positions  presented  here.  The  navigated  positions  are  tied  to 
an  absolute  frame  of  reference  using  transponder  GPS  loca¬ 
tions,  and  the  absolute  positions  are  substantiated  by  the 
GPS  positions  at  FLIP  and  the  high  negative  correlation 
between  FLIP’S  N-S  component  and  the  w'ind  speed.  FLIP’S 
position  was  driven  by  the  winds  and  tides,  but  was  con¬ 
strained  to  approximately  a  300-m  range  by  a  three-point 
mooring  The  array  position  was  driven  by  tidal  motions  at  a 
semidiurnal  period,  and  apparently  by  internal  waves  and 
coupled  surface- wave  motion  at  higher  frequencies.  Array 
high-frequency  motions  were  up  to  a  factor  of  4  larger  than 
FLIP  motions,  although  the  array  normally  maintained  a 
vertical  orientation  to  within  a  tilt  angle  of  2°  and  remained 
within  a  30-m  range  of  FLIP. 
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2  •  1011 

0  -  00-47 

3  -  15-34 

1  -  05—1 

4  -  1943 

FIG.  U  Array  shape  over  20  h.  The  positions  of  the  !2navigation  receivers  are  shown  at  5-h  intervals  over  a  20-h  period  The  two  views  are  (a)  East-West 
versus  depth  and  (b)  North-South  versus  depth.  Snapshots  every  2  nun  are  shown  beginning  at  5:31  GMT  over  a  24-min  period. 
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Abstract — After  26  years  of  operations  at  sea  with  the  Research 
Platform  FLIP,  a  355-ft-long  craft  which  has  a  draft  of  300  ft  in  the 
vertical  position,  planning  is  underway  for  a  larger,  more  capable  FLIP  II 
that  can  operate  in  seas  with  waves  of  up  to  80  ft  and  survive  more  severe 
seas.  As  our  research  in  acoustics  extends  to  studies  of  ambient  noise  and 
propagation  at  lower  frequencies,  recent  deployments  of  large  aperture 
(900  m),  multielement  (120)  arrays  tax  the  space  and  facilities  on  FLIP. 
Future  plans  include  deployment  of  240  element,  3500  m  arrays  as  well  as 
large,  powerful  sound  sources  that  further  tax  the  limits  of  current 
capabilities. 

The  original  FLIP,  while  designed  to  work  in  30*ft  waves  with  no  more 
than  a  18  percent  heave  response,  has  survived  80-ft  swells  (22  s  period) 
with  only  minor  damage.  Research  requirements  for  greater  laboratory 
and  deployment  capabilities  to  make  possible  joint  multi-disciplinary/ 
group  operations  and  the  need  to  work  in  regions  with  worse  weather 
cannot  be  met  with  the  current  FLIP  nor  is  it  economically  or  structurally 
feasible  to  modify  it  to  do  so.  Aiong  with  a  brief  description  of  FLIP,  this 
paper  lists  the  variety  of  research  conducted  from  FLIP  in  both  the 
Pacific  and  Atlantic  Oceans  and  whether  it  was  in  drifting  or  moored 
deployments  in  deep  water.  While  originally  built  for  acoustic  research,  it 
has  also  been  used  for  physical  oceanography,  meteorological,  and 
biological  studies.  A  recent  workshop  on  research  platforms  highlighted 
the  needs  of  several  different  scientific  communities  for  a  larger,  more 
capable  FLiP-type  vessel  with  a  minimal  cir  and  water  column  profile 
that  can  work  in  severe  weather. 

The  purpose  of  this  paper  is  to  invite  the  interest  and  attention  of  the 
scientific  and  engineering  communities  regarding  the  design  and  potential 
uses  of  a  larger  FLIP  II.  Current  thinking  centers  ground  420-ft-long,  30- 
ft-diameter  hull  with  a  smaller  water  plane  diameter  of  20  ft.  All  weather 
deployment  capability  of  submersibles  from  the  bottom  of  FLIP  II  (draft 
340  ft)  is  considered.  Preliminary  design  efforts  have  been  initiated  with  a 
naval  architect. 

I.  Introduction 

FTER  26  years  of  operations,  the  U.S.  Navy’s  deep  draft 
(300  ft)  Research  Platform  FLIP  (FLoating  Instrument 
Platform)  {1)  shown  in  Fig.  1  is  continuing  to  provide  a  unique 
stable  platform  from  which  to  conduct  research  at  sea  without 
suffering  the  violent  motions  that  accompany  storms  and 
rough  weather.  In  fact,  because  of  the  deep  draft  and  long 
heave  period  of  27  s,  those  working  on  FLIP  during  the  sea- 
state  six  conditions  can  enjoy  charcoal  broiling  outdoors.  The 
stability  minimizes  adverse  effects  on  deployed  instrumenta¬ 
tion  as  well  as  fatigue  for  those  personnel  on  board.  It  also 
makes  possible  experimental  research  that  can  be  done  in  no 
other  way. 

Manuscript  received  July  15, 1988;  revised  September  26, 1988.  This  work 
was  supported  by  the  Office  of  Naval  Research  under  Contracts  N 00014-80- 
C-0220  and  N00014-87-C-OI27. 
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The  355-ft-long  spar  buoy,  designed  by  the  Marine  Physical 
Laboratory  (MPL)  of  the  Scripps  Institution  of  Oceanography 
for  acoustic  research  [1],  [2],  has  been  operating  for  26  years 
since  it  was  launched  June  22,  1962,  at  the  Gunderson 
Brothers  Engineering  Company  in  Portland,  Oregon.  The  low 
heave  response  (max.  18  percent  at  17  s)  in  heavy  weather  and 
its  low  tilt  fluctuations  [2]  more  than  met  the  stability 
requirements  for  the  immediate  acoustic  research  objectives. 
In  addition,  the  design  was  intended  to  make  FLIP  a  platform 
that  would  be  useful  fo.‘  other  types  of  research  and  that,  in 
fact,  has  been  the  case  (as  seen  in  Table  I). 

The  unique  characteristics  of  a  FLIP-type  platform  (stabil¬ 
ity,  low  profile  in  water  and  air  column,  ease  of  mooring  in 
deep  ocean)  have  attracted  the  attention  of  scientists  in  several 
different  disciplines  who  are  interested  in  a  larger,  more 
capable  FLIP  II.  Physical  oceanographers,  meteorologists, 
biologists,  and  marine  physicists  studying  acoustic  propaga¬ 
tion  and  ambient  noise  in  the  ocean  participated  in  the  Platform 
Workshop  held  at  the  MPL  on  June  29-30,  1988.  Their 
specific  needs  were  spelled  out  in  considerable  detail  regard¬ 
ing  requirements  for  a  larger  FLIP.  The  workshop  was 
instigated  by  recent  requests  to  participate  in  research  in  more 
northemly  waters  and  by  the  need  to  carry  increasingly  larger 
scientific  parties  with  correspondingly  larger  amounts  of 
equipment. 

For  example,  at  MPL  we  have  been  studying  ambient  noise 
and  acoustic  propagation  in  the  ocean  over  the  last  two 
decades,  making  use  of  vertical  arrays  of  hydrophones 
deployed  from  FLIP  in  three-point  moors,  most  recently  with 
a  900  m,  120-element  array,  as  shown  in  Fig.  2.  For  studying 
low-frequency  ambient  noise  and  propagation  we  need  to  use 
large  hydrophone  arrays  to  provide  the  resolution  to  separate 
multipaths  and  to  study  the  differences  between  local  and  long 
distance  sources  of  ambient  noise,  man-made  or  natural.  The 
three-point  mooring  assures  a  minimal  flow  noise  around 
hydrophones  and  a  maximum  vertically  of  the  arrays.  In  order 
to  pursue  these  studies  in  regions  which  are  more  subjected  to 
storms,  and  with  larger  and  larger  arrays  to  study  lower 
frequencies,  we  find  a  recurring  need  for  a  larger,  more 
capable  replacement  for  FLIP. 

A  second  major  area  of  research  at  MPL  which  utilizes 
FLIP  has  been  the  study  of  internal  wave  propagation  in  the 
upper  ocean.  This  work  was  initially  performed  using  both 
fixed  and  vertically  profiling  temperature  sensors.  These  were 
mounted  on  horizontal  booms  rigged  from  FLIP  to  provide  a 
30  m  horizontal  array.  The  profiling  system  obtained  data 
from  the  surface  down  to  400  m  on  a  2-min-cycle  time.  The 
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Fig.  I.  FLIP  at  a  draft  of  295  ft  with  booms  not  yet  folded  out.  Boston 
Whaler,  on  starboard  side,  ts  used  to  transfer  equipment  and  people. 


original  profiling  system  was  developed  in  1969  and  has  been 
extensively  upgraded  several  times.  At  present,  electrical 
conductivity  is  sensed,  as  well  as  temperature  and  depth, 
enabling  the  estimation  of  the  oceanic  salinity  and  density 
fields.  Or  a  recent  cruise,  PATCHEX  (1986).  ten-thousand 
density  profiles  were  obtained  over  a  three-and-a-half  week 
period. 

The  profiling  system  has  now  been  augmented  by  MPL- 
developed  Doppler  sonars.  Scattering  high-frequency  sound 
off  zooplankton  in  the  sea.  the  sonars  are  capable  of  sensing 
ranges  in  excess  of  1 .5  km,  with  a  1  cm/s  precision  and  a  10  m 
rar.ge  resolution.  Along  with  the  narrow  beam  deep-sea 


sonars,  s  new  generation  of  fan-beam  devices  is  being 
developed  to  scatter  from  the  underside  of  the  sea  surface. 
These  devices  can  measure  surface  wave  propagation  as  well 
as  lower  frequency  surface  currents  with  a  precision  unattain¬ 
able  by  convention'’!  means. 

The  Dop-^  ’  atv.  .stic  wor’<-,  as  well  as  analogous  radar 
sensing  progiams,  piace  a  premium  on  the  pitch/roll  stability 
of  FLIP,  rather  than  the  heave  response.  It  does  little  good  to 
construct  narrow  beam  systems  and  then  mount  them  on  a 
gyrating  platform.  Increased  tilt  stability  will  be  a  design 
objective  of  the  new  platform. 

Meteorological  experiments  making  use  of  FLIP  exploring 
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TABLE  I 

TYPES  OF  RESEARCH  UTILIZING  FLIP 


Principal 

Project  Investigator 


OCEAN  ACOUSTICS/UNDERWATER  SOUND 


Phase  and  amplitude  fluctuauons 

F.  H.  Fisher 

Crustal  anisotropy 

G.  G.  Shor 

Ambient  noise 

V.  C.  Anderson 

F.  N.  Spiess 

G.  B,  Morris 

R.  C.  Tycr 

F.  H.  Fisher 

W.  S.  Hodgkiss 

J.  A  Hildebrand 

Sea  surface  noise 

V.  C.  Anderson 

Sound  propagation 

F.  N  Spiess 

J.  Northrop 

G.  B.  Morris 

F.  H.  Fisher 

N.  Booth 

D.  A.  Ramsdale 
W.  S.  Hodgkiss 
J.  A.  Hildebrand 

Eottom-bounce  propagation 

F.  H.  Fisher 

Coherent  recombination  of  multipaths 

W.  S.  Hodgkiss 
R.  Brienzo 

Attenuation  in  sediments 

R  Brienzo 

Biological  acoustic  scattering  &  taxonomy 

P.  Greer.blatt 

Sound  absorption  as  function  of  pressure 

H.  Bezdek 

Acoustic  bac.tscatter  from  surface  waves 

S.  McConnell 

H.  Medwin 

Point  versus  planar  scattering  from  density  layers 

V.  C.  Anderson 
G.  T.  Kaye 

Doppler  current  meter  (10  MHz) 

P.  Rudnick 

High-frequency  (90  kHz)  echo  sounder 

F.  Fisher 

Design  &  development  of  high-resolution 

R.  Pinkei 

Principal 

Project  investigate- 


PHYSICAL  OCEANOGRAPHY 

Storm  generated  waves 

W  H  Munk 

Interna!  waves  by  thermistor  yo-yos 

R  Pinkei 

K.  Zalkan 

Internal  waves  by  doppler  sonar 

R  Pinkei 

Langmuir  circulation  by  doppler  sonar 

R.  Pinkei 

J.  Smith 

Sea  surface  slope  distribution 

R  Pinkei 

Wave  direction  using  FLIP 

P.  Rudnidc 

Mixed  layer  dynamics/air-sea  fluxes  (MILDEX) 

R  A  Weller 

Langmuir  cell,  Ekman  circuhtions  (MILDEX) 

R.  A.  Weller 

VMCM  measurements 

R.  A.  Weller 

R.  Pinkei 

Ocean  optics  (ODEX) 

J.  J  Simpson 

Ocean  waves  (OWAX) 

C.  Friehe 

J.  J.  Simpson 

Ocean  natural  resources 

R.  Yodei 

Radar  backscatter  from  waves 

(u) 

Surface  wave  directional  spectra  (BOMEX) 

R  Davis 

L.  Regier 

Winding  profiling  between  wave  peaks  (BOMEX) 

(u) 

Air  turbulence,  MET  surface  waves  (POLE) 

R.  Davis 

C.  Friehe 

J.  Simpson 

L.  Regier 

C.  Paulson 

Turbulence  and  microstructure 

R.  B.  Williams 

air-sea  interactions  have  also  been  an  on-going  effort,  most 
notably  in  the  Barbados  Oceanographic  Meteorology  Experi¬ 
ment  (BOMEX)  in  the  Atlantic  Ocean  in  1969.  For  this 
experiment,  wind  tunnel  tests  conducted  by  Prof.  E.  Mollo- 
Christensen  at  the  Massachusetts  Institute  of  Technology 
showed  that  the  wind  profile  is  altered  by  less  than  1  percent  at 
a  distance  of  30  in  from  the  centerline.  However,  as  seen  in 
Fig.  3,  it  was  only  possible  to  support  the  servo-controlled 
wind  profiler  (constant  height  above  water)  at  a  distance  of 
about  15  m  from  the  center  of  FLIP.  A  larger  FLIP  would 
make  it  possible  to  extend  a  hydraulic  boom  farther  away  from 
the  hull  for  sensors  requiring  minimum  flow  contamination. 
Reduced  difficulties  in  holding  a  given  heading  duri’g  such 
work,  as  well  as  working  in  rougher  weather,  will  be  possible 
on  a  larger  FLIP. 

Other  notable  air  sea  interaction  experiments  are  POLE 
(1974),  MILDEX  (1983),  ODEX  (1984,,  and  PATCHEX 


(1986).  In  1982,  Dr.  B.  Weller  of  the  Woods  Hole  Oceano¬ 
graphic  Institution  measured  vertical  velocities  in  excess  of  1* 
cm/s  in  the  mixed  layer,  associated  with  particularly  intense 
Langmuir  cell  activity.  This  was  a  first  sign  of  the  enormous 
efficiency  of  Langmuir  cells  in  mixing  momentum  from  the 
wave  zone  into  the  interior  of  the  mixed  layer.  More  detailed 
studies  of  the  cells  themselves  occurred  in  MILDEX.  A  new 
experiment,  SWAPP  (1990),  is  planned  to  investigate  the 
interaction  of  surface  waves  and  Langmuir  cells.  Multiple 
laboratories  and  greater  power  on  a  larger  FLIP  would  make  it 
possible  to  plan  extended  multi-disciplinary  experiments, 
featuring  simultaneous  operation  of  oceanographic  and  me- 
ierological  sensors  as  well  as  radars. 

At  this  point  it  is  appropriate  to  emphasize  FLIP’S  utility  as 
a  platform  from  which  to  test  new  instruments  anri  develop 
new  obser  mal  techniques.  As  a  result  of  FLIP’S  extreme 
economy  .ation  it  is  not  necessary  to  convene  planning 
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Fig.  2.  Schematic  of  FLIP  in  three-point  moor  with  vertical  array  deployed. 
Sound  velocity  profile  is  also  shown. 


Fig.  3.  FLIP  in  Atlantic  Ocean  participating  in  the  Barbados  Oceanographic 
Meteorological  Experiment  (BOMEX).  Note  vertical  servo  controlled  wave 
following  wind  profiling  staff. 


meetings,  coin  acronyms,  and  appoint  steering  committees  to 
justify  funding  sufficient  to  leave  the  dock.  Early  evaluation  of 
s^nic  anemometers  and  hot-film  sensors  for  meteorology,  the 
vector  measuring  current  meter  (VMCM),  volume,  and 
surface  scattering  Doppler  sonars  were  all  performed  from 
FLIP.  Without  this  inexpensive  source  of  access  to  the  deep 
sea,  our  understanding  of  the  performance  limits  of  these 
devices  would  have  been  more  difficult  to  obtain.  A  larger 
FLIP  will  have  increased  test  capabilities  for  multiple  research 
groups,  thereby  facilitating  intercomparison  experiments. 

The  ability  to  moor  FLIP  in  deep  water,  originally 
developed  to  reduce  flow  and  strum  noise  on  hydrophone 
arrays,  has  made  it  possible  for  Dr.  Pinkel  to  extend  Doppler 


sonar  ranges.  Prior  to  this,  the  range  was  limited  by  noise 
from  a  hydraulically  driven  orientation  motor  to  keep  FLIP  on 
a  constant  heading  while  operating  in  the  drifting  mode.  In  a 
moor,  however,  there  are  such  heading  variations  as  weather 
and  current  changes  which  lead  to  the  requirement  on  a  larger 
FLIP  to  be  able  to  maintain  a  constant  heading  within  a  few 
degrees  while  in  a  three-point  moor. 

Some  of  the  operations  have  included  the  mooring  of  FLIP 
in  deep  water.  2000  fathoms  and  deeper,  in  one-,  two-  and 
three-point  moors,  a  technique  developed  at  MPL  and  refined 
over  the  years  to  become  a  routine  operation  [3],  [4].  For  these 
moorings,  10  tons  of  surplus  anchor  chain  and  a  750  lb 
Danforth  anchor  with  specially  modified  flukes  and  palms  are 
deployed  for  each  leg  of  a  1 .5-inch-diameter  power  braid  with 
a  scope  of  about  1.5:1. 

Although  designed  originally  for  operating  in  waves  of  up  to 
30  ft,  it  has  survived  80  ft  blue  water  swells  with  a  period  of  22 
s,  the  null  heave  response  period  of  FLIP  [5].  The  peak  swells 
struck  FLIP  during  the  night  and  came  to  within  2  ft  of  the 
upper  platform,  as  shown  schematically  in  Fig.  4.  No  injuries 
and  only  minor  damage  was  incurred  during  the  high  waves. 

In  the  most  recent  operation  with  a  large-aperture  vertical 
array,  120  elements  uniformly  spaced  over  900  m,  the 
acoustic  navigation  of  the  array  was  included  through  the  use 
of  transponders  on  the  bottom  which  were  surveyed  in  with  the 
GPS  satellite  navigation  system  [6J.  This  involved  the  same 
techniques  originally  developed  for  the  Deep  Tow  operations 
at  MPL  [7],  During  a  recent  operation  in  June,  1987,  with 
FLIP  moored  in  2400  fathom  water,  the  watch  circle  diameter 
for  FLIP  in  winds  of  up  to  30  knots  was  only  20  m. 

With  a  crew  of  five  to  six  and  with  scientific  parties  of  up  to 
twelve,  FLIP  is  towed  to  its  station  since  it  has  no  propulsion 
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Fig.  4.  FLIP  in  vertical  position  and  with  a  80  ft  wave  superimposed.  The 
high  waves  came  to  within  2  ft  of  the  uppermost  platform  at  night.  There 
was  minimal  damage  to  FLIP.  R.  Hasse  took  movies  of  waves  the  following 
day. 


of  its  own  other  than  a  60  hp  orientation  motor  for  the  vertical 
position.  With  funding  support  from  the  U.S.  Navy,  costs  to 
the  user  are  about  $l,400/day  exclusive  of  towing,  mooring, 
and  harbor  tug  costs.  For  Navy-related  research,  it  is 
sometimes  possible  to  obtain  sea-going  tug  services  from 
7,200  hp  T-ATF-type  fleet  tugs  such  as  the  USNS  Navajo  and 
USNS  Narragansett. 

To  the  greatest  extent  possible  the  design  of  FLIP  II  will  be 
to  keep  operating  costs  within  the  same  range  as  those  quoted 
above.  Because  of  the  larger  structure,  however,  added  costs 
will  be  incurred.  For  example,  the  current  cost  of  a  three-point 
moor  (about  $23  300)  may  be  larger  for  FLIP  II,  whose 
diameter  may  be  50  percent  greater  (30  ft)  and  whose  length 
may  be  about  20  percent  longer  (420  ft)  in  order  to  provide  the 
necessary  stability  and  workspace  spelled  out  in  the  Platform 
Workshop  Report  [8], 

A  brief  discussion  of  mooring  FLIP  at  sea  is  presented 
below,  along  with  future  plans  for  mooring  FLIP  jointly  with 
the  Ocean  Research  Buoy  (ORB),  also  developed  at  MPL,  in 
order  to  deploy  horizontal  and  vertical  arrays  simultaneously 
in  a  pseudo-Mills  cross  to  study  directional  properties  of 
ambient  noise. 

II.  Review  of  FLIP  Development 

The  need  for  a  stable  platform  arose  from  the  problem  of 
measuring  the  effects  on  long-range  sound  propagation  of 
horizontal  variability  in  the  ocean.  At  that  time,  in  the  late 
’fifties,  long  range  meant  around  30  nmi.  In  much  the  same 
way  that  fluctuations  in  the  atmosphere  make  the  stars  twinkle 
(optical  scintillation),  temperature  and  salinity  fluctuations  in 
the  ocean  cause  acoustic  scintillation.  If  temperature  and 
salinity  gradients  exist  over  substantial  regions,  they  will 
cause  bending  or  refraction  of  sound  waves  in  the  ocean.  Just 
as  a  stable  platform  is  necessary  for  a  telescope,  it  was 
necessary  to  make  use  of  a  stable  platform  to  make  similar 
measurements  in  the  ocean. 

Initial  attempts  were  made  to  do  this  using  the  Navy 


Electronics  Laboratory’s  (NEL)  research  submarine,  USS 
Bay  a  (AGSS-318).  Using  a  5  kHz  sound  source  at  a  depth  of 
300  ft  at  the  convergence  zone  range  of  about  30  nmi, 
hydrophones  on  the  two  50  ft  wings  of  the  Bay  a  were  used  to 
listen  to  the  5  ms  pulses  and  measure  the  acoustic  fluctuations 
in  phase  and  amplitude.  With  the  Baya  at  a  depth  of  300  ft,  F. 
Hale  and  H.  Westphal  of  NEL  (now  Naval  Ocean  Systems 
Center  (NOSC))  showed  the  author  how  to  find  convergence 
zones  where  the  propagation  loss  can  be  as  low  as  70  dB  at  30 
nmi  and  80  dB  at  100  nmi  [9]. 

To  measure  refraction  effects,  however,  was  more  difficult. 
This  led  to  flying  1 100  cult  Kytoons  from  the  sound  source 
ships  (USS  Rexburg  or  USS  Marysville,  also  with  NEL)  at 
altitudes  of  1OC0  ft,  with  a  200  W/s  xenon  flash  pulsing  at  1 
pps.  Launchirg  the  alloons  from  the  forward  gun  tub  was 
difficult  enough  in  25-  to  30-knot  winds  at  around  30°N  and 
120°W.  However,  to  get  a  good  bearing  with  the  scope 
comparable  to  the  accuracy  of  our  acoustic  data  proved 
virtually  impossible. 

Time  for  working  on  our  experiment  was  very  limited 
because  this  was  only  one  of  several  on  board,  which  meant 
that  if  nothing  leaked  we  might  acquire  a  few  minutes  of  data 
at  best  since  the  submarine  had  to  keep  way  on  for  control  and 
the  zone  was  only  a  nautical  mile  or  so  wide.  At  a  depth  of  300 
ft,  for  the  5  ms  pul  se  it  was  possible  to  separate  multipaths  so 
that  the  direct  path,  the  desired  one  for  measurement,  was 
unambiguously  clear. 

In  January,  1960,  I  reported  my  complaints  to  Dr.  F.  N. 
Spiess.  Director  of  the  MPL,  under  whose  cognizance  I  was 
attempting  to  do  these  experiments.  He  responded  immedi¬ 
ately  by  saying  that  A.  Vine  of  the  Woods  Hole  Oceano¬ 
graphic  Institution  had  long  been  advocating  the  upending  of  a 
submarine  in  order  to  make  a  stable  platform  for  ocean 
research.  After  quickly  determining  that  submarine  design 
didn’t  really  lend  itself  to  operating  in  a  vertical  mode,  we 
proceeded  to  investigate  the  building  of  such  a  platform  with 
the  300-ft  depth  and  a  diameter  of  20  ft  as  our  contraints  (the 
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depth  in  order  to  separate  acoustic  multi  paths  at  ranges  of 2000 
yd,  and  the  diameter  to  permit  berthing  at  the  NEL  pier). 

Dr.  F.  Rudnick  showed  how  to  improve  the  heave  response 
period  from  18  s  for  a  cylinder  of  uniform  cross  section  (same 
as  for  a  300  ft  pendulum)  to  27  s  by  a  60  percent  reduction  in 
the  water  plane  area  (12.6  ft  diameter),  as  shown  in  Fig.  5.  Dr. 
Rudnick’s  heave  calculations  were  verified  with  the  heave 
response  tests  conducted  in  the  Lockheed  missile  test  tank  on  a 
1/25  scale  model  with  an  abrupt  transition  at  mid-draft. 
Subsequent  measurements  in  the  ocean  showed  that  the  heave 
response  of  FLIP  agreed  with  the  theory  as  shown  in  Fig.  6. 

For  a  spar  buoy  of  constant  cross  section,  the  resonant 
period  is  given  by  Che  pendulum  formula  in  which  die  draft  of 
the  buoy  is  the  length.  If  one  modifies  the  pendulum  formula  to 
account  for  the  reduced  area,  it  can  be  generalized  to  the  ratio 
of  the  total  displacement  in  the  vertical  to  the  spring  constant 
of  the  '"ater  plane  area 

T=2rrt[W/gk]  (1) 

where  W  =  displacement  in  tons,  k  =  spring  constant  of 
water  plane  area  (tons/ft),  and  g  =  gravitational  acceleration 
(ft/s/s). 

Using  Rudnick’s  theory,  we  selected  an  area  reduction  of  60 
percent,  with  the  transition  to  begin  above  the  150  ft  draft  for 
the  full-scale  FLIP.  This  gave  a  period  of  T  =  27  s,  well  away 
from  the  spectrum  of  wave  energy  in  the  ocean;  the  present 
configuration  has  a  null  heave  response  at  22  s.  Using  (1),  a 
value  of  the  parameter  of  W/k  equal  to  or  exceeding  591  will 
keep  T  at  27  s  or  greater. 

The  original  configuration  of  FLIP,  after  working  with  the 
naval  architects  Rosenblatt  and  Sons  and  subsequendy  with 
Glosten  and  Associates,  is  shown  in  Figs.  7  and  8.  Model  tests 
with  1/100,  1/25,  and  1/10  scale  models  were  used  to  test 
responses  to  waves  and  to  work  out  the  transition  procedures 
for  going  from  horizontal  to  vertical  and  back.  After  testing 
five  different  hull  configurations  with  the  35  ft  sheet  metal 
model  in  El  Capitan  Lake  (near  San  Diego),  we  felt  that  we 
knew  how  to  operate  a  full-sized  FLIP. 

If  the  definition  of  an  expert  is  one  who  has  made  all  the 
mistakes  that  there  are  to  make  (and  doesn’t  repeat  them),  then 
this  is  how  we  acquired  our  expertise  in  flipping.  The  first  hull 
configuration  had  a  D-shaped  cross  section  in  the  reduced 
water  plane  area,  the  rationale  being  to  provide  a  straight  keel 
for  faster  towing.  This  hull  shape  turned  upside-down  during 
flipping  because  the  center  of  buoyancy  in  the  reduced  cross 
section  was  below  the  center  line.  A  spline-shaped  reduced 
area  rolled-over  90°,  and  a  clumsy  configuration  of  40 
controlled  rectangular  tanks  rolled-over  nearly  45°.  Finally, 
we  went  to  a  simple  cylindrical  configuration  very  similar  to 
the  configuration  shown  in  Fig.  6,  with  some  ballast  in  the 
horizontal  keel  to  insure  roil  stability  during  flipping. 

The  only  trouble  remaining  was  plunging.  That  is,  even 
though  the  model  flipped  in  the  right  way,  letting  the  water 
freely  enter  the  tanks  indicated  that  the  full-scale  FLIP  would 
plunge  40  ft  in  going  to  the  vertical.  Attempts  to  throttle  the 
flow  in  one  of  the  model  tanks  (1  ank  3  of  FLIP)  resulted  in  the 
collapsing  of  the  sheet  metal.  Finally,  with  a  rather  battered 
model  (reinforced  by  1  x  1/8  in  steel  straps)  we  convinced 
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Fig.  5.  Scbemttic  showing  decay  of  wave  motion  with  depth  for  differei.. 
frequencies  in  deep  water  and  how  corresponding  pressure  variations  affect 
a  spar  buoy  with  area  reduction  as  shown.  A  large  pressure  on  a  small  area 
at  depth  au  will  exert  a  downward  force  tint,  at  a  particular  frequency,  can 
be  opposed  by  a  smaller  pressure  on  a  larger  area  at  depth  cr2- 
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Fig.  6.  Comparison  of  theory  and  experiment  of  FLIP  heave  amplitude 
response  of  /  =  log  to)  and  measured  FLIP  horizontal 

amplitude  response  g  =  log  ixFurIXI^-). 

ourselves  that  we  knew  how  to  do  the  flipping  operation 
successfully  by  making  Tank  3  a  hard  tank  wi*h  a  reinforced 
flat  section  with  small  holes  between  the  top  and  bottom 
sections  (viewed  in  the  horizontal).  In  this  way,  we  could  take 
on  enough  water  to  start  going  to  th»  vertical  and  therefore 
minimize  the  plunging  by  throttling  the  increase  of  ballast  as 
Tank  3  experienced  the  increasing  differential  pressure  of  up 
to  90  psi. 

All  of  the  testing  of  the  model  was  documented  only  on  16 
mm  film  without  any  technical  paper  backup.  In  going  through 
each  evolution  of  hull  design,  I  would  report  back  and  confer 
with  Dr.  Spiess  before  going  on  to  the  next  configuration.  At 
the  conclusion  of  these  tests,  the  battered  and  rusting  model 
was  left  in  the  vertical  position  at  our  calibration  barge  in  El 
Capitan.  It  eventually  rusted  through  and  sank. 

Having  satisfied  ourselves  with  the  flipping  procedure  and 
the  reduced  area  configuration,  the  next  step  was  to  come  up 
with  plans  to  be  submitted  to  shipyards  for  bidding. 

It  was  around  this  time  that  I  placed  the  name  FLIP  on  the  1/ 
100  scale  model  shown  in  Fig.  9  as  a  noun  and  a  verb.  No 
better  name  for  it  came  forth  by  launching  time  and  we 
ultimately  made  the  acronym  stand  for  FLoating  Instrumenta¬ 
tion  Platform. 
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Fig.  7.  Schematic  of  FLIP'S  hull  and  note 
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Fig  8.  Schematic  showing  deck  arrangements  in  vertical  position  with  about 
500  sq  ft  per  deck. 


The  living  quarters  aboard  FLIP  suffer  from  the  fact  that 
Dr.  Spiess  and  I  originally  envisioned  going  to  sea  aboard  the 
nig  that  towed  FLIP  and  then  transferring  to  it  prior  to  the 
flipping  and  conducting  experiments,  relying  on  C  radons  and 
cots  for  our  amenides.  Our  resident  inspector,  Cdr.  E. 
Bronson,  USN  (Ret.),  at  the  Gunderson  Brothers  Engineering 
Co.  where  FLIP  was  built  in  1962  was  also  FLIP'S  first 
ofncer-rn-charge  and  overcame  our  lack  of  foresight  and 
imagination  regarding  living  quarters,  gradually  evolving  a 
way  for  the  crew  (5-6  in  number)  and  scientific  party  (8 


members  reasonably  comfortable,  more  on  occasion)  to  live 
aboard  FLIP  without  going  sdr  crazy.  Excellent  food  and  the 
ability  to  charcoal  broil  on  the  patio  in  sea-state  6  help 
immensely  in  making  life  tolerable  aboard  FLIP. 

Riding  aboard  FLIP  actually  saved  it  on  its  first  major 
expedition  :n  1963,  when  it  was  to  be  an  island  between 
Hawaii  and  A^ska  for  an  experiment  to  measure  storm¬ 
generated  waves  from  the  South  Pacific,  from  a  series  of 
island  based  stadons  [1 1].  As  it  was  being  towed  to  a  station, 
30  to  35  ft  waves  were  encountered  and  caused  a  rupture  in 
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what  is  the  engine  room  deck  in  the  vertical  position  (due  to  a 
design  constraint  for  which  I  was  responsible).  Bronson 
flipped  to  the  vertical  and  had  the  rupture  welded  shut.  They 
were  close  enough  to  the  station  so  that  they  remained  in  the 
vertical  for  the  duration  of  the  experiment  to  measure  the  small 
millimeter  size  waves  from  the  South  Pacific  storms.  From 
this  same  experiment,  Rudnick  [12]  showed  how  the  deter¬ 
mination  of  wave  direction  could  be  facilitated  by  using  a  spar 
buoy  platform  such  as  FLIP. 

Air-conditioning  and  thermal  and  sound  insulation  were 
added  after  an  operation  in  Hawaii  in  1966  in  which,  at 
breakfast  one  morning,  my  feet,  shod  in  tennis  shoes,  began  to 
feel  as  if  they  were  burning.  The  galley  deck,  directly  above 
die  engine  room  in  which  an  air  cooled  3-71  generator  was 
running,  measured  135*F  with  a  thermometer. 

For  the  electronics  laboratory  we  did  build  a  full-scale 
wooden  mock-up  in  which  we  had  twelve  6-ft-high  standard 
relay  tacks  in  four  bays.  We  flipped  this  mock-up  on  land  to 
make  sure  we’d  be  able  to  load  it  in  the  horizontal  and  work 
with  it  in  the  vertical.  Actually,  this  was  not  very  visionary 
because  it  was  colored  by  the  experience  of  loading  19-in 
racks  aboard  submarines.  The  hatch  diameter  dictated  taking 
them  aboard  one  at  a  time,  which  is  how  we  first  loaded  FLIP. 
The  possibilities  for  wiring  damage,  etc.,  become  nearly 
infinite  in  taking  apart  working  equipment  on  shore  and 
reassembling  it  on-board  FLIP  in  another  orientation.  Dr. 
Rudnick  and  E.  Squier  designed  portable  frames  that  could 
hold  assembled  bays  of  three  6-ft  relay  racks  as  integral  units 
which  could  be  loaded  through  new  hatches  to  accommodate 
large  gear.  These  loading  frames  have  wheels  on  two  sides  and 
are  rotated  before  being  loaded  aboard  FLIP  and  are  bolted 
into  position  on  their  sides  before  going  to  sea. 

The  point  of  the  foregoing  narrative  is  to  emphasize  that  Dr. 
Spiess  and  I  felt  that  we  had  a  satisfactory,  if  not  optimal, 
approach  to  a  stable  manned  platform  that  in  addition  to 
satisfying  our  immediate  acoustic  needs  would  be  useful  for 
other  scientific  projects.  By  not  working  out  the  details 
regarding  living  and  working  aboard  FLIP,  it  made  possible  an 
evolutionary  and  practical  approach  dictated  by  the  types  of 
operations  and  the  habitability  needs. 

We  now  feel  that  there  is  a  large  enough  community  of  users 
and  a  large  base  of  experience  upon  which  we  can  draw  to  do  a 
better  job  of  designing  a  FLIP  II  for  research  in  more 
inclement  conditions  where  increased  capability  and  habitabil¬ 
ity  are  critical.  Improved  habitability  will  also  be  the  key  in 
attracting  other  users  who  have  enjoyed  the  amenities  availa¬ 
ble  on  the  larger  oceanographic  vessels.  For  erample,  the 
chief  scientist  on  FLIP  shares  a  four-bunk  compartment  (with 
about  12  x  8  x  7  ft  dimensions)  next  to  the  galley,  head,  and 
lounge  (a  real  euphemism),  with  only  thin  aluminum  partitions 
separating  the  compartments.  When  I  go  to  sea  on  FLIP  I  am 
constantly  reminded  by  my  own  and  my  vocal  shipmates’ 
experiences  about  the  deficiencies  in  habitability. 

While  there  already  exists  a  growing  community  of  users 
and  prospective  users  of  FLIP  who  have  an  idea  of  what  they 
would  like  to  do  at  sea,  it  is  hoped  this  paper  will  stimulate 
further  interest  and  thinking  in  the  design  and  use  of  this  type 
of  platform. 


HI.  Workshop  Recommendations 

To  get  an  indication  of  the  FLIP  II  requirements  from  those 
who  have  used  FLIP  and  from  prospective  users,  the  Stable 
Platform  Workshop  was  sponsored  t_  the  Marine  Physical 
Laboratory  at  the  end  of  June,  1987.  The  focus  of  the 
workshop  was  to  be  the  consideration  of  a  new,  larger,  and 
more  capable  FLIP  for  future  work.  Besides  recommendations 
for  a  new  FLIP,  the  workshop  came  up  with  several  additional 
recommendations  for  different  platforms  according  to  a 
variety  of  needs  [8], 

For  those  v/ho  require  a  stable  platform  with  a  small  water 
plane  area  and  minimal  cross  section  in  the  water  as  well  as 
air,  a  larger  version  of  FLIP  is  necessary.  Representatives  of 
the  marine -oriented  disciplines  of  biological  optics,  air-sea 
interactions,  marine-air  boundary  layer  studies  near  the 
interface  and  up  to  the  cloud  base,  surface  wave  studies,  air- 
sea  flux  measurements,  physical  oceanography,  acoustics, 
biology,  and  submersible  operations  attended  the  workshop. 
In  situ  biological  sampling  and  analysis  without  the  normal 
large  time  delays  in  making  port  calls  can  be  done  on  a  stable 
platform  such  as  a  larger  FLIP  with  appropriate  laboratory 
facilities.  A  larger,  more  habitable  FLIP  also  makes  it  possible 
for  longer  term  stations. 

A  number  of  very  specific  requirements  were  spelled  out  for 
FLIP  II,  mostly  focusing  on  increased  stability,  greater 
payloads,  more  power,  low  acoustic,  electromagnetic,  and 
environmental  contamination,  greater  endurance,  and  heavy 
weather  operating  capability.  Requirements  unique  to  particu¬ 
lar  experiments  are  shown  in  part,  along  with  requirements 
common  to  many  experiments,  in  Table  II.  Foi  a  more 
complete  treatment,  the  Platform  Workshop  Report  should  be 
consulted. 

None  of  the  requirements  has  been  translated  into  payload 
or  stability  requirements,  let  alone  cost.  The  present  FLIP  is 
limited  to  about  a  12  ton  payload  and  a  moment  of  264  ft  tons. 
What  this  means  is  that  we  cannot  mount  two  heavy-duty 
winches  on  FLIP,  each  weighing  about  6  to  7  tons.  As  a  result 
of  substantial  modifications  to  FLIP  since  it  was  launched,  the 
ballast  in  the  horizontal  keel  has  been  increased  from  70  to  200 
tons  and  we  have  now  reached  our  limits  on  improving  its 
capabilities.  The  addition  of  weight  to  the  superstructure  that 
remains  out  of  the  water  in  the  vertical  has  reduced  its  stability 
in  the  vertical. 

It  will  take  a  great  -.  .'x  of  planning  to  come  up  with  a 
satisfactory  FLIP  u  that  can  meet  a  major  fraction  of  the 
requirements  stemming  from  the  workshop,  'the  possibility  of 
providing  a  deep  launch  and  recovery  capability  for  small 
submersibles— manned,  autonomous,  or  tethered— has  been 
discussed  for  some  time  and  was  also  at  the  workshop.  This 
capability  coaid  be  realistic  for  a  FLIP  II  and  would  mean  that 
such  operations  would  be  less  vulnerable  to  weather.  A  launch 
and  recovery  cradle  that  would  go  through  the  surface  on 
tracks  would  be  required. 

Another  possibility  for  improved  capability  by  reducing 
noise  for  personnel  and  experiments  is  to  enclose  the  generator 
engines  as  has  been  done  for  the  Alliance ,  the  new  NATO 
research  ship  at  the  SACLANT  Undersea  Research  Center  at 
La  Spezia,  Italy.  In  order  to  find  out  if  the  main  engines  are 
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TABLE  II 

REQUIREMENTS  FOR  A  NEW  FLIP 
(FROM  THE  PLATFORM  WORKSHOP  REPORT  [REFERENCE  8J) 

COMMON  REQUIREMENTS 

1.  Greater  rough  water  operating  capability: 

60  ft  waves,  weather  100+  ft  waves. 

2.  Reduced  tilt  response. 

3.  Increased  electrical  power  to  1000  .Kw. 

4.  Multiple  separate  laboratory  .  paces. 

5.  Weather-proofed  platforms  and  winches. 

6.  Heavier  winches  (10  ton  +  )  and  booms  (vertical  and  horizontal) 

above  water  and  booms  below  the  surface  for  simultaneous 
multiple  sensor  and  sampler  deployments. 

7.  Continued  capability  for  mooring,  singly/jointly. 

8.  Improved  habitability. 

9.  Safer  rough  weather  personnel/equipment  transfers. 

10.  Environmentally  controlled  computer  space 

11.  Instrumentation  for  and  access  to  data  for  platform  motion 

(pitch,  roll,  accelerations,  heading,  position,  speed). 

12.  Stable  (minimal  roll  and  pitch),  quiet  platform  with 

minimal  acoustic  contamination  in  .v  iter  and  air 

13.  Fixed  heading  keeping  ability,  drifting  mode  or  in  moor. 


REQUIREMENTS  UNIQUE  TO  EXPERIMENTS 
(Partial  listing ) 

Marine  Air  Boundary  Layer  up  to  Cloud  Base 

1.  Balloon  launching  area  and  hanger  for  tethered  balloon 

(approximately  15  ft  long)  and/or  fcytoon. 

Winch  for  tethered  balloon  or  kytoon. 

2.  Helium  bottle  storage  area. 

3.  Mounting  platform  for  doppler  radar  (dome  up  to  3  m  diameter) 

and  other  remote  sensors  or  accomodations  for  van  housing 
these  instruments. 

Near  Surface  Structure  of  Marine  Air  Boundary  Layer 

1.  Vertical  tower,  reaching  up  to  50  m  above  surface. 

Good  access  to  lower  10  m  without  being  subject  to 
flow  distortion  (as  in  BOMEX  experiment) 

2.  Wind  tunnel  studies  (as  Mollo-Christenscn  did  for  BOMEX). 

3.  Mounting  area,  vertical  window  for  acoustic,  laser  and 

radar  sounder;. 

Air-Sea  Flux  Measurements 

1.  Low  aerodynamic  drag  and  weather  vaning  shape. 

2.  Surface  wave  sensing  capability. 

3.  Low  radio  frequency  contamination  of  radiometric,  hot  film, 

and  other  sensitive  instruments. 

4.  Stability  better  than  present  FLIP  to  permit  eddy  correlation 

measurements  and  other  direct  measurements  such  as  that 
of  insolation  without  correcting  for  platform  motion. 


running,  one  literally  has  to  feel  the  acoustic  enclosure  to  see  if 
it  is  warm. 

IV.  Mooring  Operations  with  FLIP/FLIP  II 

Over  a  period  of  nearly  twenty  years,  MPL  has  developed 
the  capability  for  placing  FLIP  into  tight  three-point  moors  on 
a  routine  basis  in  2500  fathom  water.  The  need  'or  transilion- 
from  single-point  moors  to  three-point  moors  arose  from 
strum  and  flow  noise  contamination  of  acoustic  data  from 
hydrophones  suspended  below  FLIP.  The  mooring  operations 
for  FLIP  have  been  discussed  in  detail  by  Bronson  [3],  [d]. 
With  respect  to  a  larger  FLIP,  one  of  the  desirable  characteris¬ 
tics  would  be  its  ability  to  self-moor  or  to  be  able  to  carry 
enough  line  so  that  either  more  deck  space  would  be  available 
aboard  the  tug  or  that  a  wider  variety  of  ships  could  be  used  for 
towing. 

For  a  tight  three-point  moor,  about  80  tons  of  gear  are 
loaded  aboard  the  tug  that  tows  FLIP  to  its  station.  Each  leg  of 
the  moor  requires  three  reels  of  1.5-in  mooring  line,  each 
600C-ft  long,  a  2000-ft  leader  of  similar  line,  10  tons  of 
surplus  anchor  chain  stored  in  an  anchor  chain  bin,  a  750  lb 
Danforth  anchor  with  augmented  flukes  and  palms,  and  a 
shear  link  assembly  to  permit  recovery  of  the  line  but  not  the 
chain  or  anchor.  As  seen  in  Fig.  10,  where  the  anchor  chain  is 
flaked  out  on  deck,  there  is  not  room  for  much  else.  By  going 
to  anchor  chain  bins  we  leave  some  deck  space  for  other 
activities,  such  a:  towing  sound  sources  and  deploying 
acoustic  navigation  transponders. 

Even  with  the  large  Navy  tugs  there  is  little  .oom  left  after 
all  the  mooring  gear,  etc.,  are  lashed  down  on  the  fantail. 
Commercial  tugs  oi  mud  boats  of  5000  hp  and  at  least  185-fit 
lengths  are  generally  satisfactory.  Currently,  the  possible  use 
of  Kevlar  mooring  lines  is  being  explored  to  cut  down  on  deck 
space  needs 

The  possibility  of  carrying  its  own  mooring  lines  provides 
several  advantages  for  a  larger  FLIP.  Single  mooring  lines 
without  sections  joined  by  shackles  would  speed-up  mooring 
and  recovery  times.  There  would  be  a  substantial  decrease  on 
deck  loading  and  space  requirements  for  the  tow  ship;  about 
30  tons  and  around  640  sq  ft  for  the  ten  reels  of  1 .5-in  nylon 
line  of  a  80  000#  breaking  strength.  These  lines  could  be 
stored  in  separate  tanks  without  requiring  reels.  The  present 
nylon  lines  we  use  are  strong  enough  so  that  they  could  be  used 
to  moor  a  new  FLIP  of  the  dimensions  discussed  here.  The 
only  feature  required  to  accomplish  this  would  be  to  use 
acoustic  command  links  to  release  the  anchor/anchor  chain 
assembly  upon  termination  of  the  moor,  a  requirement  that 
doesn’t  tax  the  state  of  the  art.  This  would  make  it  possible  to 
recover  the  two  downwind  legs  simultaneously  as  long  as  the 
current  and  wind  are  in  the  same  direction  so  that  there  would 
be  no  danger  of  fouling. 

In  general,  we  now  use  either  a  commercial  tugboat  or,  if 
Navy  services  are  available  for  Navy  projects,  we  use  fleet 
tugs,  T-AFT’s  operated  by  the  Military  Sealift  Command. 
The  7200  hp  T-AFT’s  can  tow  FLIP  as  fast  as  10- 1 1  knots  in 
calm  seas;  for  rough  seas  tow  speeds  as  low  as  2  krots  are 
required,  depending  on  the  FLIP  hull  stress  (monitored  on 
line)  and  other  considerations.  A  larger  hull  with  a  50-percent 
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Fig.  9.  P.  Rotate*,  F.  H.  Fisher,  F.  N.  Sjr ess,  designers  of  FLIP,  holding 
thefirsjnoddof  FLIPmoifcfromaLouirvilleilug^rfcst  by  M.  Crouch  of 
MPU 


increase  in  diameter  means  that  stiffness  will  be  increased  by 
over  three  times,  which  in  tarn  means  that  fester  tows  in  rough 
seas  can  be  achieved  without  exceeding  stress  limits  on  the  hull 
and  that  a  more  comfortable  ride  for  the  crew  will  be  possible. 
The  high  L/D  ratio  of  30  to  1  in  the  narrow  section  of  FLIP 
means  that  it  shimmies  when  it  is  slammed  by  waves, 
something  that  would  be  reduced  greatly  witn  a  larger  FLIP. 

At  one  tune,  it  was  planned  to  deploy  the  MPL  residual 
noise  horizontal  array  connected  to  FLIP  using  the  tugboat  to 
launch  the  array  and  the  anchor/crown  line  at  the  bitter  end  of 
the  array  This  proved  to  be  too  difficult  to  do  safely,  even  in 
modest  ssss  and  with  modest  winds  inside  the  Southern 
California  bight,  let  aione  out  in  the  open  ocean.  At  Dr.  V. 
Anderson’s  suggesting  a  joint  mooring  of  FLIP  and  ORB,  it 
now  seems  feasible  and  safe  to  deploy  the  1500-tn-long,  200 
dement  residua!  noise  array  horizontally  at  the  sound  channel 
axis  depth  of  about  750  m.  The  mooting  configuration  is 
shown  schematically  in  Fig.  11.  The  idea  would  be  to  tow 
FLIP  and  ORB  in  tandem  from  the  tug  and  at  the  experiment 
site  to  place  FLIP  into  a  three-point  moor  first.  ORB  world 
pay  out  a  tether  mooring  line  tc  the  downwind  side  of  FLIP 
rad  either  drift  downwind  to  Icctaio'  or  use  thrusters  to 
maneuver  into  position.  At  the  proper  distance  from  FLIP, 
depending  on  how  long  a  horizontal  array  is  to  be  deployed 
and  at  what  depth,  the  tug  would  pick  up  tbs  two  mooring  lines 
from  ORB,  attach  the  anchors  and  chains,  and  complete  the 
five-point  moor.  A  simpler  method  would  be  to  put  FLIP  and 
ORE  into  independent  three-point  moon  if  the  tether  between 
them  proved  troublesome;  for  example,  should  traffic  go 
between  them. 

Once  the  joint  moor  is  completed,  a  wire  or  ether  messenger 
from  FLIP  is  used  to  deploy  the  horizontal  array  stored  on¬ 
board  ORB;  the  array  is  monitored  electrically  front  the 
topside  electronics  C">  ORB  during  deployment  so  that  leaks  or 


other  defects  can  be  corrected.  Weights  at  both  ends  of  the 
array  are  used  to  provide  tension  over  the  array  in  order  to 
keep  it  reasonably  straight  and  to  keep  it  below  the  tether  line 
to  prevent  fouling.  Depth  gages  in  the  array  as  well  as  acoustic 
navigation  will  be  used  to  monitor  the  position  of  the  array. 
Real-time  data  on  element  position  will  be  incorporated  in  the 
array  data  processing  to  permit  beamforming. 

For  this  type  of  operation  a  new  and  larger  FLIP  would 
make  it  possible  to  store,  deploy,  and  record  both  the 
horizontal  and  vertical  arrays  in  one  location  and  analyze  in 
real  time  the  data  from  both  through  the  same  system. 

With  a  larger  FLIP  it  may  be  possible  to  operate  in 
conjunction  with  a  research  ship  tethered  at  varying  distances 
from  it.  To  save  on  fuel,  this  was  actually  done  with  the  R/V 
Horizon  during  the  1963  expedition  for  Dr.  Munk. 

V.  Typical  Operating  Costs 

With  the  support  subsidy  of  FLIP/ORB  from  NAVSEA,  as 
with  the  research  submarine  Alvin ,  the  costs  to  a  user  are 
about  $1300-1400  per  day  for  either  platform.  There  are 
officers-in-charge  for  each  platform  and  four  other  crewmen; 
that  is,  we  have  one  crew  and  two  captains  so  that  only  one 
platform  can  go  to  sea  at  a  time  unless  other  temporary  help  is 
needed  for  simultaneous  deployment  of  both  platforms. 

These  costs  are  in  addition  to  harbor  and  ocean-going  tug 
services,  which  vary  according  to  the  market.  Currently, 
harbor  tugs  are  about  $3500  per  operation  and  ocean  tug 
services  can  be  up  to  $6000  per  day  or  else  be  essentially  free 
if  the  project  is  a  Navy  cask  for  which  Navy  T-ATF’s  are 
available.  If  mooring  is  required,  costs  for  a  three-point  moor 
in  deep  water  is  about  $23  000,  including  the  cost  for  30  tons 
of  anchor  chain,  three  shear  link  assemblies,  amortization  of 
mooring  lines,  and  augmentation  of  flukes  ard  palms  of  three 
750  ib  Danforth  anchors.  Currently,  as  mentioned  above, 
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Fig  10.  Mooring  gear  on  deck  of  tugboat.  About  80  tons  of  chain,  mooring 
line,  etc  ,  arc  necessary  lor  a  three-point  moor  in  deep  water. 


Fig.  11.  Schematic  of  a  jo;nt  mooring  of  FLIP  and  ORB  to  permit  controlled 
deployment  of  horizontal  array  in  addition  to  vertical  arrays. 


about  80  tons  of  deck  gear  must  go  on  the  tug  in  order  to  do  a 
deep  water  three-point  moor.  Costs  for  other  gear  aboard  the 
tow  ship  (winches,  cranes,  sound  sources,  etc.)  would  be  in 
addition  to  the  above. 

A  larger  FLIP  might  require  larger  anchors  and  possibly 
more  chain  in  order  to  withstand  1  knot  currents  and  high 
winds.  The  underwater  cross  section  would  be  increased  by 
about  70  percent  for  the  dimensions  discussed  in  the  next 
section.  Fortunately,  the  anchor  chain  bins  and  anchors  would 


still  not  require  much  deck  space  aboard  the  tug,  even  for  the 
increased  load  for  a  new  FLIP. 

VI.  Current  Planning 

We  are  focusing  our  thinking  on  a  new  FLIP  that  could  meet 
the  requirements  set  forth  in  Table  II  and  which  could  also  be 
berthed  in  the  same  area  used  for  FLIP.  With  a  length  of  420  ft 
and  a  maximum  diameter  of  30  ft,  it  seems  that  most  of  the 
requirements  could  be  met.  For  operating  high-resolution 
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doppler  sonars  with  a  very  narrow  beam.  Dr.  Pinkel  would 
like  to  reduce  the  tilting  motion  of  FLIP  II  to  about  half  the 
motion  we  experience  on  FLIP  in  sea-state  5  and  6. 

We  are  investigating  costs  and  designs  for  FLIP  TI  in 
conjunction  with  Glosten  and  Associates.  As  we  become  more 
"interested  in  low-frequency  acoustic  research  requiring  larger 
arrays  and  in  doing  joint  deployments  of  vertic  J  and  horizon 
tal  arrays,  especially  in  more  inclement  conditions,  a  larger 
platform  becomes  necessary.  Physical  oceanographers  have 
found  the  stability  and  size  of  FLIP  to  be  highly  desirable  and, 
for  many  of  the  same  reasons  as  for  the  acoustic  community, 
find  that  a  larger  platform  is  necessary.  At  this  point, 
therefore,  we  invite  the  interest  of  other  scientists  and 
engineers  who  wish  to  consider  the  use  of  a  larger  version  of 
FLIP  and  their  submitting  comments  and  suggestions  regard¬ 
ing  its  design. 
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